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Abstract 

We provide a rigorous numerical computation method to validate periodic, homoclinic and 
heteroclinic orbits as the continuation of singular limit orbits for the fast-slow system 

fx' = f{x,y,e), 

[F = <^g{x,y,e) 

with one-dimensional slow variable y. Our validation procedure is based on topological tools 
called isolating blocks, cone condition and covering relations. Such tools provide us with 
existence theorems of global orbits which shadow singular orbits in terms of a new concept, 
the covering-exchange. Additional techniques called slow shadowing and m-cones are also 
developed. These techniques give us not only generalized topological verification theorems, but 
also easy implementations for validating trajectories near slow manifolds in a wide range, via 
rigorous numerics. Our procedure is available to validate global orbits not only for sufficiently 
small e > 0 but all e in a given half-open interval (0, eo] ■ Several sample verification examples 
are shown as a demonstration of applicability. 

Keywords: singular perturbations, periodic orbits, connecting orbits, rigorous numerics, 
covering-exchange. 
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1 Introduction 

1.1 Background of problems and our aims 

In this paper, we consider the dynamical system in K" x K of the following form: 

\x' = f{x,y,e), 

\y' = e9(x,y,e), 


( 1 . 1 ) 


where ' = d/dt is the time derivative and f,g are C"’-functions with r > 1. The factor e is 
a nonnegative but sufficiently small real number. We shall write 0) as ( |l.l[ )e if we explicitly 
represent the e-dependence of the system. The system 0 can be reformulated with a change of 
time-scale variable as 

iex = f{x,y,e), 

\y = 5p,y,e), 

where ' = d/dr and r = t/e. One tries to analyze the dynamics of 0, equivalently ( |1.2[ ), by 
suitably combining the dynamics of the layer problem 


x' = f{x,y,0), 

y' = o, 


(1.3) 
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and the dynamics of the reduced problem 


0 = f{x,y,0), 

y = 9ix,y,0), 


(1.4) 


which are the limiting problems for e = 0 on the fast a nd t he slow time scale, respectively. Notice 
that (1.4) makes sense only on f{x,y,0) = 0, while (1.3) makes sense in whole as the y- 

parameter family of a:-systems. The meaning of the “e —>■ 0-limit” is thus different between (1.1) 


and (1.2). This is why (1.1) or (1.2) is a kind of singular perturbation problems. In particular, (1.1) 


or (1.2) is known as fast-slow systems (or slow-fast systems), where x dominates the behavior in 


the fast time scale and y dominates the behavior in the slow time scale. 

When we study the dynamical system of the form we often consider limit systems ( |1.3[ ) 

and (1.4) independently at first. Then ones try to match them in an appropriate way to obtain 
trajectories for the full system One of major methods for completely solving singularly 

perturbed systems like (1.1) is the geometrie singular perturbation theory formulated by Fenichel 
uni, Jones-Kopell m, Szmolyan and many researchers. A series of theories are established 


so that formally constructed singular limit orbits of (1.3) and (1.4) can perturb to true orbits of 


(1.1) for sufficiently small e > 0. In geometric singular perturbation theory, there are mainly two 
key points to consider. One is the description of slow dynamics for sufficiently small e near the 
nullcline {{x,y) \ f{x,y,0) = 0}. The other is the matching of fast and slow dynamics. As for 
the former, Fenichel m provided the Invariant Manifold Theorem for describing the dynamics on 
and around locally invariant manifolds, called slow manifolds, for sufficiently small e > 0. Such 
manifolds can be realized as the perturbation of normally hyperbolic invariant manifolds at e = 0, 


which are often given by submanifolds of nullcline in (1.3). As for the latter, Jones and Kopell 


originally formulated the geometric answer for the matching problem deriving Exchange Lemma. 
This lemma informs that the manifold configuration upon exit from the neighborhood of slow 
manifolds under the assumption of transversal intersection between tracking invariant manifolds 
and the stable manifold of slow manifolds. Afterwards, Exchange Lemma has been extended in 
various directions, e.g. [El [22j [33] . Combining these terminologies, one can prove the existence of 
homoclinic or heteroclinic orbits of invariant sets near singular orbits for sufficiently small e > 0. 
There are also topological ways to prove their existence for sufficiently small e provided by, say. 
Carpenter [1], Gardner-Smoller [E] and so on, by using algebraic-topological concepts such as 
the mapping degree or the Conley index [2l[25|. Such topological approaches also mention the 
existence of periodic orbits near singular orbits. 

On the other hand, all such mathematical results do not give us how large such sufficiently 
small e is. In other words, it remains an open problem whether there exist global orbits given 
by the continuation of singular orbits for a given e. This intrinsic problem has been mentioned 
in many discussions (e.g. dl)- From the viewpoint of numerical computations, if e > 0 is 
sufficiently small, ([TtJ becomes the stiff problem and numerically unstable. Although the effective 
method for computing slow manifolds is provided by Guchenheimer and Kuehn |16j , computations 
for extremely small e (e.g. close to machine epsilon) is still hard to operate correctly. These 
circumstances show that there are gaps between mathematical results (i.e. sufficiently small e) 
and numerical observations (i.e. given e) for completely understanding dynamics of fast-slow 
systems. These are mainly because there are no estimations to measure mathematically rigorous 
consequences not only quantitatively but also qualitatively. The construction of procedures which 
bridge mathematical results and numerical observations is necessary to completely understand 
phenomena in concrete dynamical systems, which is also the case of singular perturbation problems. 

Our main aim in this paper is to provide implementations for validating the continuation of 


various global orbits of (1.1) for all e S [0, eo] rigorously, where eo > 0 is a given number. In other 
words, we provide a method to validate 


(Main 1) the singular limit orbit Hq for (1.1) with e = 0, as well as 
(Main 2) global orbits near Hq for all e € (0, eo] for a given cq > 0. 
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Singular limit orbit means the union of several heteroclinic orbits in (1.31 and the submanifolds of 
the nullcline {f{x,y,0) = 0}. Global orbits mean homoclinic, heteroclinic orbits of invariant sets 
and periodic orbits. 


To this end, we provide the notion of covering-exchange, which is a topological analogue of 
Exchange Lemma. This concept consists of the following topological notions with suitable assump¬ 
tions: (i) isolating blocks, (ii) cone conditions and (iii) covering relations. These three notions have 
been already applied to validations of global orbits in dynamical systems very well, as stated in 
Section [L^ The covering-exchange is constructed by those notions, and informs us 


the existence of slow manifolds for (1.1) as well as normally hyperbolic invariant manifolds 
for (|1.3|), and 


• the existence of trajectories not only which converge to invariant sets on slow manifolds but 
also which exits neighborhoods of slow manifolds after time T = 0(l/e). 


These properties and the general consequence of covering relations yield the existence of global 
orbits. We also generalize the covering-exchange by introducing the additional notion of slow 
shadowing, which guarantees the local existence of trajectories which shadow slow manifolds with 
nonlinear structure. This notion enables us to trace trajectories which not only tend to slow man¬ 
ifolds but also stay near slow manifolds for time 0(l/e) from topological viewpoints. This concept 
is also very compatible with numerical computations, in particular, for validating trajectories near 
slow manifolds. 

The other main tool to establish our procedure is the assistance of rigorous numerics, namely, 
computations of enclosures where mathematically correct objects are contained in the phase space. 
All such computations can be realized by interval arithmetics and mathematical error estimates. 
Combination of the covering-exchange, slow shadowing and rigorous numerics in reasonable pro¬ 
cesses provide us with a method proving (Main 1) and (Main 2) simultaneously. 

Note that there are two approaches to consider singular perturbation problems; one is the 
continuation of structures from the singular limit systems (e = 0) to the full systems (i.e. e > 0), 
and the other is the consideration of full systems to the singular limit e —0. Our attitude is the 
former. In particular, we consider our problems via topological approach on the basis of geometric 
singular perturbation theory. 


This paper is organized as follows. In Section we briefly review Fenichel’s invariant mani¬ 
fold theorem and topological notions called covering relations and isolating blocks. A systematic 
procedure of isolating blocks for validations of invariant manifolds with computer assistance (e.g. 
[321 US]) is also discussed. 

In Section [^ we show how slow manifolds can be validated in given regions with an explicit 
range [0, eo] of e. One sees that our fundamental arguments are basically followed by the proof in 
Jones’ article m- Such arguments can be validated via the construction of isolating blocks and 
singular perturbation problems’ version of cone conditions (cf. [39]) and Lyapunov condition [23]. 

In Section]^ we discuss treatments of slow dynamics. First we introduce the new notion called 
the covering-exchange for describing the behavior of trajectories around slow manifolds (Section 


tracking invariant manifolds near slow manifolds in a suitable sense. This concept also solves 
the matching problem between fast and slow dynamics. We also provide a generalization of the 
covering-exchange; a collection of local behavior near slow manifolds called slow shadowing, drop 
and jump (Section |4.3[ ). These concepts enable us to construct true trajectories in full system which 
shadow ones on slow manifolds in reasonable ways via rigorous numerics. Furthermore, we provide 
a slight extension of cones, called m-cones, which enables us to sharpen enclosures of stable and 
unstable manifolds of normally hyperbolic invariant manifolds (Section |4.4[ ). The main idea itself 
is just a slight modification of cone conditions stated in Section]^ But this technique gives us a 
lot of benefits in many scenes incorporating with rigorous numerics. On the other hand, dynamics 
on slow manifolds should be considered when slow manifolds exhibit the nontrivial dynamics such 
as fixed points, periodic orbits, homoclinic orbits, etc. As an example, we discuss validations of 


4.1). This concept is a topological analogue of Exchange Lemma so that we can reasonably validate 
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nontrivial fixed points on slow manifolds (Section [4^ . In the end of Section]^ we discuss unstable 
manifolds of invariant sets on slow manifolds (Section 4.6). To deal with these manifolds 


we 


discuss the invariant foliations of slow manifolds, and translate this fiber bundle structure into the 
terms of cones and covering relations. This is one of key considerations of heteroclinic orbits in 


([13- 

In Section the existence theorems for periodic and heteroclinic orbits near singular limit 
orbits with an explicit range [ 0 ,eo] of e are presented. 


As demonstrations of our proposing implementations, we study the FitzHugh-Nagumo equation: 

{ u' = V 

v'= 6~^{cv — f{u) + w) (1.5) 

w' = €c~^{u — jw), 


where f{u) = u{u — a)(l — u), a G ( 0 , 1 / 2 ) and c, 7 , 5 > 0 . 

The existence of global orbits such as periodic or homoclinic orbits for sufficiently small e are 
widely discussed by many authors (e.g. mn]). Note that, as we mentioned before, the existence 
of those orbits for a given e remains an open question. Our proposing ideas lead a road to answer 
this question. 


Computer Assisted Result 1.1 (Existence of homoclinic orbits). Consider (1.5) with a = 0.3, 
7 = 10.0 and 6 = 9.0. Then for all c G [0.799,0.801], there exist the following two kinds of 
trajectories: 


1. At € = 0, a singular heteroclinic chain Hq consisting of 

• heteroclinic orbit from the equilibrium po near (0,0,0) G to the equilibrium qi near 
(1,0,0) G 

• heteroclinic orbit from the equilibrium pi near (0.870020061,0,0.06362) G to the 
equilibrium qo near (—0.12966517,0,0.06335) G and 

• two branches of nullcline {{u,v,w) \ v = 0 ,/(m) = w} connecting two heteroclinic orbits. 

2. For all e G (0,5.0 x 10“^], a homoclinic orbit of p’^ « po near Hq. 

The precise statements and other sample validation results are shown in Section]^ Throughout 
the rest of this paper we make the following assumption, which is essential to our whole discussions 
herein. 


Assumption 1.2. Vector fields f and g have the following form: 

m 

f{x,y,e) = fo{x,y) + ^ eV*(a:, p) -b o(e™), fo ^ 0, 

m 

9{x, y, e) = go{x, y)+ '^ Cgz{x, y) + o(e™), po ^ 0. 

i=l 


General dynamical systems depend on parameters. For example, the FitzHugh-Nagumo system 
(1.5) contains a,c,"f,S as parameters. Throughout this paper we do not care about parameter 
dependence of dynamical systems unless otherwise specified. 


1.2 Several preceding works related to global trajectories and singular 
perturbation problems with rigorous numerics 

There are many preceding works for the existence of global orbits with rigorous numerics for 
regular dynamical systems, namely, g = 0 case. For example, Wilczak and Zgliczyhski m apply 
the topological tool called covering relations to the existence of various type of trajectories in 
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dynamical systems, such as periodic orbits, homoclinic orbits and heteroclinic orbits. The essence 
of covering relations is to describe behavior of rectangular-like sets called h-sets and apply the 
mapping degree to the existence of solutions. One of powerful properties of covering relations is 
that every h-sets can be used as a joint of trajectories and that we can validate various complicating 
behavior of dynamical systems. Indeed, for example, [3S] and [55] by Wilczak validate various type 
of complex trajectories such as Shi’lnikov homoclinic solutions, heteroclinic solutions and infinitely 
many periodic solutions in concrete systems (e.g. Michelson system or Rossler system). 

On the other hand, van der Berg et. al. [34] produce the other approach of rigorous numerical 
computations of connecting orbits using radii polynomials and parametrization technique. Their 
main idea is to reduce the original problem to a projected boundary value problem in an infinite 
dimensional functional space via a fixed point argument. Their formulation involves a higher order 
parametrization of invariant manifolds near equilibria for describing stable and unstable manifolds. 
Their approach is free from integrations of vector fields. Hence one can validate various additional 
properties of invariant manifolds without any knowledges of the existence of trajectories [5]. 

As for rigorous numerics for singular perturbation problems, Gameiro et. al. m provide a 
validation method combined with the algebraic-topological singular perturbation analysis. Such 
analysis is knows as the Conley index theory [ani]. They actually apply its singular perturbation 
version mm to the singularly perturbed predetor-prey model with two slow variables. As a 
result, they prove the existence of topological horseshoe, in particular, infinite number of periodic 
orbits with computer assistance. Their results show, however, the existence of solutions for only 
sufhciently small e and hence the bound of e > 0 where the existence result holds is not given. 
When we apply the Conley index technique, it is necessary to provide an appropriate neighborhood 
of desiring orbits whose boundary transversally intersects the vector field defined by the full system 

([g. 

On the other hand, Guckenheimer et. al. [T5] discusses rigorous enclosures of slow manifolds 
with computer assistance within explicit ranges of e. They introduce a concept of computational 
slow manifolds related to slow manifolds in geometric singular perturbation theory and succeed 
validations of slow manifolds in their settings with in explicit ranges of e, while validated ranges of 
e are bounded away from e = 0, say, e G 10“^]. Note that [H] also discusses validations of 

tangential bifurcations of slow manifolds for all e G (0, eo]i where Cq is a given number, say, 10“^. 
One of the other works is a very recent one by Arioli and Koch [T] , which studies the existence and 
the stability of traveling pulse solutions for the (singularly perturbed) FitzHugh-Nagumo system. 
This study, however, focuses only for the parameter range which is not so small, say, e « 0.1 or 
0.001. In other words, the singular perturbation structure is ignored. 

In the case of singular perturbation problems, direct computations of global orbits without any 
ideas are not practical in various scenes both in the rigorous and in the non-rigorous sense. 

Direct applications of such preceding works without any modifications would yield the failure 
of operations if we try to cover not only e which is not so small (e.g. e = 10“® or 10“®) but also 
extremely small e, possibly smaller than machine epsilon. This failure is because either the stiffness 
of problems or the effect of fast dynamics. Even if it succeeds, there would be huge computation 
costs due to very slow behavior around slow manifolds. Of course, there is a matching problem 
connecting fast and slow behavior, which is generally arisen in singular perturbation problems. If 
we can overcome all such difficulties as simple as possible, the scope of applications of preceding 
concepts will dramatically extend. 


2 Preliminaries 


When we consider a fast-slow system (1.11 from the viewpoint of geometric singular perturbation 
theory following Fenichel (e.g. |10]). the central issue is normally hyperbolic invariant manifolds. 
Fenichel’s theory tells us very rich structure of normally hyperbolic invariant manifolds consisting 
of equilibria and their small perturbations. Such structures are fully applied to our arguments. 
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In the beginning of this section, we review several results about normally hyperbolic invariant 
manifolds for fast-slow systems. 

Our main methodologies to consider fast-slow systems are well-known topological tools called 
covering relations^ isolating blocks and cone conditions. These tools well describe behavior of 
solution sets as well as their asymptotic behavior. In successive sections we see that these tools 
work well even for fast-slow systems. In this section, we also review two of such tools. We also 
provide a procedure of isolating blocks suitable for fast-slow systems so that they validate slow 
manifolds, which are available to various systems with computer assistance. Cone conditions for 
singular perturbation problems are stated later. 

Finally note that readers who are familiar with these topics can skip this section. 


2.1 Fenichel’s Invariant Manifold Theorems : review 


Here we briefly review the Fenichel’s Invariant Manifold Theorems (e.g. m), following arguments 
in [TS]. The central goal of these results is the description of flow near the set Sq = {{x,y,0) \ 
f{x,y,0) = 0} with manifold structures: the critical manifold. The critical manifold Sq can be 
considered as the y-parameter family of equilibria of the layer problem (1.3). Under appropriate 
hypotheses, Sq can be represented by the graph of a function x = h{y) for y £ K, where K C 
is a compact, simply connected set. 

A central assumption among Fenichel’s theory is the normal hyperbolicity and the graph rep¬ 
resentation of Sq. 

(F) The set Sq is given by the graph of the C°° function h^{y) for y G K, where the set AT is a com¬ 
pact, simply connected domain whose boundary is an (Z — 1)-dimensional C°° submanifold. 
Moreover, assume that Sq is normally hyperbolic. 

Remark 2.1 (cf. [H]). Recall that the manifold Sq is normally hyperbolic if the linearization of 


(1.3) at each point in Sq has exactly I eigenvalues on the imaginary axis. 

A set M is said to be locally invariant under the flow generated by if it has a neighborhood 
U of M so that no trajectory can leave M without also leaving V. In other words, it is locally 
invariant if for all x € M, then (^([0,t],a:) C V implies that :^([0,t],a:) C M, similarly with [0,t] 
replaced by [t, 0 ] when t < 0 , where is a flow. 

Without the loss of generality, we can assume that h^(y) = 0 for all y G K. Then, there exist 
the stable and unstable eigenspaces, S{y) and U{y), such that dimS'(y) = s and diinU{y) = u 
hold for all y G K. With this in mind, we take the transformation x{g K") (a, 5) G ]R“+® so 

that ( 1 . 1 ) is expressed by 

a' = A(y)a + Fi(x,y,e) 

b' = B(y)b + F 2 (x,y,e) . (2.1) 

y' = eg{x,y,e) 

Here A{y) denotes the uxu matrix which all eigenvalues have positive real part and B{y) denotes 
the s X s matrix which all eigenvalues have negative real part. Fi and F 2 denotes the higher order 
term which admit a positive number 7 > 0 satisfying 


|Ti| < 7(k| + e) as \x\ -)■ 0. 

More precise assumptions for A{y) and B{y) are as follows: there exist the quantities > 0 and 
fiB < d such that 

< ReA for all A G Spec(A( 2 /)) and y G K, (2.2) 

Pb > ReA for all A G Spec{B{y)) and y G K. (2.3) 

The key consideration of slow manifolds is the following Fenichel’s invariant manifold theorem 
in terms of graph representations. 
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Proposition 2.2 (Persistence of invariant manifolds, cf. [IH])- Under the assumption (F), for 
sufficiently small e > 0, there is a function x = h^{y) defined on K. The graph Se = {(a;, y) \ x = 
h'^iy)} is locally invariant under 0 . Moreover h'^ is , for any r < -|-oo, jointly in y and e. 

Proposition |2.2| is just a consequence of the second Invariant Manifold Theorem as follows. 

Proposition 2.3 (Stable and unstable manifolds, cf. [H])- Under the assumption (F), for suffi¬ 
ciently small e > 0, then for some A > 0, 


1. there is a function a = hs{b,y,e) defined on {{b,y,e) | |&| < A, y G K}, such that the graph 
bF®(S'e) = {{a, b,y) \ a = hfib, y, e)} is locally invariant under 0 . Moreover, a = hs{b,y,e) 
is , for any r < -|-oo, jointly in y and e. 

2. there is a function b = /i„(s,y, e) defined on {(a, y, e) | |a| < A,y G K}, such that the 
graph VF“(S'e) = {{a,b,y) \ b = h„(a,y, e)} is locally invariant under 0. Moreover, 
b = hu{a, y, e) is , for any r < -foo, jointly in y and e. 


Fenichel’s theorems, Propositions |2.2| and |2.3| insist that normally hyperbolic invariant man¬ 
ifolds as well as their stable and unstable manifolds persist to locally invariant manifolds in the 
full system (1.1) for sufficiently small e > 0. In other words, slow manifolds can be realized by the 
e-continuation of normally hyperbolic critical manifolds in the layer problem (1.3). The perturbed 
manifold Se for e > 0 is called a slow manifold. One of strategies for constructing such manifolds 
is the construction of a family of isolating blocks and moving cones 1[18)1. Isolating blocks describe 
the behavior of vector fields which are transversal to their boundaries, which are well discussed in 


the Conley index theory [THIS]. Isolating blocks are reviewed in Section 2.3 


Fenichel’s invariant manifold theory informs us not only the existence of perturbed slow mani¬ 
folds but also invariant foliations of lF“(5'e) and W'^{Sf). More precisely, the following result (the 
third Invariant Manifold Theorem) is known. 

Proposition 2.4 (Fenichel fibering, cf. Theorem 6, 7 in [H])- Under the assumption (F), for 
sufficiently small e > 0 the following statements hold. For each v = Ve = (y, e) G S^, 

1. there is a function {a,y) = h1(b) for | 6 | < A sufficiently small so that the graph 

Wfiv) = {{a,b,y,e) \ {a,y) = K{b)} cWfiS,) 

forms a locally invariant manifold in the sense that (p^ pf(t,W^{v)) C W’^[g)fit,v)) holds if 
ipfis,v) G N for all s G [0,t]. Here the set (p,:^N{t,A) denotes the forward evolution of a set 
A restricted to N given by 

A) = {<Pt{t, u)\u & A and (/^^([O, t],u) C N}. 

Moreover, {a,y) = h1(b) is C’’ in v and e jointly for any r < oo. 

2. there is a function {b,y) = h^{a) for | 6 | < A sufficiently small so that the graph 

lT“(z;) = {(a, 6 ,y,e) | {a,y) = hl{b)} CW'^{S,) 

forms a locally invariant manifold in the sense that (pe,N(t,W'^{v)) C W‘^{ipe{t,v)) holds if 
ipe{s,v) G N for all s G [t, 0]. Moreover, (b,y) = h'!j(a) is in v and e jointly for any 
r < oo. 

This invariant foliation is sometimes referred to as Feniciel fibering. This fibering ensures us 
the following representations: 

W'^iA) = y lF“(u), WfiA) = y Wfiv), 

vGA vGA 

where A is a subset of slow manifolds. 
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2.2 Covering relations : review 

Our main approach to tracking solution orbits is a topological tool called covering relations. Cover¬ 
ing relations describe topological transversality of rectangular-like domains called h-sets relative to 
continuous map, and there are various studies not only from the mathematical viewpoint but also 
for applications with rigorous numerics (e.g. [31I351I3Z1IM1I1I])- In the present study, we apply 
this topological methodology to singular perturbation problems. In this section, we summarize 
notions among the theory of covering relations. For a given norm on M™, let i?rn(c, r) be the open 
ball of radius r centered at c S K™. For simplicity, also let Bm = Rm(0,1). We set = {0}, 
Bo(0,r) = {0} and dBo(0,r) = 0. 

Definition 2.5 (h-set, cf. [39l|4T]). An h-set consists of the following set, integers and a map: 

• A compact subset N c K™. 

• Nonnegative integers u(JV) and s(]V) such that u(N) + s{N) = n with n < m. 

• A homeomorphism cat : K" —>■ x satisfying 

Cat (A') = Bu(n) X Bs(^Ny 


Finally define the dimension of an h-set N by dim N := n. 

We shall write an h-set {N,u{N),s{N),cn) simply by N if no confusion arises. Let 

Ac := i?„(Ar) X i3s(Af)j 

Nc '■= dBy_(^jsr'i X i?s(A/), 

Ac . F?|^^Ar) X dBg(^j^y 

N- ■.= c],\N-), N+:=c],\N+) 

The following notion describes the topological transversality between two h-sets relative to 
continuous maps. 

Definition 2.6 (Covering relations, cf. [SHI SI])- Let N,M C R"* be h-sets with u{N) -\- 
s{N),u{M) -\- s{M) < m and u{N) = u{M) = u. f : N ^ j^dimM (;ienotes a continuous mapping 

and /c := Cm ° f o : Ac —>■ R“ x We say A f-covers M (A M) if the following 

statements hold: 

1. There exists a continuous homotopy h : [0,1] x Ac —>■ R“ x R®^^) satisfying 

ho = /c, 

/r([0, l],A-)n Me = 0, 

/r([O,l],Ac)nM+ = 0, 

where hx = h{X, •) (A G [0,1]). 

2. There exists a mapping A : R“ —>■ R“ such that 

( hi{p,q) = (A(p),0), 

A(5R,(0,1))cR“\R;(0,1), (2.4) 

[deg(A,R:, 0 ) ^0 

holds for p G Bu{0, 1), q G i?s(0,1). 

Remark 2.7. In definition of covering relation between A and M, the disagreement of dim A and 
dimM is not essential. On the contrary, the equality u{N) = u{M) = u is essential because the 
mapping degree of it-dimensional mapping A should be derived. 
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The following propositions gives us useful sufficient conditions for detecting covering relations 
in practical situations. 

Proposition 2.8 (Finding covering relations, Theorem 15 in IS]). Let N,M be two h-sets in 
K” such that u{N) = u{M) = u and s{N) = s{M) = s. Let f : N ^ K” be continuous. Let 
fc = Cm o / o : Nc —>■ K“ x Assume that there exists qo € such that following conditions 
are satisfied: 

1. Setting S{M)f = {{p,q) G ]R“ x M® | ||p|| > 1}, 

fc{Bu X {qo}) C int(S'(M);: U M^), 

UN-) = 

n M+ = 0. 

2. Define a map Aq^ : —)• by 


^qoip) ■= ^u{fc{p,qo)), 

where 7r„ : ]R“ x M® —>■ ]R“ be the orthogonal projection onto K“, 7r„(p, q) = p. Assume that 
Aqq{dBu) Bu, deg{Aqq,Bu,0) ^ 0. 


Then N d^M. 

Proposition 2.9 (Covering relation in the case u = 1, Definition 10 in [10]). Let N,M be h-sets 
with u{N) = u{M) = 1. Let f ■. N ^ continuous. Set 

:= {-1} X B;^), := {+1} x Bf^), 

S{N)^ := (-00, -1) X K®W, S{N)^ := (-tl, +oo) x IR®^. 

Assume that there exists qo G i3s(Ar) such that 

f{cN{[-l, 1] X { 90 })) C int(.S(M)^ U M U S{M)^), 
f{N)nM+ = a 

and either of the following conditions holds: 

f{N^) C S{M)^ and f{N^) C S'(M)^, 
f{N^) c S{M)^ and f{N^) c S{M)^. 


Then N M. 

We also consider covering relations with respect to the inverse of continuous maps. 

Definition 2.10 (Back-covering relation, Dehnition 3 and 4 in [S])- Let N be an /i-set. Define 
the h-set N'^ as follows: 

• The compact subset of the quadruple N'^ is N itself. 

• u{N^) = s{N), s{N^) = u{N). 

• The homeomorphism cj^t : M" —>• K" = x is given by 

CN-rix) = jiCN{x)), 

where j : x IR®(^) —>• ]R®(^) x is given by j{p,q) = {q,p). 
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Notice that = N and N’^’ = N^. 

Next, let N,M be h-sets such that u{N) = u{M). Let g : Q G K" —)■ K". Assume that 
g~^ : M —K" is well-defined and continuous. Then we say N g-back-covers M {N M) if 
mt nt 

A fundamental result in the theory of covering relations is the following proposition. 

Proposition 2.11 (Theorem 4 in [H]). Let Ni, i = 0,1, • • • ,k be h-sets such that u{Ni) = u for 
i = 0,1, ■ ■ ■ ,k and let fi : Ni ^ i = 0,1, • • • , fc — 1, be continuous. Assume that, for 

all i = 0,1, - ■ ■ ,k — 1, either of the following holds: 

TV, ^ TV,+i 
or 

N,cD{f-^) and TV, iV,+i. 

Then there is a point p € intTVo such that 

fi o /,_i o • • • /o(p) e intTVi for all i = 0, - ■ ■ ,k-l. 

If we additionally assume Nk Nq, then the point p G intAo can be chosen so that 

fk o fk-i o • • • fo{p) = P- 

We sometimes consider an infinite sequence of covering relations. To deal with such a situation, 
we define the following concept. 

Definition 2.12 (Admissibility, cf. Definition 2.4 in [3S]). Let be h-sets in K™ and 

/ : [J^i Mi —>• K"* be a continuous map. We say the index sequence C {I,-- - ,k}^ is 

admissible with respect to f if Mi. Mi.^.^^ holds of all j G Z. Similarly, we say the index sequence 

j; 

{ij}j^i, C {1, • • • , fc}^ is back-admissible with respect to f if Mi. <J= Mi-^^ holds of all j G Z. In 
this case, f~^ is assumed to be well-defined in a neighborhood of UI^i ^-nd continuous. 

Recall that the second Invariant Manifold Theorem, Proposition |2.3[ claims that the stable 
and unstable manifolds of normally hyperbolic invariant manifolds can be described by graphs of 
smooth functions. The concept of horizontal and vertical disks are useful to describe asymptotic 
trajectories in terms of covering relations for describing these situations. 

Definition 2.13 (Horizontal and vertical disk, e.g. [3Sl[Sn]). Let N be an h-set. Let bs : Rs(Af) -a- 
N be continuous and let {bs)c = cn o bg. We say that bg is a vertical disk in N if there exists a 
homotopy h : [0,1] x Rs(Af) —>■ A^c such that 

ho = {bg)c, 

hi{x) = (0,a;), for all x G Rs(Af)i 
h{t,x) G N^, for all t G [0,1] and x G 9Rs(jv)- 

Let -G N he continuous and let {bu)c = cat o We say that is a horizontal disk in 

N if there exists a homotopy h : [0,1] x B^^n) A"c such that 

Lq — (f^,z)c, 

hi{x) = (x, 0), for all x G 
h(t,x) G N~, for all t G [0,1] and x G dBu(N)- 

Combining these concepts with covering relations, we obtain the following result, which is often 
applied to the existence of homoclinic and heteroclinic orbits. 
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Proposition 2.14 (Theorem 3 in [57], Theorem 3.9 in |5S|). Let Ni, i = 0,1, • • • ,k be h-sets such 

that u{Ni) = u for i = 0 , 1 , • • • ,k and let fi : Ni ^ i = 0 , 1 , • • • , fc — 1 , be continuous. 

Let b : > Nq be a horizontal disk in Nq and let v : Rs(Affc) Nk be a vertical disk in Nk- If 

Ni holds for i = 0, 1 , • • • ,k — 1, then there exists r € such that 

ifi o fi-i O • • • o /o)(6(r)) G Ni+i, for i = 0,1, - ■ ■ ,k-2, 

ifk-i o fk -2 o • • • o /o)(&(t)) G v{Bs(Nk))- 

2.3 Isolating blocks : review and applications to fast-slow systems 

A concept of isolating blocks are typically discussed in the Conley index theory (e.g. 13121]), which 
studies the structure of isolated invariant sets from the algebraic-topological viewpoint. Central 
notions are isolating neighborhoods or index pairs in the Conley index theory, but we concentrate 
our attentions on isolating blocks defined as follows. In our case, the blocks can be considered very 
flexible from the viewpoint of not only covering relations but also rigorous numerics. Moreover, 
isolating blocks play central roles for the existence of slow manifolds (cf. [15] and Section . In 
this section, we firstly review the definition of isolating blocks. Secondly, we show a procedure 
of isolating blocks around equilibria and their alternatives for fast-slow systems with computer 
assistance, following [42j . 

2.3.1 Definition 

Definition 2.15 (Isolating block). Let N C K™ be a compact set. We say N an isolating 
neighborhood if Inv(A^) C int(A^) holds, where 

Inv(iV) := {x G N \ 93(K,x) C N} 

for a flow : M X K™ —>■ on K”^. Next let B C K™ be a compact set and x G dB. We say x an 

exit {resp. entrance) point of B, if for every solution cr : [—(5i,d2] through x = cr(0), with 

i5i > 0 and (52 > 0 there are 0 < ei < i5i and 0 < €2 < 62 such that for 0 < t < 62, 

a{t) ^ B (resp. a{t) G mt{B)), 


and for —ei < t < 0 , 

a{t) ^ dB (resp. a{t) ^ B) 

hold. (resp. i3®“*) denote the set of all exit (resp. entrance) points of the closed set B. We 

call and R®"* the exit and the entrance of B, respectively. Finally B is called an isolating 

block if dB = 5®==** U R®"* holds and R®’^‘‘ is closed in dB. 

Obviously, an isolating block is also an isolating neighborhood. 


2.3.2 Construction around equilibria via rigorous numerics : a basic form 

There is a preceding work for the systematic construction of isolating blocks around equilibria 


|42| . Here we briefly review the method therein keeping the fast-slow system (1.1) in mind. The 
first part is the review of the preceding work [15]. As the second part, we discuss the analogue of 
arguments to fast-slow systems. One will see that such procedures are very suitable for analyzing 
dynamics around invariant manifolds. 

Let AT C be a compact, connected and simply connected set. Consider first the differential 
equation of the following abstract form: 


x' = f{x,X), xG 


A G AT, / : K” X AT 


(2.5) 


which corresponds to the layer problem (1.3). For simplicity, assume that / is C°°. Our purpose 


here is to construct an isolating block which contains an equilibrium of (2.5). 
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Let po be a numerical equilibrium of (2.5) at Aq S RT and rewrite (2.5) as a series around 
(Po,Ao): 

x' = fx{po,Xo)ix - Po) + (2.6) 


where fxiPoj^o) is the Frechet differential of / with respect to x-variable at (po,Ao). f{x,\) 
denotes the higher order term of / with 0(|a;—poP + |^~Ao|). This term may in general contain 
an additional term arising from the numerical error f{po, Aq) ~ 0 . 

Here assume that the n x n-matrix fx(j>o, Aq) is nonsingular. Diagonalizing fx{po, Aq), which is 
generically possible, (2.6) is further rewritten by the following perturbed diagonal system around 

z'j = PjZj + fj{z,X), j = (2.7) 


Here pj e C, z = {zi, - ■ ■ , Zn) and f(z, A) = (fi(z, A), • • • , fniz, A))^ are defined hy x = Pz + po 
and f{x, A) = P(/(z, A)), where P = {Pij)ij=i^... is a nonsingular matrix diagonalizing fx{po, Aq) 
and is the transpose. 

Let iV C K" be a compact set containing pq. Assume that each fj{z,X) has a bound 
in N X K, namely, 

|/j (z. A) \ X = Pz +po e N,X G k'^ C [S~,S^]. 

Then z' must satisfy 





< Zj < Pj 



Vz with x = Pz+pqGN, VA G at. 


For simplicity we assume that each pj is real. We then obtain the candidate of an isolating block 


B in z-coordinate given by the following: 





n 


r 5+ 

5~] 



1=1 

Bj = [Zj , z+] := 

3 

Pj' 

3 

P-3 _ 

if Pj > 0, 

(2.8) 



\ 5- 

5+' 




Bj = [zj,z^] := 

_ 3_ 

5 

_ 3_ 

if Pj < 0. 

(2.9) 



Pj 

Pj _ 


Remark 2.16. In the case that 

Pj is complex-valued for some i 

, fx{Po,Xo) contains the complex 


conjugate of pj as the other eigenvalue. Without the loss of generality, we may assume pj = 
aj + Pj+i = Pj = aj — y/—iPj, I3j ^ 0. To be simplified, we further assume that pj and 

Pj+i are the only complex pair of eigenvalues of fx{po,Xo). The general case can be handled in 
the same manner. The dynamics for Zj and zj+i is formally written by 

z' = pjZj + fj(z,X), 

Zj+i = Pj+iZj+i -b fj+i{z, A). 

Now we would like to consider real dynamical systems. To do this we transform the above form 
into 

w'j = ajWj + (3jWj+i + fj{w, A), 

= -PjWj + ajWj+i + fj+iiw, A) 

via Q = -V^)’ ^ QiVj^yj+iV and (f^Jj+i)'^ = Q(/j,/j+i)^, where 

w = {wi,--- ,w„) is the new coordinate satisfying Wi = yi for i ^ j,j + 1. Let rj{w,X) := 

\//l {w, A)^ + fj+i{w, A)^ and assume that rj{w, A) is bounded by a positive number rj uniformly 
on N X K. Our aim here is to construct a candidate of isolating block and hence we assume that 
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• the scalar product of the vector field and the coordinate vector 


{wj,Wj+i) ■ (ajWj + l3jWj+i + fj{w,X),-l3jWj + ajWj+i + fj+i{w,X)) 


has the identical sign and the above function never attain 0 on | (wj, Wj+i) \ } 

{h, > 0 ). 

With this assumption in mind, we set the candidate of isolating block in i- and i + 1-th coordinate 

Bjj+i ■■= |(wj,wj+i) I 

Its boundary becomes exit if aj > 0 and entrance if Oj < 0. Finally, replace Bj x -Bj+i in the 
definition of B ((2.8) and (2.9|) by Bj^+i. 

A series of estimates for error terms involves N and it only makes sense if it is self-consistent, 
namely, {po} + PB C N. If it is the case, then B is desiring isolating block for (2.5). Indeed, if 
/ij > 0 , then 

2 ^L,.=z -<0 and z '|^,^^+>0 

hold. Namely, the set {z & B \ Zj = z^} is contained in the exit. Similarly if pj < 0 then 




and 


\z,=z+< 0 


hold. Namely, the set {z G B \ Zj = Zj} is contained in the entrance. Obviously dB is the union 
of the closed exit and the entrance, which shows that B is an isolating block. 

Once such an isolating block B is constructed, one obtains an equilibrium in B. 

Proposition 2.17 (cf. [H]). Let B be an isolating block constructed as above. Then B contains 
an equilibrium of (2.5) for all A € AT. 

This proposition is the consequence of general theory of the Conley index ([23]). Note that 
the construction of isolating blocks stated in Proposition 2.17 around points which are not nec¬ 
essarily equilibria implies the existence of rigorous equilibria inside blocks. With an additional 
property such as uniqueness or hyperbolicity of equilibria, this procedure will provide the smooth 
A-parameter family of equilibria, which is stated in Section]^ 

When we apply these ideas to the fast-slow system we only consider the fast system 

x' = f{x,y,e). Let K C M} he as above, eq > 0 and a; = po be a numerical zero of f{x,yo,0) at 
y = yo ^ K■ In this case we set e = 0 for computing numerical zeros. Via the procedure in the 
above, the fast system x' = f{x, y, e) can be generically formulated by the following form; 


a' = Aa-h Fi{a,b,y,e), b' = R 6 -|-^ 2 ( 0 , 6 , p, e). 


( 2 . 10 ) 


Here A is a rt x u-dimensional diagonal matrix each of whose eigenvalues has positive real part. 
Similarly, R is an s x s-dimensional {u -\- s = n) diagonal matrix each of whose eigenvalues has 
negative real part. Fi and F 2 are higher order terms depending onpo and po- Equivalently, writing 
( 2 . 10 ) component-wise. 


a) =/i“aj-hFij(a;,p,e), Rep“ > 0, j = l, 
bb = ^i]bj + F 2 j{x,y,e), Re/iJ < 0, j = l,---,s. 


Let V C M" be a compact set containing pq. As before, assume that every and is real and 
that each Fij., i = 1,2, ji = 1 , • • • ,u, j2 = 1 , • • • ,s, admits the following enclosure with respect 
to IV X AT X [0, Co]: 

{Pt,jiix,y,e) I X = Py+Po e N,y e AT, e S [0,eo]} C (2.11) 
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Define the set C by the following: 


. .a ’ , ,a 

, [&-, 6 +]:= 

..b ’ ..b 

[ 1^3 1^3 \ 


[ h ^^3\ 


Dc =X{ [a^ > X n ’ 4 ] •= 

i=i t=i 

( 2 . 12 ) 

If {(poiyo)} + PDc C N X K holds with an appropriate linear transform P on then this 

procedure is self-consistent. Under this self-consistence, we immediately know that 


a' > 0 V(a, b, y, e) £ Dc x [0, eq] with aj = a'^, 

a'j < 0 V(a, b, y, e) £ DcX [0, eq] with aj = aj, 

b'j < 0 V(a, b, y, e) £ DcX [0, eq] with bj = b^, 

b'j > 0 V(a, b, y, e) £ DcX [0, eq] with bj = bj. 


Remark that these inequalities hold for all e G [0, eq]. This observation is the key point of the 
construction not only of limiting critical manifolds but of slow manifolds for e £ (0,eo], which is 
stated in Section [S] 


Definition 2.18 (Fast-saddle-type block). Let Dc C be constructed by (2.12 1 . Assume that 
K — Bi C M.K We say Dc, equivalently D := c^{D) via a homeomorphism cd, a fast-saddle-type 
block. Moreover, set 


DI~ := {{a,b,y) £ Dc \ Oj = af, j = !,■■■ ,u}, 

Df’+ := {ia,b,y) £ Dc \ bj =bf, j = 1, • • • ,s}, 

Dc ■= {ia,b,y) £ Dc\y £ dBi}, 

Df’- := c-^\Df’-). := c^i(Df-+), D* := c^i(D^). 


We say Df^ (equivalently D-f^ ) the fast-exit of D and D[^'^ (equivalently D-f’^) the fast-entrance 
ofD. 


Remark 2.19. Obviously D is an h-set, but the integer u{D) and s{D) in Definition 2.5 are not 
necessarily equal to u and s, respectively. We do not have any assumptions for D 
definition. Indeed, D is not necessarily an isolating block in the sense of Definition 2.15 
why we omit the word “isolating” from the definition of D. 


in the above 
That is 


This construction can be slightly extended as follows. Let u 32=1 ••• s a sequence 

of positive numbers. Defining 


4:=n 

i=i 


a, ,a+]x 


Is- ■= 


1=1 

<57 


h-j h-j 



^2,32 





(2.13) 


we can prove that Dc is also a fast-saddle-type block if PDc -£ {(po,i/o)} C N x K holds. We 
further know 


o' > 0 y{a,b,y,e) £ Dc x [0,eo] with Oj £ [a+,a+], 

a'j < 0 V(a, 6 , y, e) £ Dc x [0,eo] with Oj £ [a~,aj], 

b'j < 0 V(a, b, y, e) £ Dc x [0, eq] with bj £ [b^,b^] 

b'j >0 V(a, b, y, e) £ Dc x [0, Eq] with b^ £ [bj,b~]. 

This extension leads to the explicit lower bound estimate of distance between £)7A and slow 
manifolds, which is stated in Section 
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2.3.3 Construction around equilibria via rigorous numerics 
approach 


the predictor-corrector 


Here we provide another approach for validating fast-saddle-type blocks. In previous subsection, 
fast-saddle type blocks are constructed centered at {(x,y) | y € K}, where iF C KMs a small 
compact neighborhood oi y G K. All transformations concerning eigenpairs are done at a point 
{x,y). Fig. [i}(a) briefly shows this situation. On the other hand, we can reselect the center of 
the candidate of blocks so that blocks can be chosen smaller. As continuations of equilibria with 
respect to parameters, the predictor-corrector approach is one of effective approaches. We now 
revisit the construction of fast-saddle type blocks with the predictor-corrector approach. 

Let {x,y) be a (numerical) equilibrium for (1.3), i.e., f{x,y,Q) « 0, such that fx{x,y,0) is 
invertible. Let AT be a compact neighborhood of y. The central idea is to choose the center as 
follows instead of {x,y)'. 


( dx 'N 

( s + - y),y] = {x- fx{x,y)~^fy{x,y){y - y),y) , (2.14) 

where x = x{y) is the parametrization of x with respect to y such that x = x{y) and that 
f{x{y),y,0) = 0, which is actually realized in a small neighborhood of y in K* since fx{x,y) is 
invertible. See Fig. [^(b). Obviously, the identification in (2.14) makes sense thanks to the Implicit 
Function Theorem. 




Figure 1: Fast-saddle-type blocks with (a) basic form, and (b) predictor-corrector form. 

(a): The center point x of blocks in x-coordinate is fixed. The blue rectangle shows a fast-saddle- 
type block at y G K. This procedure can be realized just by extending error bounds in (2.11) 


to y-directions. Flowever, the “higher order” term Fij^ (x, y, e) may contain the linear term with 
respect to y. This may cause the increase of error bounds. 

(b): The center point of blocks is moved along the tangent line (or plane for higher dimensional 
systems) at {x,y). The blue rectangle shows a fast-saddle-type block atyGK. In principle, size of 
blocks will be smaller than (a) since the linear term with respect to y in error term is approximately 
dropped. 

Around the new center, we define the new affine transformation T : (z, w) i—)■ {x, y) as 

{x, y) = T{z, w) := {Pz-\-x- fx^x, y)~^ fy{x, y)w, w-\-y) . (2.15) 

where P is a nonsingular matrix diagonalizing Dx{x,y). Over the new (z,ri;)-coordinate, the fast 
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system x' = /(x, y, e) is transformed into the following: 
z' = P~^ (x' -f /a; 

= P~^ (^f(x,y,e') + ef^ ^fyg(x,y,e)J 

= P~^ (^fx{Pz - fx ^fyw) + f{z,w,e) + efx ^fygix,y,e)^ 

= Az + P~'^ (^-fyW + f{z,w,e) + efx ^ fyg{Pz + x- fx ^ fyW,w + 
= Az + F{z, w, e), 


where fx = fx{x,y) and Dyf = fy{x,y), and A = diag(/rj,--- ,g°). The function 

f{z, w, e) denotes the higher order term of / with 0(|2p -|- |r(;|). Dividing z into (a, b) corresponding 
to eigenvalues with positive real parts and negative real parts, respectively, as in (2.10), we can 
construct a candidate set of fast-saddle-type blocks as in (2.12). By the similar implementations 
to Section 2.3.2 we can construct fast-saddle-type blocks centered at (2.14) for y £ K. 

Note that the higher order term f{z, w, e) contains the linear term of w as fyW with small errors 
in a sufficiently small neighborhood K of y. This fact indicates that, in principle, size of blocks 
becomes smaller than those in Section 2.3.2 as seen in Fig. [i}(b). This benefit is also useful for 
the following arguments. 


3 Slow manifold validations 


In this section, we provide a verification theorem of slow manifolds as well as their stable and 
unstable manifolds. Our goal here is to provide sufficient conditions to validate not only the 
critical manifold Sq but also the perturbed slow manifold of (l.l)e for all e € (0, cq] in given 
regions. 

Recall that Fenichel’s results (Propositions 2.2 2.3) assume normal hyperbolicity and graph 
representation of the critical manifold Sq. These assumptions are nontrivial, but very essential to 
prove the persistence. Our verification theorem contains verification of normal hyperbolicity and 
graph representation of Sq. 

The main idea is based on discussions in |18j . For technical reasons, we use a multiple of e as 
the new auxiliary variable. We set e = ga and cr := cq > 0, where eo is a given positive number. 
We add the equation g' = 0 to (2.1). Furthermore, we consider the following system instead of 
(2.1) for simplicity: 


a' = Aa + Pi(x, y, e) 
b' = Bb + F2(x,y,€) 
y' = eg{x,y,e) 

.v' = o 


(3.1) 


Here A denotes the u x u matrix which all eigenvalues have positive real part and B denotes the 
sxs matrix which all eigenvalues have negative real part. Note that matrices A and B are (locally) 
independent of y. This formulation is natural when we the construction of fast-saddle-type blocks 
stated in Section [2)3] is taken into account. 


Let M be a fast-saddle type block for (1.1). Section 2.3 implies that the coordinate represen¬ 
tation, Me, is given by (2.12) (or (2.13)), which is directly obtained from the system (3.1). A fast- 
saddle-type block M has the form (2.12), which has a-coordinate, 5-coordinate and y-coordinate 
following (3.11. With this in mind, we put following notations. 


Notations 3.1. Let tTo, tt;,, tTj,, iTa.b, T^a,y and be the projection onto the a-, 5-, y-, (a, 5)-, 
(a, y)- and (5, y)-coordinate in M, respectively. If no confusion arises, we drop the phrase “in M” 
in their notations. 

We identify nonlinear terms Fi(x,y,e), F 2 (x,y,e) and g{x,y,e) with Fi{a,b,y,e), ^ 2 ( 0 , 5, y,e) 
and y(a, 5, y, e), respectively, via an affine transform x(g K") >->■ (a, 5) S K“+®. 
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For a squared matrix A with Spec (A) C {A € C | ReA > 0}, Aa > 0 denotes a positive number 
such that 

Ayi < ReA, VA S Spec(Al). (3.2) 

Similarly, for a squared matrix B with Spec(i?) C {A G C | ReA < 0}, fj,B < 0 denotes a negative 
number such that 

/i_B > Re/i, V/i G Spec(i?). (3.3) 

Finally, let dist(-,-) be the distance between A and B, A, B C given by dist(A, B) = 

iTlij:^A,yeB \x - J/|. 


The following assumptions are other keys of our verification theorem, so-called cone conditions. 

Assumption 3.2. Consider (3.1). Let N C he a fast-saddle-type block such that the coordi¬ 

nate representation Nc is actually given by (2.12), and z = {x,y,e). 

Define cta“ = cta“(z), cta^ = crA“(-z), = ob“(2:), <tbj = crM^iz), ag^ = Og- 

tTgu(z) be maximal singular values of the following matrices at z, respectively: 


(z) and ag^ = 



: A^z) 

crA“ 

: A“(z) 

(Tbj 

: B“(z) 

OBJ 

: B“(z) 


: 5i(2) 


: 52"(^) 


dFi 

da 


iz) 


: u X u-matrix. 


dFi. , aBi , , dFi. . 

w*" 


30, , 


.• s X u-matrix, 


30, , 30, , 30, , 
^(*) 


{z) I : I X u-matrix, 


% 

da^ 


dg, . dg dg 

36*"* 3»<"* 


: u X {s + I + l)-matrix. 


: s X (s + I + 1) -matrix, 


I X {s + I + l)-matrix. 


Assume that the following inequalities hold: 


Aa - sup cta“ -b sup ctaj > 0, (3.4) 

Aa + \fJ-B\ - {supOA" + sup (TAJ -b SUp Ob” -b SUp CTbj + O (sup (Tg- -b SUp tTgj)} > 0, (3.5) 


where A a and pLs are real numbers satisfying (3.2) and 1(3.3\ ), respectively, and the notation “sup” 
means the supremum on N x [0, eo]. 

Assumption 3.3. Consider (3.1). Let N C he a fast-saddle-type block such that the coordi¬ 

nate representation Nf, is actually given by (2.12), and z = {x,y,e). 

Define ctaj = crAj(x), ctai = crA‘{z), obj = crBj(2 :), ctb| = crB|(^), Ogj = crgs(z) and ags = crg|(z) 
be maximal singular values of the following matrices at z, respectively: 


(dF \ 

(Taj : A^( 0 ) = ( ^ ^ s-matrix, 

fdF, dFi dF, 

"‘X A2(-’) = (aA*"> 1G*"’ W*"’ 

fdF2 \ 

OBJ : Bi(z) = ( ~^i^) j ■ ■s X s-matrix, 

(dFo dFy dFo 


: u X (u + I + 1)-matrix. 


: s X {u + I + 1)-matrix. 




■ = ( ^(^) ) ■ ^ X s-matrix. 


db 


dg, , dg 


dg, 


■ 52(-2*) = ^(-2*) ^(-2) : I X {u +I + l)-matrix. 


da 


dy dt] 


17 








Kaname Matsue 


Rigorous Numerics for Fast-Slow Systems 


Assume that the following inequalities hold: 


I/tbI - SUpOBJ + SUpCTB" > 0, (3.6) 

Aa + ImbI - {supUAJ -I- SUPOAI -I- SUpOBf S" SUp CTb| + O (sup tTg* -I- supcTgi)} > 0, (3.7) 


where A a and jiB are real numbers satisfying (3.2) and 1(3.3\ ), respectively, and the notation “sup” 
means the supremum on N x [0, eo]. 

Definition 3.4 (Cone conditions). We say a fast-saddle-type block N satisfies the unstable cone 
condition for (l.l)e if Assumption |3.2| holds. Similarly, we say a fast-saddle-type block N satisfies 
the stable cone condition for (O). if Assumption |3.3| holds. They make sense only when e G [0, eg]. 


By following the proof of Theorem 4 in 
result in this section. 


we obtain the following theorem, which is the main 


Theorem 3.5. Consider (3.1). Let N C be a fast-saddle-type block such that the coordinate 
representation Nc is actually given by (2.12) with HyN = K C 


Then 


1 . if we assume 


Aa- 


sup crA“ 


sup CTA“ + sup CTb? 


> 0, |mb|- 


SUPOBJ + 


sup cta" -k sup CTb; 


> 0, (3.8) 


there is an l-dimensional normally hyperbolic invariant manifold SoCiN in for the limit 
system such that S'o n A is the graph of a smooth function depending on y. 

2. if the stable cone condition is satisfied, there exists a Lipschitz function a = hs{b,y,e) defined 
in Bg X Bi X [0, eg] such that the graph 

VF®(S'e) := {(a, b,y,e)\a = hs{b, y, e)} 

is locally invariant with respect to (3. l]j . Moreover, hs{b,y,e) is C” for any r < oo on 
Rj X R) X [0,eo]. 

3. if the unstable cone condition is satisfied, there exists a Lipschitz function b = hy^[a,y,e) 
defined in Bu x Bi x [0, eg] such that the graph 


W^iSf) := {(a, b,y,e)\b = hu{a, y, e)} 


is locally invariant with respect to Ij3.1\ ). Moreover, hu{a,y,e) is C” for any r < oo on 
BuXBiX [0,eg]. 


As a consequence, under stable and unstable cone conditions, := VF®(S'e)nVF“(S'e) can be defined 
by the locally invariant l-dimensional submanifold inside N for any e G [0, eg]. 


Proof. Our proof is based on the slight modification of the proof of Theorem 4 in [18] for our current 
setting. The proof consists of three parts: (i) normal hyperbolicity and graph representation of 
So, (ii) existence of the graph representation of lF®(5'e) and (iii) coincidence of the graph of such 
a derived function with all of kF®(S'e). We shall prove them by tracing discussions in (TS]. Readers 
who are not familiar with the invariant manifold theorem in singular perturbation problems very 
well can compare our discussions with Chapter 2 in |18| . Here we shall derive discussions only 
for the stable manifold lT®(S'e), while discussions for the unstable manifold can be also derived 
replacing matrices in Assumption |3.3| by those in Assumption |3.2| and time reversal. Remark that 
an assumption for Nc implies that all discussions are done in terms of (2.12|. Without the loss of 
generality, we may assume N = Nc. 


(i): Normal hyperbolicity and the graph representation of Sq. 
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First of all, consider the special case e = 0. In this case, the slow variable y is just a parameter. 
Inequalities in ( |3.8[ ) yield the existence of local Lyapunov functions, which is a consequence of 
arguments in |23| with a little modifications as follows. Fix y G K. Define a function Ly : K” —>■ M 
by 

Ly{a{t), b{t)) := -|a(t) - + \b{t) - 

where (oy, by, y) is an equilibrium for (3.1) with e = 0 in at y. The existence of {ay, by, y) follows 
from Proposition |2.17| 

Then 

^(a(t), 6(t)) \t=o= -^{a - Uyf ■ {a - ay) - (a - Oy)^ • ^(a - ay) 


dt 

+ ■ (.b-by) + (b-by)^ ■ ^(b-by) 


dt 


UyJ yu UyJ-r\U Uy J Uy J 

= -{Aa + Fi{a,b,y,0))'^ ■ {a - ay) - {a - ay)'^ ■ (Aa + Fi{a,b,y,0)) 

+ {Bb + Fi{a,b,y,0))'^ ■ {b - by) + {b - by)'^ ■ {Bb + Fi{a,b,y,A)). 

Let Ofe be the fc-dimensional fc-vector and Oki,k 2 be the (/ci, fc 2 )-zero matrix. Now we obtain 
-{Aa + Fi{a,b,y,0))'^ ■ {a - ay) 

T 

b - b„ 

0i 







Ou,S 

Os,s 

Ol^s 

Ol^s 

Ou,l 

Os,l 

Oi,i 

Ou 

{zi)\a - ay 

^ + OA 5 (.zi)|a - ay 

\\b-by\ 


<—A/i|a —Uyl + cta“ (5i) |a — Uy I -I- 


2 ^A'^izi) 2 


o-Oyl +\b-by\ ) 


2 \l“ “‘yl W Uyl 

via the Mean Value Theorem, where zi = (di,bi,y,0) € N x [0, eo]- The same inequalities hold for 
(a — Uy)^ • {Aa + Fi{a, b, y, 0)). Similarly, we obtain 

{Bb + F 2 {a, b, y, 0))^ • {b - by) 


b-by 

0i 


< 6yp + aas{z2)\b - 5yp -h crB|(^2)|a - ay\\b - by\ 

< yB\b - by\'^ + aaAS2)\b - by\^ + + \b - by]"^) 


(a ay^ 

L L 

1 

/ Ou,U 


Ol^u 

Oi.ii 

b-by 

0i 

V 0 ^ 


^ + ^(^ 2 ) 


^(- 2 ) 


Os,l 

Oi,i 

Oi,i 


via the Mean Value Theorem, where Z 2 = ( 02 , 62 , y, 0) £ N x [0, eo]- The same inequalities hold for 
{b — by)’^ ■ {Bb + F 2 {a, b, y, 0)). Summarizing these inequalities, we obtain 

1 dL,, 


2 dt 


-{a{t),b{t)) |t=o 


^ I —^A + Oa“(5i) -|- 


0-A^{zi) + CTMAh) 


a — fly I + I y-B + 00= {Z2) + 


OAjC^i) -I- Ob; ( 22 ) 


\b-by 


< -Aa + sup crA“ -f 
V Afx[0.£o] 


fJ.B + sup OBJ 
. Afx[0.£o] 


SUPArx[0.£o] + SUp^><[o ep] ^ 


SUPArx[0,£o] + SUPAfx[0,£o] ' 


\b-by 


The right-hand side is strictly negative unless {a,b,y) = {ay,by,y), which follows from (3.4) and 
(3.6). Obviously, dLy/dt = 0 if and only if {a,b,y) = {ay,by,y). Therefore, Ly is a Lyapunov 
function. This observation leads to the following facts: 
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{a,b,y) = {ay,by,y) is an equilibrium for (1.3) which is unique in N with fixed y. 


The linearized matrix 


oa 


^(- 2 ) 


(^ 1 ) 


is strictly negative definite for all Z\,Z2 £ N. 


hyperbolic fixed point for the layer problem (1.3). 


5+#(- 2 ). 

This implies that (a, b, y) 


= {ay,by,y) is a 


These observations hold for arbitrary y £ K. Thanks to the Implicit Function Theorem, we 
can construct the graph of a smooth j/-parameter family of such hyperbolic fixed points, which 
is S'g. In particular, Sq has the structure of normally hyperbolic invariant manifold with graph 
representation. 


(ii): The graph representation of IT®(S'e 


We perform the modification to deal with slow directions. We choose a set K whose interior 
contains K so that its boundary is given by the condition v{y) = 0 for some C°° function v{y) 
and that v{y) satishes Vv{y) 7 ^ 0 for all y £ dK. The function v(jj) is assumed to be normalized 
so that S/i){y) = Uy is a, unit inward normal to dK. Let p{y) be a C°° function that has the 
following values 

^1 ify£iF^ 

0 if 2 / £ iF 


piy) = 


We now modify our system (3.1) by adding the term Sp{y)ny, where <5 is a positive number which 
remains to be chosen. 

We add an equation for the small parameter e. Following Jones [H], we use a multiple of e 


as the new auxiliary variable, e = ija, and append the equation rj' = 0 to (3.1), as noted in the 
beginning of this section. We then obtain the system 


a' = Aa + Fi{x,y,e), 
b' = Bb + F2{x,y,e), 
y' = W9{x, y, e) -b Sp{y)ny, 

.i = 0 - 


(3.9) 


As in (3.1) it is understood that a: is a function of a and 6 , which is already realized in Section 

[231 


and e is a function of 77 . If Theorem |3.5| is restated with e replaced by rj and proved in that 
formulation, its original version can easily be recaptured by substituting e back in. 

We define the new family of sets N by 

N := TTa,bN X k X [OA]. 


Define the set 


Fs := {(a, 5 , y, 77) | (pAt, (a, b, y, 77)) £ N for all t > 0 }, 


where is the flow generated by (3.9). In this part we shall prove that Fj is the graph of a 
function of (&, 77,77), say, a = hs{b,y,r]). Set C = and denote the crossing section of N 

at C = 6 


N^:= 


{(a,0 £^}. 


^In the lecture note by Jones m> riy denotes outward normal vector. We should remark that it is wrong. Indeed, 
we determine the immediate exit and entrance of a block later. Our claim here is that the exit is {(a, | a E d{7raN))} 

and does not contain dK. To this end, Uy should be the inward normal vector. If we prove the existence of W^(S'e), 
riy is chosen as the outward normal vector. 

^This modification guarantees that the vector field has the inflowing property with respect to {N,dK) (e.g. 
ElIlQl). Namely, at any point p E dK, the vector field at p goes inside N. This property plays an important role 
to prove the existence of center-stable manifolds via the graph representation. In the case of W"^(Se) or generally 
center-unstable manifolds, we modify the vector field so that the vector field has the overflowing property with 
respect to (N,dK). Namely, at any point p E dK, the vector field at p goes outside N or tangent to dN. In our 
case, it is sufficient to choose the unit vector Uy as the unit outward normal vector to dK and to choose u in the 
similar manner to W'^(5e). 
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We shall show that there is at least one point (a, () G N for which (a, Q) G N for all t > 0. 
To achieve this, w e us e the Wazewski Principle. If the exit is closed in dN and N is an 

isolating block for (3.9), then the mapping W : iV —>■ TV®’'** given by 


W{z) = (fe{T z = {a,b,y,r]), T (z) = sup{t > 0 | (^<,([0, t], z) n = 0}, 


is continuous. It is the general consequence of isolating blocks (see EH, Chapter 22. In [TH], such 
N is called Wazewski set). Since N is of saddle-type, the flow is repelling on the fast-exit N-f’~ for 
all 77 G [0,1] and is attracting on the fast-entrance TV'^’+ for all 77 G [0,1]. It is actually satisfied by 
choosing K D K sufficiently small, if necessary. The rest are the 77 -direction and the 77 -direction. 
Note that an appropriate choice of K enables us to construct local Lyapunov functions Ly in (i) 
for all y G K ■ 

For points on TV with y G OK, we know 


iy', ny) = e{gix, y, e), Uy) + d{ny, Uy) 


since p = 1 holds on dK, where (•,•) denotes an inner product on appropriate vector spaces. 
Setting TVf^ := sup^yjlpl, |T> 5 |}, the above can be estimated by 


{y',ny) > 5-eoMjy > 0 , 


if i5 > eoMj^. Since we can choose 6 arbitrarily, then it can be achieved. Therefore, dK n iV is a 
part of entrance. Finally, in the case 77 = 0 or 1, both of these sets are invariant since 77 ' = 0 holds 
everywhere and thus render neither the entrance nor the exit. 

The exit TV®’*** is then seen to be {(a,C) I o G d{TTaN)}. By our construction of N, for any 
the set is a ball of dimension u. Suppose that Tg H = 0. All points in then go to TV®*'** 
in finite time and hence, restricting Wazewski map to the crossing section, we obtain a continuous 
mapping 

W : N®”‘. 

If we follow this by the projection 7 ra(a, i^) = a, we see that tTq o IF maps a it-dimensional ball 
onto its boundary, while keeping the boundary fixed. This contradicts the well-known No-Retract 
Theorem. There is thus a point in T^ n N^. Since C is arbitrary, this gives at least one value for a 
as a function of {b,y,rf). We shall put it hs{b,y,ri). 


(iii): Uniqueness of IF®(S'g) as a graph of a = hg{b,y,T]). 


In the sequel we show that the graph of the above derived function is all of F^. At the same 
time, it will be shown that the function is, in fact, Lipschitz continuous with Lipschitz constant 
1. A comparison between the growth rates in different directions will be derived in the following 
proposition in our setting. Let {ai{t)Xi{t)), i = 1, 2, be two solutions of (3.91, set Aa := a\{t)—a 2 {t) 
and AC := C 2 (t) — Ci(^)- Further we define 


M{t) :=|Aa(t)p-|AC(t)|^ 


Proposition 3.6 (cf. Lemma 2 in [E]). Under assumptions of Theorem 3.5, if M(t) = 0 then 
M' ft) > 0 holds as long as the two solutions stay in TV, unless Aa = 0. 


Proof of Proposition \3.6[ Discussions in our proof are based on the proof of Lemma 2 in |18j with 
arrangements in our setting. We will make on each of the quantities (Aa, Aa) etc. The equation 
for A a is 

Aa' = A(a 2 - oi) -I- Ui(a 2 , 62 , 7 / 2 ,?? 2 o) - Ai(ai, 61 , 7 / 1 , 7710 -), 


which we rewrite as 

Aa' = AAa -I- AFi, 


(3.10) 
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where AFi = ^ 1 ( 02 , 627 2 / 2 , » 72 o) — F’i(ai, 5i, yi, T/itr). Since N is the product of h-sets associated 
with an affine transformation, thanks to the Mean Value Theorem, one can derive the following 
equality: 


dFi ~ dFi ~ dFi ~ dFi 

AFi = -^{ai,bi,yi,fii)Aa+—{ai,bi,yi,Tii)Ab+^-{ai,bi,yi,fji)Ay+—{ai,bi,yi,fji)Ar] 


for some points (oi, 61 , yi, yi) = 
We can estimate (Aa, Aa)' 
result, we obtain 


: Zi € N. Such an estimate makes sense since N is convex. 

= 2(Aa', Aa) by taking the inner product of (3.10) with Aa. As a 


( FiV riF' F)P riF 1 

(Aa', Aa) = (Aa)^ • |AAa + -^{5i)Aa -k -^{zi)Ab + -^{zi)Ay + -^{zi)Ar] j . (3.11) 

We also obtain the following equalities in the same manner: 

( r)F r)F r)F r)F 1 

(A6',A6) = (A6)^ • |RA6-k -^{z2)Aa+ -^(z2)Ab+ -^{z2)Ay+ -^(z2)At ]j (3.12) 

as well as 

{Ay', Ay) = cr(Ay)'^ • ||^(53)Aa-k ||(z3)A5-k ^{z 3 )Ay + |^(53)Ay| (3.13) 


for some Zi G N, i = 2, 3. Obviously (Ary', Ary) = 0 since ry' = 0. 

Recall that our final objective is the estimate of M'{t) = 2 {(Aa', Aa) — {{Ab', Ab) + {Ay', Ay))}. 
To this end, we estimate (3.11), (3.12) and (3.13) as the quadratic form associated with non-square 
matrices. For example, the equality (3.11) can be rewritten by the following matrix form: 

:= (Aa)'^AAa -k (Aa)^A"Aa -k (Aa)^A“AC, 


(Aa', Aa) = {Aa)^AAa + (Aoj’^A" 


/ Aa\ 
Ab 
Ay 
\ArjJ 


where A“ is the u x it-matrix given by {dFi/da){zi) and A 2 is the u x {s + I + l)-matrix given by 

Denoting || • II 2 the (matrix) 2-norm, general theory of linear algebra yields 

/ Aa\ 


(Aa)^A" 


Ab 

Ay 

\Arj/ 


<|Aa|(||A}|! 2 |Aa| + ||A“|| 2 |AC|) 


= crA 5 ‘( 5 i)|Aap -kcrAj(zi)|Aa||AC|, 

where oa" = cta'‘(5i), r = 1, 2, is the maximal singular value of A“ at zi G N stated in Assumption 


3.3 We obtain the following estimate of (3.11): 


(Aa',Aa) > (Aa - suptTA“)|Aap - supcrAj|Aa||AC|, 

N N 


(3.14) 


where Aa is a positive number satisfying (3.2). In the similar manner we also obtain the estimate 
of {Ab',Ab) and {Ay', Ay) by 

{Ah',Ab) < yslAfop -k |A 6 |(supCTB“|Aa| -k supcrajIACI), 

N N 

{Ay', Ay) < ( 7 |Ay|(supcrg 5 .|Aa| + sup | AC|), 


N 


N 
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where denotes a negative number satisfying (3.3). Functions cti" = ob" = crBj(5), 

(Tgu = crgj‘(5) and CTgu = Og^(z) acc maximal singular values of the matrix valued functions defined 
in Assumption |3.2| at 5, respectively. 

In particular, at a point satisfying M{t) = 0, we further obtain 

(Aa',Aa) > (Aa-sup(ta“ - supcrA“)|Aap, 


to show 


N 


N 


{Ab',Ab) < (/is-I-supcrB“-I-sup(TBj)|Aap, 

N N 

(Ay', Ay) < cr I Ay |(sup CTg-1 Aa I -h sup cTg-1ACI) 


N 


N 


M'{t) > 


Aa + |/.tB I - < sup CTA” -I- sup tTA“ + sup (Jb“ + SUp (78“ + O SUp (Tg- + SUp (Tg- 
\ N N N N \ N N 


|Aa|^ 


which proves our statement. Note that the coefficient of |Aap can be positive \i K Z) K \s 
sufficiently small, thanks to cone conditions. □ 

We go back to the proof in (iii). We have already shown that the set F^ contains the graph 
of a function denoted by a = hs(b,y,rf). Suppose that Fg contains more than one point with the 
same value of 6, y, y. There would then be oi and 02 so that both (oi, &, y, y) and ( 02 , b, y, y) lie in 
Fg. At t = 0, we have |Aa| > |A<^|. By Proposition |3.6[ 


|Aa(t)| > |AC(t)| 


holds for all t > 0. In the estimate (3.14) we can then replace |A()| by |Aa| to obtain 


{|Aap}'> (Aa-^)|Ac 


A positive number /3 can be chosen as sup^ cta" +sup^ ctaj- See (3.4). From which it can be easily 
concluded that Ao grows exponentially, which contradicts the hypothesis that both points stay in 
N for all t > 0. 

The same argument can be used to show that hg is also Lipschitz. If (ai,(Ci) and (02, C2) are 
both in Fg and |a2 — ai| > IC2 — Ci|) then |a2 — ai| grows exponentially, which contradicts the 
hypothesis again that both points lie in Fg. We then have shown that the set Fg is the graph 
of a Lipschitz function. In particular, the Lipschitz function a = hs{b,y,e) is well-de fined . This 
manifold is IT®(S'e) when y G K is restricted to K, in which the modified equation (3.9) agrees 
with the original one (3.11. Smoothness result follows from the same arguments as □ 


A direct consequence of graph representations of IT®(S'e) as well as IF“(S'e) is that they are 
homotopic to flat hyperplanes inside Ii-sets. In terms of the theory of covering relations (Section 
2.2), IT®(5'e) is a vertical disc and IT“(«S'e) is a horizontal disc in a given fast-saddle-type block. 


Corollary 3.7. Let N G be a fast-saddle-type bloek for (3.1) such that the coordinate rep¬ 

resentation Nc is actually of the form (2.12 ). Assume that N satisfies (3.4\ ) and ( fg.5| ). Let b = 
hu{a, y, e) be a funetion defining W^{Se) in Nc. Then there is a homotopy H : [0,1] x Ac x [0, cq] —>■ 
Nc X [0, eo] satisfying 

Ho{a,b,y,e) = H{0,a,b,y,e) = (a,/i„(a, y, e), y, e), Hi{a,b,y,e) = (a,0,y,e). 


Similarly, if Nc satisfies (S.f^ and (3. 7[ ) and a = ha{b,y,e) denotes a function defining ^^(S'c) 
in N, then there is a homotopy H : [0,1] x Ac x [0, eg] —>■ Ac x [0, eg] satisfying 


Ho{a,b,y,e) = H{0,a,b,y,e) = (/ig(6, y, e), 5, y, e), Hi{a,b,y,e) = (0,6,y,e). 


23 


















Kaname Matsue 


Rigorous Numerics for Fast-Slow Systems 


Proof. When Nc satisfies (3.41 and (3.5), define a map H on [0,1] x W x [0, cq] —>■ by 

H{X,a,b,y,e) := (a, (1 - X)hu{a,y,e),y,e). 


By tlic contr^iCtion of hut ^ contiiiuoiis. Since convex, then H nrnps [0,1] x Nc ^ [ 0; Cq] 
into Nc X [0,eo]. Thus H is our desiring homotopy. The proof of the case that Nc satisfies (3.6) 
and (3.7) is similar. □ 


Remark 3.8. We do not need to assume the normal hyperbolicity of the critical manifold Sq 
and the graph representation of Sq in Theorem |3.5[ while the Fenichel’s classical theory needs to 
assume them. Indeed, Part (i) in the proof of the theorem shows that this assumption can be 


validated by a computable estimate (3.8). 


Remark 3.9. It is well known that is not uniquely determined in general. 


Let D and D be fast-saddle type blocks given by (2.12) and (2.13), respectively. Theorem 3.5 


says that, under cone conditions, the slow manifold is contained in D. Obviously is also con¬ 
tained in D since D C D. Moreover, is uniquely determined in D. If do < dist{d{Tra,bD), TVa^bD), 
then our observations imply that the distance between and dD in fast components is greater 
than do- For example, the distance of dDc and in ui-component is greater than minlryf]^}. 
Summarizing these arguments we have the following result, which is essential to Section [4.3| 


Corollary 3.10 fFast -sadd le-tvpe blocks with spaces). Consider (2.10). Let D,D C be fast- 
sad dle-ty ve blo cks fo r (2.10) such that the coordinate representations Dc and Dc are actually given 


by (2.12) and (2.13), respectively, for a given sequence of positive numbers u j 2 =i 


Ass ume that stable and unstable cone conditions hold in Dc. Then the same statements as Theorem 
3.5 holds in Dc. Moreover, the distance between dDc and the validated slow manifold is estimated 


by 


dist(7ra(5i)c),7ra(VF^(S'e))) > min {ryX}, 

j — lc" ,U 

dist{TTb{dDc),nbiW^{Se))) > min {rjfA fore€[0 ,eo]. 

j-l,... ,s 


The main feature of our present results is that slow manifolds as well as their stable and unstable 
manifolds in Fenichel’s theorems are validated in given blocks with an explicit range e € [0, eo]. Our 
criteria can be explicitly validated with rigorous numerics, as many preceding works (e.g. [SillST])- 
We end this section noting that all results in this section holds if we replace y G M. and M-valued 
function g hy y gM’’ and M^valued function g, respectively, with / > 1. 


4 Covering-Exchange and dynamics around slow manifolds 


In Section we have discussed the validation of slow manifolds. In this section, we move to the 
next stage; the behavior near and on validated slow manifolds. There are mainly two cases of 
the exhibition of slow dynamics. One is monotone and the other is nontrivial in the sense that 
dynamics on slow manifolds admit nontrivial invariant sets such as fixed points or periodic orbits. 

First we consider the dynamics near slow manifolds. If, for sufficiently small e > 0, a point 
qc G is sufficiently close to slow manifolds, the trajectory through q^ moves near slow man¬ 

ifolds spending a long time. However, the precise description of trajectories off slow manifolds is 
not easy since they clearly have an influence of fast dynamics x = f{x,y,e). In order to describe 
such behavior as simple as possible, we introduce an extension of covering relations: the covering- 
exchange (Section |4.1[ ), so that it can be appropriately applied to singular perturbation problems. 
This extension is a topological analogue of well-known Exchange Lemma. In the successive subsec¬ 
tions, we provide a generalization of covering-exchange: the slow shadowing and m-cones (Sections 
4.2 - 4.4). These concepts enable us to validate slow manifolds with nonlinear structures as well as 


their stable and unstable manifolds in reasonable ways via rigorous numerics, keeping the essence 
of covering-exchange property. 
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Next, we consider the nontrivial dynamics on slow manifolds, namely, the presence of nontrivial 
invariant sets. We can use ideas discussed in Section 12.31 to validate nontrivial invariant sets even 
on slow manifolds. 

Finally, we discuss the unstable manifold of invariant sets on slow manifolds. The invariant 
foliation structure with respect to slow manifolds is essentially applied to validating unstable 
manifolds of invariant sets in terms of covering relations. 


4.1 Covering-Exchange with one-dimensional slow variable 

First we consider the case that the vector field near slow manifolds is monotone. In such a 
case, geometric singular perturbation theory has an answer which describes dynamics around slow 
manifolds in terms of locally invariant manifolds; Exchange Lemma (e.g. da EDI [HIES]). 

Exchange Lemma solves matching problems between fast dynamics and slow dynamics. More 
precisely, if a family of tracking invariant manifolds {Se}£>o have transversal intersection Eq Ht 
lF®(5'o) at e = 0, then a trajectory for a point G describes the match of fast 

and slow trajectories near the slow manifold in the full system for sufficiently small e > 0. 

However, Exchange Lemma requires the assumption of transversality between a tracking (locally 
invariant) manifold and the stable manifold of a slow manifold, which is not easy to validate 
in rigorous numerics. Alternatively, we consider the topological analogue of statements in (C®- 
) Exchange Lemma. 


Every trajectories in fast-saddle-type blocks leaves them through exits. Note that our desiring 
trajectories are the continuation of chains consisting of critical manifolds and heteroclinic orbits 
in the limit system (1.31. One thus expects that singularly perturbed global trajectories leave 
fast-saddle-type blocks through fast-exits. To describe our expectations precisely, we define the 
following notions. Eor our purpose, we restrict u = u{M) for h-sets to m = 1 unless otherwise 
noted. 


Definition 4.1 (Fast-exit face). Let M be an (n-l- l)-dimensional h-set with the following expres¬ 
sion via a homeomorphism cm : 


Me = cm{M) = BuX BsX [0,1] 

with u = u{M) = 1 and s{M) = s -|- 1. We say an h-set a fast-exit face of M if = u, 

s{M°‘) = s and there exist an element a G dB^ and a compact interval Kq C (0,1) such that 

cm{M°^) = {a} X RJ X Kq. 


See also Fig. 

Definition 4.2 (Covering-Exchange). Consider (1.1with fixed e > 0. Let N C be an h-set 
with u{N) = u = 1 and M be an (n -|- l)-dimensional h-set in with u{M) = 1. We say N 
satisfies the covering-exchange property for Q. with respect to M if the following statements are 
satisfied: 

(CEl) M is a fast-saddle-type block with the coordinate system 

Me = Hi X X [0,1], a G Hi, {bi, ■ ■ ■ ,bs) G Bs, y G [0,1]. 

(CE2) Stable and unstable cone conditions (Definition |3.4[ ) for ( |l.l[ )g are satisfied in M. 

(CE3) For q G {±1}, q ■ g{x,y,e) > 0 holds in M. We shall say q the slow direction number. 

(CE4) Letting be the flow of (1.1 L, there exists T > 0 such that N ^ M. 
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(CE5) M possesses a fast-exit face c with the expression 

cm(M“) = {a} X Bg X aedBi, [y~, y^] C {0 , 1) 


such that 

sup(7ryV?e(T, N) n M) < y~ if g = -|-1 
inf (7ryi^£(r, TV) n M) > ?/■*■ if g = —1 

holds. 

In this situation we shall call the pair (N,M) the covering-exchange pair for 

A brief illustration of the covering-exchange property is shown in Fig. (CE2) implies the 
existence of the slow manifold as well as a limiting normally hyperbolic critical manifold Sq at 
e = 0, as shown in Theorem |3.5| Combining with (CEl), one sees that M is repelling in a-direction 
and attracting in 6-direction. (CE3) means that slow dynamics in M is monotone. Remark that 
the monotonicity is assumed not only on but also in the whole region M. (CE4) topologically 
describes the transversality of the stable manifold kF®(S'e) of the slow manifold Se and the 6-set 
MT,N). 

For a fast-saddle-type block M of the form in Definition 
Bi X Bs X {1} and = Bi x Bs x {0} (equivalently and M^’+) the slow exit and the 

slow entrance of M, respectively. Similarly, in the case of g = —1, = Bi x Bg x {0} and 

= Bi X Bg X {1} be the slow exit and the slow entrance of M, respectively. 


4.2 


with q = -1-1, we say M^’ = 



Figure 2: Illustrations of the Covering-Exchange pair and a fast-exit face. 

The rectangular parallelepiped M is an h-set with u{M) = 1 and s{M) = 2. Arrows colored by red, 
blue and black show the fast unstable direction, fast stable direction and slow direction of (3.1) at 
corresponding boundaries, respectively. The black curve is an image ipe{T,N) of another h-set N 


for some T > 0. This figure describes N ^ M. Since the y-component of vector field in M is 
monotone and the fast component of the vector field on dM is already validated concretely, then 

(pe{T,N) keeps the covering relation N M for all t > T until N exits M from its top. A 

fast-exit face is colored by green, in which case m = 1, s = 1, T = 1. 


Define the Poincare map : M —>■ dM in M by 

P^{z) := ipg{tg{z),z), tg{z) := sup{t | tp^{[0,t],z) C M}. 

If no confusion arises, we use this notation for representing the Poincare map in a set M. 
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Notations 4.3. For y G [0,1] and any set A C let Ay := AD {y = y}, A<y := An {y < y} 

and A>y := An{y> y}. 


Lemma 4.4. Fix e > 0. Let (N,M) be a covering-exchange pair for (1.1)^ with the slow direction 
number g = -|-1 and the fast-exit face with C , where 


cm(M°’^“) = {1} X Rs X [y ,y+], [y ,?/+]c(0,l). 


Then, for each y G the restriction of the Poincare map to is a 

homeomorphism. Moreover, there exists an h-set M G M such that M M®’''*. 

Proof. Via a homeomorphism cm, we may assume that M = Bi x Bg x [0,1]. 

Since M is an isolating block for then Pf"^ is continuous on M, which is the consequence of 
standard theory of isolating blocks miiiisi]. Consider the backward Poincare map Pf^’ : M -G 
dM given by 


P^' {z) ■■= Te{T{z),z), r(z) := inf{t < 0 I (/7e([t,0],z) C M}. 

This is obviously well-defined via the property of flows. Since M is also an isolating block for the 
backward flow, then the exit of M in backward time is closed and hence P^’~ is also continuous. 
The correspondence between points on (Pj'^)“^(M®’'‘*) P {y = y} and their image is one-to-one 
and onto under Pf'^. This property holds since M G {y = y} is a. crossing section from (CE3). 
By restricting M to M>y and redefining Pf^ and Pf^’~ for M>y, their continuity yields that 
(pM)-i(^®xit)n{y = ig compact. Therefore, the restriction of Pf^ to (P/^)“^(M®^‘*)n{?/ = y} 
is a homeomorphism. 


Next consider P^-^ instead of Pf^. Since M C is a fast-saddle-type block with u{M) = 1 
satisfying stable and unstable cone conditions, then the stable manifold bF®(S'e) divides M into 
two components, where is the slow manifold in M. Let M+ be the component of M \ W^{S,:) 
containing M®’'**. Now we can construct an h-set with desiring covering relation as follows. 

Consider a section ^ = {(a, b, y) \ (a, b, y) G N}. Let G Pi be such that lV®(5'e)n{(6, y) = 

(b,y)} = {(ag,S, y)}. Since P^-’' M+ —)• P^^-^ {Mj^) C dM is homeomorphic, g{x,y,e) > 0 
holds and 

Pf^*(l,6,y) = (l,6,y) S M^'", P^>-%ay,b,y) G M*-" 


hold. Proposition 3.6 implies that zq,zi,Z2 G M with zq G kF®(5e), zi,Z 2 G C“(zo)j \T^a{zi — ZQ)\ > 
\T^a{z 2 — zo)| and T^y{zi) = TTy{z 2 ) satisfy TTy{P^-'‘ (zi)) < TTy{P^-'‘{z 2 )). Here C'^{z) denotes the 
unstable cone with the vertex z: 


C^{z) = {{a,b,y) \ |a-7ra(z)P > \b - TTb{z)\^ + \y - ny{z)\^}. 

Thus there exist points {a^,b,y) G M_^^ such that P^-’' {a^,b,y) G dM G {y = y^} with 
ay<a+ <a7^ ^ 

Let : Ps —>■ Pi be given by . Thanks to the continuity of P^-'‘, both '0=*= are 

continuous. Finally, define My by 

My := {{a,b,y) \ ip~^{b) <a< 'ij}~{b),bG PJ. 


It is indeed an h-set with 

■= {i.a,h,y) \ a = -0=^(6), 6 G PJ, Mf = {{a,b,y) \ ■0+(6) < a < tp~{b),bG dBg}. (4.1) 
Slight modifications of yield the correspondence 

'Ky{P^-^{ili~{b),b,y)) G {y,y~), 'Ky{P^-^{'il)"^{b),b,y)) G {y^,1) for bG Bg. 
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We rewrite the corresponding h-set as My again. 


Note that n = 0 by definition of Also, P^-^My) n = 0 holds 

since is on The arbitrary choice of qq G Bg leads to the statement in Proposition 


2.J 


In particular, My 


'*>« . - p“ 

=4> j\^exit Qf course. My also holds since Pt 


M>j 


is just a 


restriction of P,™ in M> 


Moreover, the union M := Ui/g[o y] satisfies M Proof of the case g = —1 is 

similar. □ 

The following proposition is the core of the covering-exchange property. 


Proposition 4.5. Let {N,M) be a covering-exchange pair for (1.1)^ with a fast-exit face M®’'**. 
Then there exists an h-set M G M such that 


vAT,-) ~ P’' 

N ’ M = 


M' 


exit 


Proof. Without the loss of generality, we may assume that M = Bi x Bg x [0,1] and = 

{1} X Bg X \y~ ,y'^ ] via a homomorphism cm, where [y~,y~^] C (0,1). Also, let q = -|-1. 


By Lemma 4.4 we can construct an h-set depending on y G [0,y ) which Pj^-covers M™^. 


Let My be such an h-set. Note that such an h-set can be constructed for all y G [0,?/ ). Now 
choose y± G [0,y“) so that 

< inf TTy[ip{T, N) H M) < sup7ry((^(P, N) n M) < y+, 

and define a set M by M := lJyG[y y+] monotonicity of g and homeomorphic correspon¬ 
dence of the flow imply that M is actually an h-set with M^ = Uj/efS- w+1 ^ where M^ are 

pM 

given by (4.1). The covering relation M M®’'** immediately holds from the proof of Lemma 


Moreover, N ^ M also holds by the construction of M and the assumption N ^ M. 


4.4 


□ 


Proposition 4.5 shows that the covering-exchange property enables us to track solution orbits 
near slow manifolds. This is just a consequence of Proposition 2.11[ 

Remark 4.6. Our implementations of covering-exchange property automatically solve the match¬ 
ing problem in singular perturbation problems. Indeed, let {N, M) be a covering-exchange pair 
with a fast-exit face M®’'**. In most cases, A^ is a family of tracking invariant manifold whose 
behavior is mainly dominated by fast dynamics. On the other hand, the behavior of points in M 
obtained in Proposition |4.5| is mainly dominated by the slow dynamics until it leaves M. If such a 
point is also on (^e(T, N), then we can prove the existence of a point G N which goes to M along 
the fast dynamics in the first stage, stays M dominated by the slow dynamics until it leaves M 
through M®’'** in the second stage and goes outside M along the fast dynamics again in the third 
stage. The trajectory through q^ is actually dominated by the full system (l.l)e. Such a behavior 
is thus considered as the match of dynamics in different time scales. 

One of other key points of the covering-exchange is that the choice of u-dimensional variables 

is changed during time evolutions around slow manifolds. Indeed, the covering relation N ^ 
M typically corresponds to validation of codimension p connecting orbits between saddle-type 
equilibria in the limit system. In such a case, the u = u{N) = p-dimensional variable are chosen 
from ?/-variables, while the u-dimensional variables of M are chosen from a-variables: fast-unstable 
directions, which is natural from the viewpoint of codimension p connecting orbits. Finally, during 

__ pM 

the covering relation M M®’'**, the u-dimensional variables retnrn to the part of y-variables, 
since M™^ is regarded as the initial data of other codimension p connecting orbits in the limit 
system. Consequently, we observe the “exchange” of t he choice of u-dimensional variables during 

pM 

the covering relation M M™^ in Proposition 4.5 In this sense, we call this phenomenon the 
covering-“exc/iarif^e”. Needless to say, one can describe such matching property without solving 
any differential equations in M. 
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4.2 Expanding and contracting rate of disks 


Assumptions of the covering-exchange in the previous subsection require that each branch of slow 
manifolds should be validated by one fast-saddle-type block. However, slow manifolds may have 
nontrivial curvature in general. It is not thus realistic to validate slow manifolds by one block with 
computer assistance without any knowledges of vector fields. To overcome such computational 
difhculties, we provide several techniques for generalizing conditions in the covering-exchange. 

As a preliminary to next subsections, we consider the expansion and contraction rate of disks 


in a fast-saddle-type block. Consider the system (3.1). Note that A in (3.1) is the u x rt-diagonal 


matrix such that all eigenvalues of A have positive real part, and that H in (3.1) is the sx s-diagonal 
matrix such that all eigenvalues of B have negative real part. Let N C a fast-saddle-type 

block satisfying stable and unstable cone conditions for fixed e > 0. Without the loss of generality, 
via homomorphism ctv, we may assume that 


N = Bu X Bs X [0,1]. 


Here we also assume that g{x,y,e) > 0 holds in N. Two cross sections x Bg x {0,1} then 
become the slow entrance and the slow exit of (pe, respectively. 

Definition 4.7. An a-Lipschitz unstable disk is the graph of a Lipschitz function tp ■ > 

Bg X [0,1] with Lip(^/>) < a. Similarly, an a-Lipschitz stable disk is the graph of a Lipschitz 
function ijj : Bg —>■ Bu x [0,1] with Lip('i/;) < a. 

Next we consider the expanding and contracting rate of stable and unstable disks. 

Lemma 4.8. Let N C be a fast-saddle type block satisfying stable and unstable cone con¬ 

ditions. Let D^{r) and D‘^{r) be a 1-Lipschitz unstable disk and a stable disk of the diameter r, 
respectively, contained in N. Let T > 0 be fixed. 


1. Assume that T], Zl“(r)) C N. Then for all zi,Z 2 S D'^{r) the following inequality 

holds: 

|7ra{zi(t) - Z2{t))\ > e^"”*|7ra(2:i - 22)1 for t G [0,T], 


where A„iin = — (supoA” -ksupoA”) > 0 given by (3.4). Here Aa denotes a positive 

number satisfying 


2. Assume that ipg{[—T,0],D‘‘{r)) C N. Then for all zi,Z 2 G D^(r) the following inequality 
holds: 

k&(zi(-t) - Z 2 i-t))\ > - Z 2 )\ for t G [0,T], 


where /Tmin = Tb + (supcrB® +supctb=) < 0 js given by (3.6). Here ys denotes a negative 
number satisfying (SJ 


The above inequalities imply that the expanding rate of D^{r) is uniformly bounded below by 
Aa and that the contracting rate of D^{r) is uniformly bounded above hy ys- 


Proof. 1. Assume first that zi,Z 2 G D"(r). The unstable cone condition with our assumption 
implies that Z 2 G C“(zi) and Z 2 (t) G C"(zi(t)) hold for all t G [0,T]. In particular, |7ra(zi(t) — 
Z 2 {t))\ > \TTb^y{zi{t) - Z 2 {t))\ holds for all t G [0,T]. Let Aa(t) := 7ra(zi(t) - Z 2 (t)), Ab(t) := 
TTb{zi(t) — Z 2 (t)) and Ay(t) := 7rj,(zi(<) — 22 ( 1 )). The first variation equation in a-direction yields 

{Aa',Aa) > A^IAap - (supoA^IAap + suptTAj|Aap), 

which is obtained by arguments in the proof of Proposition |3.6[ Here we used the fact that 
|Aa(t)| > (|A&(t)p -I- |A?/(t)|)^/^ holds in C'“(zi(t)) for all t G [0,T]. This inequality leads to 
(|Aap)' > 2Amin|Aap, and hence |Aa(t)| > exp(Aniint)|Aa(0)|. 


2. Next assume that zi,Z 2 G D‘^{r). The stable cone condition with our assumption implies 
that Z 2 € C"*(zi) and Z 2 {—t) G C®(zi(—t)) hold for all t e [0,T], where 

C^z) = {{a,b,y) \ \b - Trb{z)\^ > \a - Tra{z)\^ + \y - ny{z)\^}. 
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In particular, \Tri,{zi{—t) — Z 2 {—t))\ > |7ra,y(zi(—t) — Z 2 {—t))\ holds for all t S [0,r]. 

The backward first variation equation in 6 -direction yields 

i^|A6|^ > IfisW^bf - (supCTB5|A6|^ -I-sup(TB5|A6|^), where i = -t. 

Here we used the fact that |A 6 (—1 )| > (|Aa(— -I- |Ay(—holds in C'®(zi(—t)) for all t e 
[0,T]. This inequality leads to (|A 6 p)' > 2|/imin||A6p, and hence |A 6 (— 1 )| > exp(—/rniint)|A6(0)|. 

□ 


Note that Lemma [T 8 | holds even when u ^ 1. 


4.3 Slow shadowing, drop and jump 

Under suitable assumptions, slow manifolds as well as limiting critical manifolds are represented 
by graphs of nonlinear functions. In such cases, it is not easy to validate slow manifolds by one 
block of fast-saddle-type. It is thus natural to construct enclosure of such manifolds by the finite 
sequence of fast-saddle-type blocks, which will be easier than validation by one h-set. However, the 
union of finite fast-saddle-type blocks is not generally convex and hence we cannot apply arguments 
of the covering-exchange directly. In particular, we have to consider the effect of slow exits and 
slow entrances of the union. Nevertheless, the behavior in y-direction is much slower than the 
behavior in (a, 6 )-direction for sufficiently small e. It is thus natural to consider that slow exits 
and entrances cause no effect on validation of trajectories sufficiently close to slow manifolds for 


(1.1). Here we construct a sufficient condition to explicitly guarantee such expectation. The key 


feature consists of two parts: (i) the comparison of expansion and contraction rate with the speed 
of slow vector field, and (ii) abstract construction of covering relations around slow manifolds. Our 
proposing concept named slow shadowing exp resses both features. Slow shadowing can incorporate 


4.5 


the drop corresponding to N 




M and the jump 


with the covering relations in Proposition 

pM 

corresponding to M M®’'**. As a consequence, the concept of covering-exchange is generalized 
to the finite union of fast-saddle-type blocks. 

The center of our considerations is a sequence of fast-saddle-type blocks {A^j} satisfying all 

indicates that the stable manifold W^{Se)j 
Cone conditions also yield that these 


assumptions in Theorem 3.5 In this case, Theorem 


3.5 


is given by the graph of a Lipschitz function in each Nj. 
manifolds are patched globally in IJ Nj. 

Lemma 4.9. Let {Nj}Tj,Q be a sequence of fast-saddle-type blocks satisfying all assumptions in 


Theorem 3.5 Assume that, for all j, Nj^Nj.^.l 0 and that each section (A^ n Aj_|_i)y contains a 


unique point of the critical manifold Sq. Then, for all e G [0, eg], validated stable manifolds W^{S,:)j 
and in blocks Nj and Aj+i, respectively, coincide with each other in the intersection 

Nj n Nj^i for j = 0, • • • , ruj — 1. The similar result holds for IT“(S'e). 


Proof. The same arguments as Theorem 3.5 for (ai, 6 , y) G Nj n Aj+i n W‘^{Se)j and ( 02 , 6 , y) G 
N,nN,+,nw^{s,)j+i with a fixed ( 6 , y) yield the result. □ 


Consider two fast-saddle-type blocks Ni and N 2 . Assume that each Ni is constructed in the local 


coordinate (( 0 ^, bi),yi) whose origin corresponds to (x^, y^) such that (1.1 1 locally forms (3.1). Two 


coordinate systems {((oi, 6i), y^) | i = 1,2} are related to each other by the following commutative 
diagram: 

((ai, 6 i),yi) SEllhj. ((a 2 , 62 ),y 2 ) 


Pix7i=, 




P2 


(4.2) 


/ _ _ \ translation , _ _ n 

(x-xi,yi) -^ {x-X2,y2) 
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where Pi and P 2 : K" —)■ M" be nonsingular matrices determining the approximate diagonal form 
(3.1) around {xi,yi) and {x 2 ,y 2 ), respectively, and Im ■ K"* —> R™ denotes the identity map on 
K™. The map T 12 = o (Ti 2 )c o c^Vi denotes the composition map given by 

(( 02 ,^ 2 ), j/ 2 ) = iTi 2 )c{{ai,bi),yi) = ((T!r 42 )c(ai, &i), (T’y.i 2 )cj/i) 

:= (P2 X Ii)~^{{Pi{ai,bi),yi) + {xi - X2,yi - m)}- 


We assume 

(551) u{Ni) = u{N 2 ) = u. 

(552) iVi n A ^2 7^ 0- For given sequences of positive numbers ^ ^1 ... g, fc = 1, 2, {N^c 


is constructed by (2.13) with the sequence {Vi’^Yj^Li,... ,u,j 2 =i,-- ,s- y-component of Nk 
is given hy Try{Nk) = [yk,yk] with 

% < (^ = i> 2 ), y^G(yr>y^]: 

where tTj, is the projection onto the slow variable component. 


(SS3) All assumptions in Theorem 3.5 are satisfied in iVi and A 2 . Moreover, q ■ e{x,y,e) > 0 
holds with the slow direction number q G {±1} in iVi U iV 2 . 


(SS4) For fc = 1, 2, let Ak and Bk be diagonal matrices representing (|2.10[) in Nk, and ak and Pk 


be real numbers satisfying (3.2) for Ak and (13.31) for Bk, respectively. Also, let 


ra,k-= min Ta := min{ra,i, ra, 2 }, Ta.fe := diam(7ra(fVfe)), 

n^k ■= min n := mm{rb,i,rb,2}, fb,fc := diam(7rf,(fVfe)), 


■= sup tg{x,y,e), 
Nk X [0,€o] 


Jow 


£ 5 ( 2 ;,y,e), Cfc :=max{|e7|,|ej°“'|} 

NkX[ 0 ,eo] 


(4.3) 


Afe := ak- sup (Ta“ -I- sup < 
VAffcX[ 0 ,eo] ATfeXiO.eo] 


, y-k ■= Pk+ { sup (Tb= + sup ( 
yAffcX[ 0 ,eo] ATfcXfOjeo] 


See Assumptions 3.2 and 3.3 for the definition of oa" , oa" , obj and (Tb=. 


(SS5) Fix positive numbers da,db G (0,1) arbitrarily. Let Zl“ C Ni and C N 2 be families of 
disks given by 


(F>i)c = I IJ BYq,dbrb) \ n {Ni)c, {02)0= \ U Bu{q,dara) ^ n {N 2 )c 

/ \qew-is,)„ 


via cni and cn 2 , respectively. 
Then 

holds for all y G 'Ky(Ni) n 'Ky(N 2 ). 


{D'Yy 




(4.4) 


Since both Pi and P 2 are linear maps, the covering relation (4.4) is just a transversality of 
rectangular domains relative to an affine map. 

Definition 4.10 (Slow shadowing pair). Let x G (0,1] be a fixed number. Assume (SSI) - (SS5). 
We say the pair {A^i,7V2} satisfies the slow shadowing condition (with the ratio x OLi^d the slow 
direction number q) if the following inequality holds: 


1 , (ra.k - ra 

max — log - 

Afc V dara 


1 , f rb,k - rt 

’ l^fcl V dbXb 


h 

< X —, 

e/c 


fc = 1,2. 


(4.5) 
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We shall say the pair {-/Vi, A^ 2 } the slow shadowing pair {with x the slow direction number q) if 
{iVi, A^ 2 } satisfies the slow shadowing condition (with q). We call the number y the slow shadowing 
ratio. 

In practical computations, we set sequences of positive numbers and 

k = 1,2, as identical positive numbers, which make settings in practical computations simple. The 
assumption (SS5) admits a sufficient condition for validations in terms of cones. We revisit the 
condition in the end of Section 14.41 


The slow shadowing condition can be generalized to a sequence of fast-saddle-type blocks as 
follows. 

Definition 4.11 (Slow shadowing sequence). Consider a finite sequence of fast-saddle type blocks 
{d^i}iLo- We shall say the sequence {Ni,Xi}^o of blocks and positive numbers satishes the slow 
shadowing condition (with {xi}^o ^'^d q £ {±l}j if, for * = 0, • • • ,m — 1 and hi satisfying (SS4), 
each pair {Ni, iVi+i} is a slow shadowing pair with an identical slow direction number q in Definition 


the slow shadowing sequence (with q). 

The slow shadowing ratio y gives us a beneht in practical computations. We mention this point 
concretely in Section |6.2[ 


4.10 where Xi is the slow shadowing ratio for {A^i,7Vi+i}. We shall say such a sequence 


The core of the slow shadowing is that all disks transversal to stable and unstable manifolds of 
slow manifolds rapidly expand and contract, respectively, so that covering relation on a crossing 
section can be derived. In what follows we only consider the case q = -1-1. The corresponding 
results to = — 1 can be shown in the similar manner. 

Proposition 4.12 (Slow shadowing). Let {fVi,IV 2 } be a slow shadowing pair with the ratio y. 
Then there exist h-sets Mi C {Ni)y = Ni H {y = y} with y £ [Vi ,yi — h] fl [yi^ — yh, — h\ and 
M 2 C {N2)y+^h = d^ 2 tl{y = y + yh} such that 

Ml ‘ M2. 


Proof. For the simplicity, we assume y = 1. All arguments below are valid for general y. 

Let Se be the slow manifold validated in A^i U N 2 . For simplicity, we write Ni = {Ni)y_^_j^, 

Ni^< = (iVi)<,+^ and Pf = 

First note that Ck is an upper bound of the absolute speed in y-direction. This implies that 
any point in {Ni)y which arrives at Ni through the orbit in A^i takes at least time h/ii. Also note 
that, IF'’(5'e) and IF”(S'e) can be represented by families of 1-Lipschitz stable and unstable disks, 
respectively, which follows from cone conditions. 

Let My{Si) C {Ni)y be the (5i-neighborhood of VF®(5'e) {y = y} in {Npy, namely, 

My{Si) = {z= {a,b,y) £ Ni \ dist(z, IT^(5,)) < di}. 

Also, let M-,y{5i) be a section of My{6i) at {b,y) = {b,y), namely. 


= {z= (a, b, y) £ My{5)}. 


Obviously, all points in AT-.g is contained in the unstable cone C'^{zo{b,y)), where zo{b,y) is the 
unique point in IF®(5'e)n{(6, y) = {b, y)}. The local positive invariance of the unstable cone implies 
that zi{t) £ C'^{zo{t;b,y)), where {zo{t-,b,y)} is the solution orbit with zo{0-,b,y) = zo{b,y) and 
Zi{t) is the solution orbit with zi(0) £ M_,g((5i). This relationship holds for all t > 0 until zi{t) 
arrives at dNi. Lemma 4.8 (1) then yields 


T^aizi{t) - zo{t;b,y))\ > e^^*\Tra{zi{0) - zo{0;b,y))\, t > 0. 
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The slow shadowing condition implies that the unstable disk expands through the flow, 

and all points on the boundary arrive at before the time T = h/ei if i5i > daTa- In this case, 

P^^{My{Si)) has an intersection with N(’~ n {a = a} for all d G dBu- This observation holds for 
arbitrary b G Bg. 

The preimage Mi = {P})~^{P}{My{Si)) n tV^’<) n {Ni)y thus satisfies 

sup dist( 7 ra(z), 7 raVF'*(S'£)) < daVa- 

zeMi 


Note that is equal to Ni = {Ni)y_^_j^. 

Lemma 4.13. The set Mi is an h-set. 

Proof. We may assume that Ni = B^ x Bg x [0,1] via a homeomorphism Mi is contained in 
My{5) for some (5i G (0, daVa). For each b G Bg, consider the section M-.y{Si) = {{a,b,y) \ a G -B«}. 
Let Oj, y be such that (oy y,b,y) = zo{b, y), which is uniquely determined for each b G Bg. Note that, 
for each b G Bg, {ay y,b,y) = zo{b,y) is contained in Mi. Since Pf |(iVi) 5 : {Ni)y -^£^((^ 1 ) 5 ) C 
i9iVi < is homeomorphic, g{x,y,e) > 0 holds in Ni and 

Pg{d, b, y) = (a, b, y) G nI'^ for d G dBu, 

then for all 0 G dBu = there exists ay-{9) G Bu such that Pf{ay-{6),b,y) G NiCi n 

{a = 9}, which is uniquely determined by the property of flows. Notice that \Tra{ay -{9),b,y) — 
'^aio-y y,b,y)\ monotonically increases along the flow for all 9 G dBu. Since Pf is continuous, then 
ayy{9) depends continuously on b G Bg and 9 G S^~^. As a result, we have a continuous graph 
'ijji,y ■ Bg X 5 '““^ —>• Bu given by 0 ) = cn,^y{9). Then Mi is given by 

Ml = {(o, 6, y) I a = (1 - X)ab^y + \tpi,y{b, 9), b G Bg, 9 G A € [0,1]}, (4.6) 

which is an h-set. □ 

We go back to the proof of Proposition |4.12| The next interest is 

62 = sup dist(7rt,(z),7r6(IF“(S'£))). 

zePpMi) 


Consider the behavior of sections Mi-y, ■= {z = {d,b,y + h) G Ni \ b G Bg} n Pf{Mi) in backward 
flow. Since each Mi-y. is a stable disk with Lipschitz constant less than 1, then each point zi G Mi-a 
is contained in C‘‘{zo{d, y + h)), where ZQ(d, y + h) is the unique point in VF“(S'e) H Mi-^a- Lemma 
4^(2) yields 


kh(zi(t) - zo{t-,d,y + h))\ > e'"i*|7rb(zi(0) - zo{ 0 -,d,y + h))\ for t < 0. 

The slow shadowing condition implies that the stable disk Mi-^y. expands in 6 -direction through the 
backward flow, and its boundary arriv es at n(’'^ before the time T = h/ti if 62 > drt. Note that 
Ml C {Ni)y is an h-set given by (4.6), which implies that there exists some t G [—T, 0) such that 
the image ipe(t, Mi,a) must have intersections with A^i n {6 = 6 } for all 6 G dBg. This observation 
holds for arbitrary d with (a, b,y + h) G Ni. Therefore, 62 should be less than dbry. 


The same arguments enable us to construct an h-set M2 C (A^2)y+^ in the same way as Mi 
under assumptions of ( |4.5[ ). In this case, our assumptions and the definition of slow shadowing 
imply sup,,gjvf^ dist( 7 ra(z), < daVa. 
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Boundaries of Mi and M 2 are given by 

M^i = {(Oj 2/) I a = (1 ~ ^)^b,y + ^)) b S Bg, 0 G A G [0,1]}, 

M 2 = {{ajb^y + h) | a = (1 — A)af, y_|_^ + ^4’2,y+hib, 0), b G Bg, 0 G A S [0,1]}, 
Mf = {{a,b,y) \ a =-ipi^y{b,0), bGBg}, M^ = MinNl'^, 

M 2 = {{a,b,y + h)\a = 'il^ 2 ,y+hib), bGBg}, M 2 += M 2 0 iV/’+, 


where 1/^2 y+h ■ Bg x ^ B^ is the map associated with M 2 , which is constructed in the same 
way as 



as families of horizontal and vertical disks in and Z3|, respectively. In particular, from (SS5), 
M 2 n {Di)-^’~ = 0, where {Di)f'~ = D'l n n(’~ . This disjointness yields Pg{M^) n M 2 = 0. 
Similarly, P,^(Mi) n = 0, which yields Pl{Mi) n M+ = 0, where (£>|)^’+ = £>2 ^ . 


The rest of assumptions obviously holds if we choose go S Bg so that {a,qo,y) G Sg for some 
a G Bu- The property of degree obviously holds since Pg : {Npy {Nl)y_^_f^ is just a composite of 
uniformly contracting and expanding maps in corresponding directions. □ 


can be generalized to a slow shadowing sequence {IVi}™ 1 , which is straight¬ 
forward. 


Proposition 4.12 


The same arguments as Proposition|4.5|yield the following result: the covering-exchange : drop. 


Proposition 4.14 (Covering-Exchange : Drop). Let {Ni,N 2 } be a slow shadowing pair with the 

ratio X- Assume that there is an h-set N such that N {Ni)<y holds for some T > 0, where 

y G {Vi ^yt — b]<^ [yf ~ xb, yt ~b]- Then there are h-sets M C {Ni)<y and M 2 C {N 2 )y_^_^fi such 
that 


N M 




M2. 


Proof. As in the proof of Proposition |4.12l we may assume y = 1. 


Let Ml C {Ni)y and M 2 C (A^ 2 )y+ft be as in Proposition 


4.12 


and M := (pi^'^-®)-i(Mi). 


Obviously N M holds since N (A^i) <y. By the construction of M, the covering 




relation M 


<y + h 


M 2 also holds. 


□ 


We provide the slow shadowing when a fast-saddle type block adm its a fast-exit face, which 

^ pM 

corresponds to the covering relation M M®’'** in Proposition 


4.5 


the covering-exchange : 


jump. We restrict the unstable dimension m to 1 in the current considerations. 

Proposition 4.15 (Covering-Exchange : Jump). Let {Ni, N 2 } be a slow shadowing pair with u = 
1. Also, let C be the fast-exit face of N 2 andy G [yf ,y^ — h]r][yf —xh,y 2 —h]. Assume 

that dist(At|’'‘*, [y = y -\- xh}) > y/i. Then there are h-sets Mi C (A^i)y and M 2 C (iV 2 )y_|_^^ such 
that 

<y + xb p^2 


Ml 


M 2 


Nexit^ 


Proof. As in the proof of Proposition |4.12[ we may assume y = 1. The proof consis ts of two 

M 2 . Part (i) is a consequence of Lemma 


, p 

parts : (i) M 2 ^ and (ii) Mi ‘ 


4.4 


with 
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additional property of M 2 C {N 2 )y^h. The assumption dist(-/V|^‘*, {y = y + h}) > h and slow 
shadowing condition allow M 2 to satisfy 

sup dist(z, W%Se) r]{y = y + h}) < daVa- 

XGM2 


Therefore, the same arguments as Proposition 4.12 
Note that n M 2 = 0. See also Fig. [^(e). 




yield Mi 


M 2 , the statement (ii). 

□ 


For the convenience and correspondence to Definition 4.2 we introduce the following notion. 

Definition 4.16 (Covering-exchange sequence). Let N C be an h-set and be a 

sequence of fast-saddle-type blocks. Assume that 

(Slo-w shado-wing) {Niy,^Q^ is a slow shadowing sequence with u = 1. 


(Drop) N 




(A"f )<y holds for some T > 0, where y is given in Proposition 


4.14 


(Jump) {N^" is a slow shadowing pair with a fast-exit face of satisfying 

assumptions in Proposition |4.15] 

Then we call the triple {N, }j”o’ covering-exchange sequence. 

Obviously, the case jjq = 1 is nothing but the notion of covering-exchange pair. Remark that, 
in the current setting, covering-exchange sequences are always assumed to be defined with u = 1. 


Throughout the rest of this paper, the bold-style phrases Drop and Jump denote the corre¬ 
sponding descriptions stated in Definition |4.16| 


4.4 m-cones 

We have discussed in Section|^that it is systematically possible to construct fast-saddle-type blocks 
as well as cone conditions. However, such blocks are generally too small compared with validation 
enclosures of trajectories, if we try to validate covering-exchange sequences. Moreover, when we 
solve differential equations with a fast-exit face as initial data in this situation, solution orbits 
will hardly move in the early stage because the vector field is close to zero. This phenomenon 
causes accumulation of computation errors (e.g. wrapping effect) and extra computation costs 
(e.g. memory or time). In particular, there is little hope to validate covering-exchange sequences. 
Such difficulties can be avoided if we find a large fast-saddle type block directly. A direct approach 
would be finding crossing sections which form boundaries of a large fast-saddle type block. It is 
not realistic to find such sections via interval arithmetics because vector fields are nonlinear and 
we have to consider the effect of slow dynamics. In many cases, direct search of blocks would be 
based on trial and error, which is not systematic. Our aim in this subsection is to provide an 
appropriate method to overcome difficulties with respect to solving differential equations. 


In Section]^ we have constructed cones of the form C" = {|o — ao| > |C ~ Co|} and C"* = 
{|6 — 6o| > \v — vqW with a vertex (a, b, y) = (oq, bo,yo), where C = {b, y) and v = (a, y). See Fig. 
[^(a) for the illustration of cones and the unstable manifold of a saddle fixed point. One expects 
that, in a neighborhood of (7“, for example, (7“ can be extended to (7 like Fig. B(b)- 

On the other hand, if the unstable manifold is sufficiently regular, one also expects that such 
an extended cone can be sharper. More precisely, (7“ can be extended to the union of (7“ and 
a collection of m-cones = {|a — oo| > m|C — Co|} for some m > 1 keeping isolation, which is 
shown in Fig. ®£c). The same expectation is valid for (7^ 


Proposition 


3.6 


means that the difference of two solutions in a fast-saddle-type block N along 
dynamics is restricted by the moving cone (7“ = {{a,b,y,r]) \ M{t) > 0}. The analogue of this 
argument in the case of sharper cones is derived below. 
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(d) 



(e) 


4-12 is described by the 


Figure 3: Slow shadowing, drop and jump : schematic illustrations. 

Note that all figures here show the projection onto the {a,b)-space. 

^(«l)<5+S 

Slow shadowing ; Covering relation Mi ‘ M 2 in Proposition 

process “(a)^(b)^(d)”, where : N —>■ ON is the Poincare map in N. A set colored by blue in 

(a) denotes Mi. The Poincare map Pe maps Mi to Pe C described 

by the red set in (b) and (d). The slow shadowing condition admits the choice of M 2 in {N 2 )y.^.y 
drawn by the blue set in (d). 

3<~1><S+S 


Covering-Exchange : Drop 

is described by the process “(c)- 




M 


M 2 in Proposition 


4.14 


Covering relation N 
(b)^(d)”. In (c), the set colored by red denotes (pe{T,N) and the 


blue one denotes M. The Poincare map maps M to (M) C {Ni)y_^_f^ described 

by the red set in (b) and (d). The slow shadowing condition admits the choice of M 2 in {N 2 )yyh 
drawn by the blue set in (d). 




Covering-Exchange : Jump ; Covering relation Mi 




M2 


pN2 


in Proposition 


4.15 is described by the process “(a)^(b)^(e)”. A set colored by blue in (a) denotes Mi. The 

Poincare map maps Mi to (Mi) C described by the red set in (b) 

and (e). The slow shadowing condition admits the choice of M 2 in {N 2 )yyy drawn by the blue set 
in (e). The h-set M 2 Pf^^-covers C iV/’~ 

rectangle which admits horizontal red arrows. 


The fast exit Nf’ is drawn by the edge of white 
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Assumption 4.17. Consider (3.1). Let N C be an h-set, z = {x,y,e) and fix m> 1. 

Define = a^.{z), agl = = a-„(z) and = 

(T^(z) be maximal singular values of the following matrices at z, respectively: 


(dF\ \ 

: Ai( 2 :) = ( ——(z) j ; u x u-matrix, 


A 2 (z) = m 
Bi(z) = m 


dFr 

ay'*' 


X (s -|- 1 -|- 1)-matrix, 


dF2 

da 


iz) 


: s X u-matrix. 


(dF2, , dF2, . dF2 

^2iD=[FFriz) -ir-izj 


db 


dy 


drj 


; s X (s -|-1 -I- l)-matrix, 


'Si“ 


S2 


dg, 


gi{z) = m —(z)J ; 1 x u-matrix, 

^ If If ^^0 ^ ^ + ^ + l)-matrix. 


Assume that the following inequalities hold: 


Aa — (i 


- supo-A? +supcta; 


>) > 0 , 


Aa + Imb I - (sup CTa- + sup cta- + sup crgl -h sup + a (sup -t sup ) | > 0, 


(4.7) 

(4.8) 


where A a and hb are real numbers satisfying (3.2) and (3.3), respectively, and the notation “sup” 
means the supremum on N x [0, eo]. 


Proposition 4.18 (Unstable m-cone). Consider (3.1). Let N C be an h-set and fix m > 

1. Assume that inequalities |i(. 7j ) and in Assumption 4-17 hold. Then, letting a function 

M“’'"(t) := |Aa(t)p — m^|A^(t)f, M^''^ (t) > 0 holds for all points in N satisfying M“’'"(t) = 0 
unless Aa = 0, where ( = {b,y,r]). 


Proof. Do the same arguments as the proof of Proposition 3.6 replacing M{t) by M“’"*(t). □ 


Assumption 4.19. Consider (3.1). Let N C be an h-set, z = {x,y,e) and fix m> 1. 

Define = ^.(z), = oSl(z), = a^.{z), = a^,{z) and = a^.iz) 

be maximal singular values of the following matrices at z, respectively: 


(dF \ 

: Ai(z) = TO ( j : ux s-matrix, 

m . t : fdFi. , dFi. . 5Fi , . 


; tt X (m -I- 1 -I- 1)-matrix. 




:.iz)=^ 


db 

-1 


: s X s-matrix. 


dF,, , dF,, . dF,, , 


s X (u + 1 -I- 1)-matrix. 


91 


9i 


, . f dg , A 

gi{z) “ ^ ( ^u) j -lx s-matrix, 

^ (If If If ^^0 ^ ^ ^ + l)-matrix. 


Assume that the following inequalities hold: 
lysl - (supogl supcr£|) > 0, 

Aa + Imb I - (sup (T^s + sup a^s + sup ogi -f sup crgl -h sup -h sup cr^) >0, 


(4.9) 

(4.10) 
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where and fj-B are real numbers satisfying (3.2) and (3.3), respectively, and the notation “sup” 
means the supremum on N x [0, eo]. 

Proposition 4.20 (Stable m-cone). Consider (3.1). Let N € be an h-set and fix m > 1. 

Assume that inequalities (4-9) and I14-10 ) in Assumption 4-17 hold. Then, defining a function 
:= |A6(t)p — (t) < 0 holds for all points on N satisfying = 0 

unless Ab = 0, where v = (a, y, rf). 


Proof. Do the same arguments as the proof of Proposition 3.6 replacing M{t) by M^''^{t). □ 


Definition 4.21 (m-cone conditions). We shall call inequalities (4.7) and (4.8) in Assumption 


4.17 the unstable m-cone eondition in N. Similarly, we shall call inequalities (|4.9[) and (4.10) 


in Assumption |4.19| the stable m-cone condition in N. When these conditions are satisfied, the 
unstable m-cone and the stable m-eone with the vertex z = (ao,&Oj2/o) (in fh® (a, 6, 2 /)-coordinate) 
are given as follows, respectively: 

C(n{z) := {{a,b,v) \ |a-aoP > m^d^- + ly-yoP)}, 

:= I I&- > m^da-OoP -b |?/- J/oP)}- 

Validations of m-cones themselves are in fact independent of the construction of fast-saddle-type 
blocks discussed in Section [2.3| Moreover, the choice of m can be arbitrary as long as corresponding 
cone conditions hold. 


An implementation of m-cones in unstable direction is the following. 


1. Prepare a fast-saddle-type block N such that is given by (2.13) satisfying the unstable 
cone condition. Via a homeomorphism cjv we may assume that N is represented by 

U S 

N = l[[a- , a+] X l[[b- , b+] X [0,1]. 
i=i i=i 


2. Choose a fast-exit face iV®’''*. For example, set 

to-l u 


^exit ^ 




- „+i 


e{±}, 


i=i 


i=io-i 


i=i 


where [y , y"*"] C [0,1]. Let £ > 0 is a given number and be an h-set given by 


= l[[a-,d+]x 


n 




bj -, bj H - 

X 

y 

i 

+ 

+ 

1 ^ 

■' m ■' m 



m 

m 


(4.11) 


[«, ,a+] = 


K- 


+ £] 


.K 


if j = jo and * = +, 
if j = jo and * = -, 
otherwize. 


See Fig. |^(c). 


3. Verify the unstable m-cone condition, (4.7) and (4.8), in . 


The following lemma is a consequence of discussions in Lemma 4.18 and arguments in Theorem 


Lemma 4.22. Let N be a fast-saddle type block satisfying the stable cone condition, and be 

given by (4.11). Assume that the m-unstable eone eondition is satisfied in . Then any points 
on leaves the set 

N U := {(a, C) = (a, b, y) G | |a - ao| > m|C - Col, (ao, Co) e iV®-‘} 

under the flow through (Cififloy ;= n {a = dy}, * G {±}. The sign * is exactly the location 
of 
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TO-cones in stable direction is constructed in the similar manner. 


1. Prepare a fast-saddle-type block N such that is given by (2.13) satisfying the stable cone 


condition. Via a homeomorphism cjv we may assume that N is represented by 

U S 

N = [][a-, a+] X l[[b-, b+] X [0,1]. 
i=i i=i 

2. Choose a face of the fast-entrance As an example, set 


io-i 


iv-‘=nK”-<]x n *e{±}, 

i=i 1=1 l=lo-i 

where [y~ ,y'^] C [0,1]. Let £ > 0 is a given number and be an h-set given by 


=n 

1=1 


a, - ,a- H- 

•' m ra 


X 

1=1 


y -, 2 / + — 

m m 


(4.12) 


+ ^] J = ■^0 * = +, 

[^ 7 ’= 'I [^7 “ J = -^0 * = -. 

[6“,67] otherwize. 


3. Verify the stable m-cone condition, (4.9) and (4.10), in V^^°. 


The following lemma is a consequence of similar arguments to Lemma 4.22 


Lemma 4.23. Let N be a fast-saddle type block satisfying the stable cone condition, and 
he given by { Assume that the m-stable cone condition is satisfied in N U . Then any 

points on leaves 

N U := {{b, v) e \v = {a, y), \h - bo\ > m\v - i^ol, {bo, i^o) £ 

under the backward flow through {Cff'^)* := Cffl° n {6 = ^ '^be sign * is exactly the 

location of . 

Definition 4.24 (m-cones). We call the set of the form an unstable m-cone of N with the 
length £. Similarly, we call the set of the form a stable m-cone of N with the length 1. 

Using m-cones, we can extend an h-set after construction of a small fast-saddle type block. Such 
an h-set keeps isolation in radial direction, thanks to m-cone conditions. Moreover, we can obtain 
a priori estimates of trajectories far from equilibria, keeping their accuracy as well as possible 
by adjustments of m and £. Of course, we do not need to solve differential equations to obtain 
such estimates. This technique reduces extra computation costs mentioned at the beginning of 


this subsection. Furthermore, in layer problems (1.3), m-cones immediately yield continuation of 


stable and unstable manifolds of critical invariant manifolds. Indeed, the unions N U and 
N are h-sets satisfying m-(un)stable cone conditions. The unstable (resp. stable) manifold 

is hence represented by a horizontal (resp. vertical) disk in V U (resp. N U See [39] 

for details. 


Remark 4.25. Arguments involving Propositions 4.18 and 4.20 still holds replacing y G 
K-valued function g by y and K*-valued function g, respectively. 


and 


In what follows we show applications of m-cones to fast-slow systems. In the case of fast-slow 
systems, we have to care about movements of trajectories in slow direction. With the help of slow 
shadowing condition, we can construct a sequence of covering relations for m-cones. For technical 
reasons we restrict u and s to 1, in particular, n = 3. This is exactly the case in our verification 
examples, Section]^ Consider (3.1) again. 
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fcj 


(d) 


Figure Unstable m-cones. 


(a). A fast-saddle-type bloek N containing the slow manifold Se validated in Theorem 3.5 In this 
figure, a black ball corresponds to slow manifold Se and a blaek curve corresponds to its unstable 
manifold VF“(S'e). The accuracy of the stable (resp. unstable) manifold (resp. W^{Se) 

is measured by the size of the fast-entrance iV-f’’*' (resp. the fast-exit N^’~). In general, blocks 
are small and the flow stay near those blocks for small t, which cause the accumulation of various 
computational errors. 


(b). A candidate of extended cones. One expect that our validated cones stated in Theorem 3.5 


can be locally extended. However, one can imagine that it is quite too large for the enclosure of 
tF“(S'e) ifW^{Se) is sufficiently smooth. The extension of cones drawn here is thus quite coarse 
for smooth manifolds. 


(c). Validation of m-cone condition. It is sufficient to verify m-cone condition on a rectangular 
domain Vf) given by ( 4 . 11). The set V)( is colored by green. The orange region is the m-cone with 
a vertex in N. m-cone conditions imply that all points in an m-cone stay the m-cone until they 
leave Vf. 


(d). An unstable m-cone Of of N, which is the union of m-cones with vertices on . The 
union NUCf, is colored by yellow. All trajectories through leaves V)) through In 

general. Of H dV)( is far from slow manifolds. This fact helps us with validations with reasonable 
computation steps and accuracy. See also Section\6^ 
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Definition 4.26 (Depature time, Arrival time). Let TV be a fast-saddle-type block such is 
given by (2.131 with m = s = 1. Also, let be a set of the form (4.11) with u = s = 1 satisfying 
the unstable m-cone condition. Define the departure time Tdep in by 


2dep — 2dep(f^m) 2 


1 


d{a‘^) 


r(ra+t) 




(4.13) 


where da is a given number in (SS5), = diam(7ra(iV)), is given in (4.3), 5 > 0 is a sufficiently 

small number, 


'^min(h) I— A.4 


sup crAi^(z)+ sup a^^{z) 

<zGV^r\{a<d} V^nla^o,} 


and Aa > 0 is a real number satisfying (3.2). Note that the unstable m-cone condition implies 
that Amin (a) > 0 for a G {dara,ra + i)- 

Similarly, let be a set of the form (4.12) with m = s = 1 satisfying the stable m-cone 


condition. Dehne the arrival time T^rr in by 


Tarr — Tarr(Kn) / 

^ Ju 


(dbThY P-niin{b) 


(4.14) 


where dt is a given number in (SS5), rt = diam(7rb(TV)), ri, is given in (4.3), (5 > 0 is a sufficiently 
small number. 


Minin(&) ■= Id-B + 


sup erg*. ( 2 :)+ sup a^s(z) 
.,zev^n{b<b} z6Fin{f)<b} / 


and fiB < 0 is a real number satisfying (3.3). Note that the stable m-cone condition implies that 
Minin (b) < 0 for 6 G {dbrf,, rf, £). 

Remark 4.27. In practical computations, we use the following upper bounds of departure and 
arrival times: 


1 


r{rc,+ef 


1 


2 J (daVa)^ Amin(n) tt 




1 rirb+ef 


2 J{dbrhY Mmin(^) b^ 


min \ 

“ J-O 
T-1 

~ ^ ^ '^min [ 
i-0 

T-1 


, (j + l)^a 
T 

(j + l)^a 

T 


-1 


d{a?) 


log 


d{b'^) 1 


i=o 

T-1 


dbrb 


[3 + ^)Lb 
T 


fd I + 1 f 

2_^ Mmin ( dbrb + -=- I log I 

1=0 ^ ^ 


(dara+zd^y n 

Tdgrg + (j + l)La 

Tdgrg -\- jLg 

l-id.rb+^-^y < 62 ) 

{dbrb+^^y 

Tdbrb (j -f l)Lb 


Tdbrb + j Lb 


where Lg = fg — dgrg and Lb = rb -\- £ — dbrb- Usually, the bigger the number of partitions 
T is, the smaller right-hand sides of above inequalities are. We may set the difference between 
summations with different T’s as S in (4.13) and (4.14). 

The following results show that m-cones generate additional covering relations to validate global 
orbits for (1.1 )e. 

Proposition 4.28. Let N be a fast-saddle-type block for |'5'.l| ) in with u = s = 1 which forms 

N = [a",a+] X [6",6+] x [yo,yo] 


and actually given by (2.13). Also, let Vff, be a set of the form (4-11) with u = s = 1 which 
contains N. Assume that the unstable m-cone condition is satisfied in Vf^. Let Tdep = Tdep(V0() 
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be the departure time in V^. Define 

e+:= sup eg{x,y,e), e~ := irif eg{x,y,e) 




V^x[0,eo] 


:= {a*} X [& ,b+]x[y ,y+] with [y ,?/+]U[y -t e Tdep,y~^ + e+T^ep] C [y^ ,y^], 


(6-“ := n {«*} X 


7- ^7-4- ^ 

b - ,b+ + — 


[y -He+Tdep,y’^-He Tdep]^ , (4.15) 


where * G {±} is ehosen so that either a = a —£ora^= +£ holds. Then 

holds. 

Proof. Let zi = (ai,^i) G lF®(5'e) C N and Z 2 = ( 02 , C 2 ) G be such that Ci = C 2 and 

that 'Ky{zi) = Try{z 2 ) G [y~,y~^]. The unstable m-cone condition and Lemma 4.18 imply that 


Z 2 {t) G C'))j(zi(t)) for alH > 0 until Z 2 {t) arrives at where (Gif)* is given in Lemma 4.22 

The difference Aa = 02 ( 1 ) — ai(t) = TTa{z 2 (t)) — 'Xa{zi(t)) satisfies 

(|Aap)' > 2Ainin(a2)|Aap, 


which follows from the same argument as the proof of Lemma 4.8 This implies that the solution 
orbit Z 2 {t) with 22 ( 0 ) = Z 2 arrives at (Clf)* at time t < Tdep and that 

n{y = y+})) C (y+ + e-Tdep,y+ + e+Tdep), 
n{y = y~})) C (y~ + e~Tdep,y~ + e+Tdep). 

Thus n {y = y+}) n = 0 holds. Note that d(Clf)^^'* consists of 


(Clf) 


'U \exit,— _ 


= Clfn {a*} X 


b- —,&+ + — 
m m 


{y+ + e+Tdep}^ , 


(C“ )-+>+ = dCf, \ {a = a*Ga*} n {h*} x 


h- -,6+ 


[y + e+Tdep, y+ -H e Tdep] ^ ■ 


Moreover, (C“ )^*’+ n = 0 holds by Proposition |4.18| Then ^ (C")®’'‘‘ holds 

by Proposition 2.8 with qp G n lT“(S'e). □ 


Definition 4.29. We shall call given in (4.151 the fast-exit face of Clf. 

The following statements also hold from the same arguments under the backward flow. 


Proposition 4.30. Let Ni and N 2 be fast-saddle-type blocks for ]3.1\ ) in with u = s = 1 which 
form 

Ni = [a)-,a+] X [bf,b+] x [yf ,yt], N 2 = [af,at] x [bf,b+] x [yf ,yt] 

such that {Ni,N 2 } is a slow shadowing pair with the slow direction number q = -|-1 and h given 
below. Also, let Vff be a set of the form (f.l2) with u = s = 1 containing Ni. Assume that the 
stable m-cone condition is satisfied in Vff. Let Tarr = Tarr(V(^) be the arrival time in Vff. Define 

e+:= sup eg(x,y,€), e" := irif eg(x,y,e), 

V,^x[0,£o] Fi;x[0xo] 

iV™* := [af ,at] X {b+} X [y-,y]. with [y",y] C [yf ,yt]- 


Assume that there is an h-set Nq and Tg > 0 such that Nq ’''AA’ ’ H {y G [y , y]} and that 
y~ < inf7rj^((^e(To, A)) -|- 

^ ^arr 5 sup 7 ry((^e(To, A)) -H e+ 

^arr V- (4.16) 

Let h > 0 be such that y -|- h < y+. 

Then there exist h-sets Mi C (Ai U Cff)<y and M 2 C iN 2 )yd-h such that 


No 


Mi‘‘ M2. 
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Proof. First let is the slow manifold validated in U -/V 2 . Also, let Mi C {Ni)y and M 2 C 


4.12 


Note 


(A^ 2 )y+h be h-sets such that Mi ' M 2 , which are constructed in Proposition 

that Ml contains n {Ni)y. Let M := C Ni and M®"* := Mi n 

Clearly the set Mi contains lF®(S'e) n [Ni)<y. 

Next, let wi = {bi,vi) £ M®"* and W 2 £ W^{Se) be such that W 2 £ Cf^{wi). The stable m-cone 
condition and Lemma 4.23| imply that W 2 {—t) £ C!^{wi{—t)) for all t > 0 until W 2 [—t) arrives at 
9(7^ under the backward flow. The difference A6 = & 2 (—t) — bi{—t) = TTb{w 2 {—t)) — TTi,{wi{—t)) 
satisfies 

4(|A6n > 2 A.„i„( 62 )|A &|2 

at 


by the same argument as the proof of Lemma 4.8 where t = —t. This inequality implies that the 


backward solution orbit W 2 {—t) with ^ 2 ( 0 ) = ^2 arrives at at time t < Tarr and that 


7ryiiP^-)-\M-’-^ n{y = y})) C {y - e+Ta„,y - e-T,„), 


n {y = y-})) C (y- - e+Parr,y- - e-Tarr). 


Define an h-set Mi := M U [P?"^) where 


M+ = (Ml n {6 = 6±}) U (Mn {y = y,y-}) U (pf™)-i(M®"* n {y = y,y-}) 


Mj- = dMi \ M+. 


Note that (^e(To, Aq) n M^ = 0 follows from (4.16) and Nq 


Ve{Tn,-) 




relationship ipg{To,NQ ) n Mi = 0 immediately follows from Nq 
the covering relation Nq ^ Mi also follows from Nq ^ C {y £ [y“,y]}. 


Cm n {y £ [y ,y]}. The 
C^ n {y £ [y",y]}. Then 


□ 


Slow shadowing sequence with extended m-cones determine the generalized sequence of covering- 
exchange sequence in Definition |4.16[ 


Definition 4.31 (Covering-exchange sequence with extended cones). Consider (l.l)e in 
A C be an h-set with u(N) = 1 and {A|}^"g be a sequence such that 


Let 


{NfJy^Q is a slow shadowing sequence for (1.1 with rt = s = 1; 


• {A^” A|™} is a slow shadowing pair for (l.l)e with the exit A™'* C (A^")'f’ . 

Also, let be the unstable m'^-cone of Nf" and C^s be the stable m®-cone of A°. Assume 
that 


• all assumptions in Proposition 4.28 are satisfied with iV™'* and = C“ 

• all assumptions in Proposition 


4.30 


are satisfied with Nq = N, Ni = A® ^ and = C; 


2-1 


Then we call the collection (A, {A™*, the covering-exchange sequence with 

extended cones . 

Another benefit of m-cones is that we can validate the assumption (SS5) in slow shadowing in 
terms of m-cones. In the following proposition, we do not any restrictions on u and s. 

Proposition 4.32 (A sufficient condition of (SS5)). Consider two fast-saddle-type blocks Ni and 


N 2 constructed in the local coordinate {(ai,bi),yi), where (1.1) locally forms Il3.1\ ) in each coordinate 
with the commutative diagram 


). Assume that both Ni and N 2 satisfy the unstable mu-cone 
condition and the stable mg-cone condition for mu, mg > 1, and that (SS1)-(SS4) are satisfied. 
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Independently, consider two h-sets Z?“ and given hy 


{Dr)c 


Bu{I>,ra) X Bs 


(^,dbrb + 



:= Bu 


^ 0 ; da^a 



X Bs{0,rb) 


via cni and cn^, respectively, where Vg, rb, d a,db denote positive numbers given in (SS4)-(SS5). 
Write the nonsingular matrix P = P^^Pi in {4-^) into the block form 


P = 



Pl 2 \ 

P 22 J ’ 


where the first row (pn P 12 ) aets on the u-dimensional vector ai and the second row (p 2 i P 22 ) nets 
on the s-dimensional vector bi. Finally assume that 


D 


U 

1 



b 


s 

25 


and that 


m, ^ < 


ml ■ min - CT 21 P 2 - 

2 pl^ -I- m 2 cr 2 ip 2 


(4.17) 


where p, and pi are the maximal and the minimal eigenvalue of the matrix Pa, respectively, and 
aij denotes the maximal singular value of Pi j. Then (4-4^ holds. 

Proof. First we set 


{D'fb = dBg{0, Vg) xBbO, dbVb + — , {Df)t = R„(0, rg) x dB, 0, dbVb + — 


{b^b: ■=dBu(0:dgrg + —) xBfO,rb), {b^t ■= Bg (o,dgrg + ^] xdBfO,rb) 
\ ms J \ ms J 


corresponding to the h-set structure. 

We embed T,s,i 2 {bi) and b^ on (A^i)y and (^ 2 )^, respectively, so that the origin of both 
Tx,i 2 {bf ) and is the point g S S'g n (A^i n N 2 )y. We write the embedded sets as Tx^i 2 {bf ) and 
D 2 again, respectively. Without the loss of generality, we may assume that the representation of 
q in ( 02 , b 2 , j/ 2 )-coordinate is given by <7 = (0,0, y). Note that, in this case, Df H bf is a family of 
horizontal disks in £)“ and I?| n D 2 is a family of vertical disks in D|. 

Define two sets £>“ and i)| as follows. See also Fig. 

(.D“), := (.D“), U U (C“Jz)n{z = (ai,5i,2/i)||ai|>r4), 

Z&(D^)a 

ibbc--={bbc'J U (C'4(4n{z = (a2,62,2/2)||62|>r4). 

0 G(b|)+ 


The unstable m„-cone and the stable TOg-cone conditions imply that Df H F C Df and H F C 
Z)|, where Df and are given in (SS5) and F = {y = y}. Indeed, is given by a tube 
Bs{W'^{Sbcidbrb) in (IVi)c- The unstable manifold W'^{q) is contained in the unstable TO^-cone 
C^^(q). If IT“(S'e)c runs on the range {|a| < rg} in a-direction, the enclosure of W^{Ss)c in 
6 -direction is bounded by {| 6 | < rgjmu}. The similar arguments yield that the enclosure of the 
tube Bs{W^{Ss)c, dbVb) H {|a| < Tq} in 6 -direction is bounded by {| 6 | < dyrb + rg/mb}. Relations 
between unstable manifolds through q £ Se and (q) thus yield that Df n F C bf . The similar 

^ '^x 12 ■' 

arguments to D 2 yield the second. In particular, the covering relation '^ 02 , 62^2 

implies (SS5). Note that this argument is independent of g S S',, n IVi n IV 2 , in particular, ny{q). 


By the construction of Df and D|, 

T^aubiDf 


it is sufficient to prove the following two conditions to verify 
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Figure 5: Illustration of Proposition 


(a) : The set Df H F with u = s = 1, where F = {y = y}- The red curve represents the slice of the 
unstable manifold VF“(S'e) by F. The green tube represents the set Df H F. Definition of Df and 
the unstable mu-cone conditions indicate the inclusion Df H F C Df. Note that Df H Df consists 
of a collection of horizontal disks. 

(b) : The set Z)| n F with u = s = 1, where F = {y = y\. The blue curve represents the slice of 
the stable manifold VF®(S'e) by F. The orange tube represents the set nF. Definition o/-D| and 
the stable mg-cone conditions indicate the inclusion H F C -D|. Note that H consists of 
a collection of vertical disks. 


(c) : Covering relation Df _D| and {4.11). The inequality {4-l'7} indicates that two cones out¬ 
side Df and D| are sufficiently sharp satisfying iTTa 2 ,b 2 ° {Tx,12 X Il)){Df \ Df)nTTa2,b2{D^ \ Df) = 


(d) : Resulting covering relation (4-4'l- 
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• {'^a 2 ,b 2 O iTx ,12 X Il)){Di \ Di) n T^a 2 ,h 2{^2 \ ^2) ~ 

The former is just one of our assumptions. The latter concerns with intersections between snbsets 
of cones. Since an m-cone C^(z) is a collection of lines inside C^(z) through z, it is sufficient to 
consider the location of lines through base sets. 

Now we consider the image of the unstable m„-cone Cc = Uz6(b'‘)“ under Tc = 

{Tx, 12)0 X Ii- Since Tc is the nonsingular affine map, any lines are mapped into lines via Tc- 
Choose a point (012,612,^2) G ■ TheboundaryaCcisasubsetofU(^^^ j,^^_j^^)g(^„)-{|aii- 
ai2p = mj(|6ii - 612P -I- |yi - j/2n}- The transformation fc : {au,bii,yi) (02^, 62*, j/i) yields 

|a2i — 022^ = bii(aii — 012) + ^12(611 — 612)1^ 

= m^(|p2i(aii - 012) +^22(611 - 612)1^ -I- \yi - y2\‘^)- 


The triangular inequality yields 

|pii(aii — 012 ) + ^ 12(611 — 612 ) 1 ^ 

< |pii(aii — ^12)1^ + 2|pii(aii — 012)11^12(611 — 612)1 + bi2(6ii — 612)^ 

< 2 (|pii(oii — oi2)P -I- |pi2(6ii — 6i2)P) 

< 2 (pi^|oii — oi2p -I- cri2|6ii — 6i2p), 

|P2i(aii — 012) + ^22(611 — 612)1^ -I- \yi — j/2p 

> (|P2i(aii — 012)1 — 1^22(611 — 612)1)^ -I- \yi — 2/2^ 

= |p2i(aii — oi2)P — 2|p2i(aii — 012)11^22(611 — 612)1 + b22(6ii — 612)^ -I- \yi — 2/2P 

> |P2l(oil — Oi2)P — 2(T21p2|ail — O12II611 — 612I + p2^|6ll — 612 ^ + \yi — ?/2p 

> —cr2ip2(|aii — oi2p -I- |6ii — 6i2p) + p2^|6ii — 6i2p -|- \yi — 2/2p- 


Thus 


2 (pi^|oii — O12P -I- cr^2|6ii — 612P) 

> ml {-cr2i^(|oii - 012P -h |6ii - 6i2n -hp2^|6ii - 612P -h \yi - 2/2P} , 

further, 

(2pr^ -P m„(J2ipi) 

> - 021^- |6ii -612P -P \yi -y2p| 

> ml ■ min ^l,p2^ - cr2i^- {l^n “ + Ivi “ 2/2p} , 


equivalently, 


laii — ai2| > 


l,p2^ - Cr2lp2 - 


2pi -Pm2CT2lp2 


{|6ll - 612P -p \yi - y2p} 


Let ^ 

my min(l,p22-cr2ip^- 

771,2 ._ _1_ " 2 

2p];^ -P my2ip2 

The above inequality indicates that any line in I?“ is mapped into another one included in C^(z) 
for some z. By our construction, we know that any (half) lines in TDf \ I?“ lie on ones with 
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vertices in Tx^i 2 Di n (£)| \ Such vertices can be also chosen as points on Tx, 12 (^ 1 ) , 

which are disjoint from Remark that T'x,i 2 -D“ H = 0 follows from the covering relation 

For any point 2 : S , 12 - 0 “ n i)|, we easily know that C^(z) n C^^(z) = {z} if {mmsY > 1. 
This implies that any lines in TD'l \ Z)“ are disjoint from C^^(z) with z € T'xp 2 -D“ n D^. This 
fact holds for any 2 S T'x,i 2 -C>“ n D^. 

On the other hand, for any point w G there is a point 2 £ Tx,i 2 -D“ 0 such that 

C^^{w) C C^^{z), which follows from the property of cones and the structure of T^^uDi H 

from the covering relation £)“ D^. Consequently, we know that the set i)| \ is contained 

in the union of C^^(z) with 2 £ T'x,i 2 -D“ 0 D^- 

Combining these observations, we obtain T{D'^ \ _D“) n \ = 0. □ 

4.5 Invariant sets on slow manifolds 

Next we consider invariant sets on slow manifolds. When e > 0, the dynamics on slow manifolds 
can generally exhibit non-trivial dynamics. In other words, slow manifolds can have nontrivial 
invariant sets, such as equilibria, periodic orbits and so on, for slow dynamics. Such sets play 
a key role for describing nontrivial global solutions such as singularly perturbed homoclinic or 
heteroclinic orbits. 

Consider the slow manifold in a fast-saddle type block M with stable and unstable cone 


conditions. Theorem 3.5 implies that is represented by the graph of a function x = h'^iy) 


smoothly depending on y and e including e = 0. Substituting h'^{y) into (1.11, one sees that 


the j/-equation will decouple from the a;-equation. Since 5'^ is parameterized by y, the resulting 
decoupled equation 

y' = 

describes the dynamics on S^. After time rescaling t = tle we obtain 


y = 9ih%y),y,e), '=-^, 


(4.18) 


to show that the dynamics on slow manifolds is reduced to the regular perturbation problem. Taking 
the limit e —>■ 0 in (4.18) we obtain the dynamics on the critical manifold S'o C {/(x,y, 0) = 0}: 


y = 9{h°{y),y,0). (4.19) 

When we want to study dynamics on slow manifolds, such as the existence of fixed points and their 


stability, (4.18) as well as (4.19) are in the center of our considerations. 


In general, what we know about is just the fact that it is the graph of a function x = h'^iy), 
and its concrete description is quite difficult to obtain. As for the normally hyperbolic critical 
manifold So = {x = h^{y)}, we know that it is the subset of the nullcline {f{x,y,0) = 0}. From 
such a fact, we can study the detail of S'o such as the differential of as a y-function via implicit 
function differential equation fx{h^{y),y,0){dh°/dy) -I- fy{h°{y),y,0) — 0. On the contrary, the 
function h^{y) does not possess such simple and useful properties. For example, f{h’^{y),y,e) = 0 
does not necessarily hold. In general, h'^{y) is the solution of the following nonlinear partial 
differential equation: 

dh'^ 

<^y^{y)9{h%y),y,e) = f{h%y),y,e). (4.20) 

dy 

In order to calculate the differential dh'^/dy for studying the stability of fixed points on slow 
manifolds, for example, we have to solve the equation (4.20) rigorously. Nevertheless, from the 


viewpoint of rigorous numerics, it will be more reasonable regarding h'^ as hP with small errors (in 
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C’’-sense) than solving (4.201 directly. With this in mind, we consider (4.18) as the differential 
inclusion 

y G I e G [0,eo], G M}, 

where M is a fast-saddle type block containing S^. The right-hand side possesses h^{y) as the 
representative of the enclosure. 

As an example of dynamics on via rigorous numerics, we consider the validation of fixed 
points for (4.18) on as well as their stability for e G (0, eo]- Thanks to regular perturbation 
relationship between (4.18) and (4.19), it is natural to consider that the fixed point of (4.18) on 
S'e will be close to that of (4.19) on S'q. We then rewrite (4.18) as 

y = ff(h°(y),y,0) + {g(h'^(y),y,e) - g(h°(y),y,0)} . 

Using the Mean Value Theorem, the error term possesses the following expression: 

9 (h^{y), y, e) - 9 {h°{y),y, 0) = g{h''iy), y, e) - g(h°{y),y, e) + g{h°{y), y, e) - g{h°{y), y, 0), 

y, e) - g{h°{y),y, e) = ^(i, y, e)(y{,y) - h°{y)) for some x G M, 

9 ih°{y),y,e) - g{h°{y),y,0) = y, e)e for some e S [0,eo]. 

The difference between h‘^{y) and h^{y) can be estimated by \{h'^{y) — < Pxi- Here 

Pxi corresponds to the size of the Uth component in x-coordinate of a fast-saddle-type block M 
obtained by rigorous numerics. It is precisely computable. Such an estimate can be done because 
both h'^{y) and h^{y) belong to an identical h-set. Finally we obtain the following estimate: 

{9{h^{y),y,<^) - 9{h°{y),y,0)} G ||(i,y,e)/32,J-l,l] -b |^(/i°(y),y,e)[0,eo]. (4.21) 

Let yo be the numerical zero of g{h^{y),y, 0) = 0. The principal part 9 (hP[y), y, 0) can be then 
expanded by 

g(/i°(y), y, 0) = gih°{yo), yo, 0) -b {gxh°y + 9 v)y=y^ “ 2/o) + h.o.t. 

The Jacobian matrix Jq = g^hyG-gy at e = 0 can be explicitly calculated via the Implicit Function 
Theorem to obtain 

9xhy gy = ~gxifx) fy + ffy 

One can explicitly calculate eigenvalues of Jq, which will be close to those of the Jacobian matrix 
of Jj = dg/dy for e > 0. On the other hand, the chain rule yields 

and it contains dh'^/dy arisen in ( 4.20[ ). We additionally need to the Frechet differential of error 
terms of the form g{h'^{y), y, e) — (dh’^/dy) ■ {y — yo), when we try to validate the unique existence 
or hyperbolicity of equilibria via cone conditions |39| or Lyapunov conditions |23j . for example. 
Instead, we apply a topological tool with Jq to validation of invariant sets on slow manifolds. 


Here we construct an isolating block on slow manifolds. Using enclosures like (4.21) and eigen¬ 
values of Jo, we can compute the rigorous enclosure of vector fields on crossing sections on slow 
manifolds. 

Consider the case I = 1, which is our original setting. Assume that we have an e-parameter 
family of slow manifolds {S'e}gg(o,eo] as well as the critical manifold So, which can be validated by 
the method in Section Then one would have obtained a fast-saddle-type block M surrounding 
{S'jjggjo.eo]- schematic illustration of such a procedure is shown in Fig. Let yo G Sq he a 
numerical zero of g{x,y,0). One expects that there exists an equilibrium y^ on Se for sufficiently 
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small e > 0. We calculate the enclosure oi g{h'^{y),y,e) on {y = yo ± 6} for small d in such a 
situation. It can be done by 

gih^y), y, e) \y=y „±5 = gih°{y),y, 0) + {gih'^iy),y, e) - gih^iy), y, 0)} 

= 9ih°{yo), yo, 0) -f {gxhl + 9v)y^y^ {y - yo) + h.o.t. 

G { 9 ilT°{yo),yo, 0 ) ± JoS + h.o.t.} + + ^{h°{y),y,e)e, 

where (3x > 0 can be determined by the size of M. Remark that the enclosure of all terms in 
the right-hand side are rigorously computable and that they make sense for all e G [0, eo]- The 
principal term is zLJqS. 9{h^{yo),yo,0) is very close to 0 since j/o is assumed to be the numerical 
zero of g{h°{yo),yo,0). If 


9 {h''{y),y,€) \y=yg-s > 0 and g{h''{y),y,e) \y^y^+s < 0 (4.22) 

are validated on M x [0, cq], then the set := {ih'^iy),y, e) G M x [0, eg] \ y € [yo — <5, yo + ^]} is 
an isolating block on Se for y = y(h'^(y), y, e) with the attracting boundary. Similarly, if 


9 {h^{y),y,f^) \v=yo-s <0 and g{h%y),y,e) \y^y^+s > 0 (4.23) 


are validated on M x [0, Cq], then the set = {(h‘^(y), y, e) G M x [0, eg] | y G [yo — d, yo + d]} 
is an isolating block on for y = g{h^{y),y,e) with the repelling boundary. In both cases, the 
general Conley index theory yields that Inv(i?e) ^ 0 for all e G (0, e]. Inv(i?g) is actually the 
e-parameter family of invariant sets on slow manifolds. The direction of vector fields corresponds to 
the stability of Inv)!?,.). Moreover, Proposition 2.17 shows the existence of an equilibria in Inv(R£) 
for (4.18), which leads to an equilibrium in the full system (O- Note that isolating blocks can be 
constructed by the same implementations as Section [2.3| 


Remark 4.33. In terms of the Conley index theory, constructed sets {i3e}gg[o,eo] singular 
isolating neighborhoods. Indeed, for e = 0, Nq consists of the curve of hyperbolic equilibria and 
hence it is not an isolating neighborhood. On the other hand, for e G (0, e] it is an isolating 
neighborhood for y = y(h'^(y),y,e). 

The general theory of the Conley index can be referred in e.g. [3 [13131]. In particular, |1S| 
gives the explanation of singular perturbation version of the Conley index. More advanced topics 
for singular perturbation problems are shown in [HI HI US] for example. 


Remark 4.34. All validation of vector fields such as (4.22) and (4.23) make sense only on slow 
manifolds {S'e}£g[o,eo], since the decoupled vector field (4.18) makes sense only on slow manifolds. In 
other words, estimates of vector fields discussed here does not always give reasonable information 
of vector fields off slow manifolds. Calculations of vector fields discussed here and Section 14.11 
should be thus considered independently. 

A simpler way to study the dynamics on slow manifolds will be to focus on the dynamics (4.19), 
namely, neglect the error term g(lf{y),y, e) — y(h°(y), y, 0). We then study the dynamics (4.18) on 
slow manifolds via the regular perturbation theory. The essence of this idea is the same as ours, 
but such results make sense only for sufficiently small e > 0. One of our main purpose here is to 
study the slow dynamics for all e G (0, eo] and hence such a simpler idea is not sufficient to our 
study. 


4.6 Unstable manifolds of invariant sets on slow manifolds 

We discuss the unstable manifold of invariant sets in the full system. Combining covering-exchange 
property with graph representations of locally invariant manifolds, the standard consequence of 
covering relations. Proposition 2.14[ yields the existence of connecting orbits between slow man¬ 
ifolds for (l.l)e. However, it does not always mean the existence of connecting orbits between 
invariant sets for (l.l)e. This difference comes from the fact that the meaning of slow variable y 
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Figure 6: Invariant sets on slow manifolds. 

(a). Critical manifold Sq (e = 0, colored by ash), and perturbed slow manifold (e G (0, eo]) 
colored by black) via Theorem 3.5 The yellow rectangle with red arrows describes a fast-saddle- 
type block with the fast-exit. The critical manifold So is a subset of nullcline f(x,y,0) = 0. In 
general, Se exhibits a nontrivial structure such as equilibria (drawn by a black ball). 


(b). Validation of vector fields on Se (e G [0, eo])- One cannot detect the position of the slow 
manifold in advance. In order to detect vector field on such invariant manifolds, one put a 
section S (colored green) having an intersection with and validate the vector field on S via 
rigorous numerics. This contains information of the vector field on S D Sg. 


(c). An illustration of isolating block containing an equilibrium on S,,. Analysis in (b) contains 
the information of rigorous vector field on S^, which gives an isolating neighborhood colored by 
red. General topological arguments such as the mapping degree or the Conley index guarantee the 
existence of nontrivial invariant sets. In the case of illustration described here, there exists a fixed 
point on Se, namely, an equilibrium of the full system (1.1). 


changes between e = 0 and e > 0. When we apply covering relations to connecting orbits between 


equilibria for (1.1 )o, namely (1.3), a natural approach is to prove the existence of orbits connecting 


points on the stable and unstable manifolds of equilibria at a certain y. This immediately yields 
validation of a connecting orbit at y = y. On the other hand, if e > 0, a parameter y becomes 
a time dependent variable and slow manifolds do not always consist of a family of equilibria or 
general invariant sets. In fact, if unstable manifolds of normally hyperbolic invariant manifolds are 
perturbed for e > 0, then the negatively invariant manifolds (e.g. unstable manifolds) of invariant 
sets are drastically collapsed in general. In particular, dimensions of h-sets corresponding to un¬ 
stable manifolds of invariant sets change. The standard arguments of covering relation is not thus 


sufficient to validate connecting orbits in (l.l)e. 


To overcome this difficulty, we apply further structure of slow manifolds, called Fenichel fibering, 
to validate (un)stable manifolds of invariant sets on slow manifolds. 

A direct consequence of discussions in the previous subsection is stated as follows. This result 
can be generalized in arbitrary dimensions, but we omit the detail because we do not make such 
generalized arguments in this paper. 


Lemma 4.35. Let N be a fast-saddle type block for (l.l)e satisfying .stable and unstable cone 
conditions, and C intM be a subset of 1-dimensional slow manifold in N. Assume that is 
an isolating block on the slow manifold whose boundaries are repelling (cf. ). Then a trivial 

collection {56} is admissible with respect to (pfit, ■) for all sufficiently small t > 0 (cf. Definition 

IjI ). 

Similarly, assume that is an isolating block on the slow manifold whose boundaries are at¬ 
tracting (cf. !j4-22)). Then a trivial collection {iSe} is back-admissible with respect to (pfit,-) for 


all sufficiently small t > 0. 

Here we regard the flow ip„ in the full system as restricted to fi in both cases. This restriction 
makes sense since slow manifolds are locally invariant. 
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Since 1^ is 1-dimensional, if 1^ is repelling, the covering relation 




/e is nothing but 


le C ipe(t,Ie)- Similarly, if is attracting, the back-covering relation jg nothing but 

f 

It C (pt{—t,It). For simplicity, we write fc-covering relations N N for an h-set N and a 
continuous map / by |a^ -/vj . Let |-/V fv| be defined by the similar manner. 

V > /c ^ J k 

Lemma [4. 35| gives the other proof of existence of an invariant set in 1^. Note that this invariant 
set is actually invariant for the full system (l.lb. 

A key essence of Fenichel fibering is cone conditions, which is the same as invariant manifold 
theorem. On the other hand, our statements in Theorem |3.5| involve cone condition which are 
computable via rigorous numerics. Theorem |3.5| thus yields the following fiber validations. 


Corollary 4.36. Consider ]3.1\ ). Let N C be a fast-saddle-type block such that the coordinate 
representation Nc is actually given by (2.12) with TTyN = K C Assume that stable and unstable 
cone conditions are satisfied in N and let St be the validated slow manifold. Then for each v € intiFe, 
the following statements hold. 

1. There exists a Lipschitz function (a, y) = hs{b,e) defined in a neighborhood of Trb(v) in Bg 
and e G [0, eo] such that the graph 

:= {(a, b, y, e) \ (a, y) = hs{b, e)} 
is locally invariant with respect to 13.T(. 


2. There exists a Lipschitz function {b,y) = hu{a,e) defined in a neighborhood of Tra{v) in 
and e G [0, eo] such that the graph 

VF“(u) := {(a, b, y, e) | (b, y) = hu{a, e)} 

is locally invariant with respect to 


These functions are well-defined in a neighborhood of each v, which means that, for example, W'^(v) 
may have an intersection with dN \ NI’~ . In other words, hu is not always a horizontal disk in 
N in the sense of Definition \2.lS[ 

In the case of fast-slow systems, the property that becomes a horizontal disk in N is non¬ 
trivial, which is due to dynamics in slow direction. 

From now on we relate the fiber bundle structure of slow manifolds to covering relations. Fiber 
bundle structure of IT“(/e) for A C leads to construct a deformation retract r\ of unstable 
manifolds so that the fiber bundle IF“(Je) r>,-covers a fast-exit face. Namely, the following lemma 
holds. 


Lemma 4.37. Consider (3.1). Let N C he a fast-saddle-type block with u = 1 such that 

the coordinate representation Nc is actually given by Nc = Bi x Bg x [0,1]. Also, let TV®’'** be a 
fast-exit face of N with 

^ 1^1 ^ (iV®’'^*)- = {y±} C [0,1] wtth 0 < y_ < y+ < 1 

and Ig C intA^ be a connected subset of 1-dimensional slow manifold in N. 

Assume that, for any z G fi, the Lipschitz function (b, y) = 6^ ^(a) is a horizontal disk in N. 
We also assume that, at z± G dfi with z+ > z_, TTybuffil) > y+ and TVybuffil) < y_ hold. 

rio&* g 

Then /g where rxo ^ is the deformation retract given by 

{rx o b^cfia) := (A + (1 - A)a, b^JX + (1 - A)a)), A G [0,1], z G (4.24) 

If TTaN^^'*^ = {—1}, then the same statement holds by replacing r\ in f4-^4 ) by 

{rx o blc){a) := (-A + (1 - A)a, b^^fi-X + (1 - A)a)), A G [0,1], z G fi. (4.25) 
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Proof. Our assumptions imply that r\ o bu%{Bi) n = 0. Note that graphs of bu% lie on the 
unstable manifold hF“(/e), which implies rx o 6^ H (iV®’'‘*)+ = 0 for all 2 G h- The rest of 
assumptions in Proposition |2.8| or |2.9| is easy to check. □ 

If we additionally assume that /g is a singular isolating neighborhood, then we can describe 
unstable manifolds of invariant sets. 


Proposition 4.38. In addition to assumptions in Lemma \4.3'l[ we assume that is an isolating 
block on the slow manifold. Then there exists a point z G Inv(/£) such that W'^{z) O TV®’'** ^ 0. 


Proof. Since 1^ is isolated, Lemma 


4.35 


k > 1 and all sufficiently small t > L) it 


yields the covering relation < I, 




for arbitrary 


/<; is repelling. Similarly, if is attracting, the covering 

relation j. holds for arbitrary A: > 1- 

We only deal with the attracting case. By Proposition |2.11| there exists a point Z}. G 1^ such 
that 

forallA:' = 0,l,--- ,A:, IF“(zfe) n TV®’'** ^ 0. (4.26) 

Since 1^ is compact, then the sequence {zk\k>i contains a convergence subsequence. Let Zoo be 
the limit. We can choose a sufficiently small > 0 such that te, 0), z) fl /e = 0 holds for 
all z G die, since 1^ is an isolating block. Note that die means the boundary with respect to 
the slow manifold in TV. Obviously we can choose a limit point Zoo from (4.26) with t = te. 
The property of te,Zoo and (4.26) thus imply g}e{[—kte,Q\,Zaa) C le for all A: € N, which yields 
oo, 0], Zoo) C le. The inclusion oo), Zoo) C le is obvious since le is attracting. As a 

consequence, Zoo G Inv(/e). By continuity of the graph 6^ ^ with respect to z, then W“(zoo) is 
well-defined and W“(zoo) O TV®’'** ^0. □ 


The validation of A being an isolating block on slow manifolds can be done by arguments in 
Section [T5| On the other hand, if we apply Lemma [4.37[ we have to know the location of fibers 
IF“(z) for all z G le. The unstable (TO-)cone condition indicates that, for z G le, the fiber W'^(z) 
is included in the unstable cone C^(z). We can apply this cone to validate assumptions in Lemma 
|4.37[ Namely, the following lemma holds, which can be validated by rigorous numerics. 


Lemma 4.39. Consider (3.1). Let TV C K"’*'* be a fast-saddle-type block with u = 1 such that 
the coordinate representation TV^ is actually given by TV^ = Bi x Bg x [0,1]. Also, let TV®’'** be a 
fast-exit face of TV with 


^aTV®’'** = {!}, (TV®"**)- = {j/±} with 0<y_<y+<l 


and le Cl intTV be a connected subset of 1-dimensional slow manifold in TV. Assume that, for some 
m > 1, the unstable m-cone condition is satisfied in TV and that, at z± G die, the unstable m-cones 
Cfn(.z±) satisfy the following properties: 

(C“(z±))n9TVcTV^’-, (4.27) 

inf(7ry(C'“(z+))n{a = l}) > 2 /+, sup(7ry(C'“ (z_)) 0 {a = 1}) < ?/_. (4.28) 

Then all assumptions in Lemma jf.STj hold. 

Proof. The unstable cone condition indicates, as mentioned, that the graph of the Lipschitz func¬ 
tion b(^ e is contained in C(((z) for all z G h. Since cones have the identical shape for z G le and A 
is 1-dimensional, then the condition {{a,b(^ e) I a S Bi} C N for all z G le obviously follows from 


4.27). Inequalities Trybufe{l) > y+ and TTybu^eil) < y- are direct consequences of (4.28). 


□ 


5 Topological existence theorems of singularly perturbed 
trajectories 

Now we are ready to state our verification theorems. In this section we derive existence theorems 
for periodic and heteroclinic orbits which can be applied to 
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5.1 Existence of periodic orbits 

First we state the existence theorems for periodic orbits near singular orbits. 


Theorem 5.1 ( Exis tence of periodic orbits.). Consider (1.1). Assume that there exist p € N, 
eo > 0 such that (1.1) admits the following e (g [0, Cq])- parameter family of sets in ]R"+^ (see also 
Fig. \^(a)): 

S^: j = Q, ■ ■ ■ , p: a fast-saddle-type block satisfying the covering-exchange property with respect 
to defined below and for > 0. 

Ti'. j = 0, • • • , p: a fast-exit face of S(. Identify with Tf. 


Then the following statements hold. 

1. When e G (0, Eq], Q admits a periodic orbit {rg(t) | t G M} with a sequence of positive 
numbers 

0 < < tg < < • • • < = Te < oo 

such that 


r,(r,) = r,(o) g 

r,{^,F+^])cSt\ r,(ti+i) G (j=0,--- ,p-l) and 

rd[t^,7;])c5°. 


2. When e = 0, (1.1) admits the collection of heteroclinic orbits {F-^ = (t)}tgR}^^Q with 


^ ^ ‘^0 /fe : % > 0 ) = 0 , 

Xy^ (t) = qj G with f{qj,yj,0) =0 (j = 0, • • • , p - 1), 
lim Xy (t) =ppG S(( with f{pp, pp, 0) = 0, 

c—>• —oo 

lim Xy (t) = qp e Sq with f{qp, ?/„, 0) = 0 and 

t —>-+oo ^ 

{j = o,...,p) 


bridging 1-dimensional normally hyperbolic invariant manifolds C 5g and C 
for j = 0, • • • ,p — 1. In case that j = p such a heteroclinic orbit connects Sq C Sq and 
5° C 5°. 


Proof. First consider the case e G (0, cq]. By Proposition 4.5 and our assumption, there exist h-sets 
S( C Si such that we get the sequence of covering relations 


j-p ^0 




vATa,-) 


S} 




SC 


pp 


J-f, 


where Pf : S( —?> dS( is the Poincare map of S^. Then our statement is just the consequence of 
Proposition |2.11[ The case e = 0 is just the consequence of Proposition |2.14| All arguments with 
respect to horizontal and vertical disks are valid from Corollary |3.7[ □ 


Remark 5.2. Standard Exchange Lemma says that the period when trajectories stays near slow 
manifolds in 0(e)-distance is 0(l/e) (actually 0(exp(—c/e)) for some c > 0, discussed in [TU]'). 
The period is thus expected to be 0(l/e) for e G (0,eo]. 

Theorem |5.1 1 can be generalized as stated below, replacing covering-exchange pairs by covering- 
exchange sequences. 


Corollary 5.3 (Validation of periodic orbits with slow shadowing). Consider (1.1). Assume that 
there exist p G N, eg > 0 such that (1.1) admits the following e (g [Q, Cq]) -parameter family of sets 
in (see also Fig. 
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(a) 


Figure 1: Illustration o/{5^ , in Theorems 



5.1 


and 


5.6 p = A. 


Rectangular parallelepipeds (colored red and yellow) are fast-saddle-type blocks with fast-exit faces 
(colored by green) satisfying the covering-exchange property. Vertical black arrows show the slow 
vector field. Bold black curves outside blocks describe heteroclinic orbits {F^j^^g. Black lines 
inside blocks correspond to limit critical manifolds. Their union generates a singular limit orbit 


Fg. Theorem 5.1 claims that there exists a family of periodic orbits {r£}jg(g gg] near Fg. 


/ NldV-'^r-- ,p 

J 2=0,••• 

below. 


a sequence of h-sets which forms a covering-exchange sequence with F) ^ defined 


j = 0, • • • , p: a fast-exit face of iVe. Identify ^ with Ff. 


Then all statements in Theorem 5.1 hold with S) = u: =0 


Proof. We only consider the case e G (0, eo]- By Propositions 4.12 4.14 4.15 and our assumptions, 
there is a sequence of h-sets with C which admits a sequence of covering 

relations 


Ff 




p0,0 p0,l 


pO,mg-l 




pO,mo 


F\ 


VeiTn, 


M}’° 


Mf'' 


FP 


where , is the Poinc are map in with approprl 

ate choices oi y = yj^i and h = hj^i following Propositions 
is just the consequence of Proposition|2.11[ 


4.12 


4.14 


and 


4.15 


Then our statement 

□ 


Note that points on the slow manifold S) are contained in all h-sets {MF) in the proof of 
Corollary 5.3 This observation implies that true orbits in the full system 0) shadow trajectories 
on S{ via covering relations (cf. [H]). This fact gives us a suggestion to describe true trajectories 
for fast-slow systems with multi-dimensional slow variables from the viewpoint of shadowing. 

Corollary [5^ can be further generalized as stated below, at least in the case u = s = 1, replacing 
covering-exchange sequences by covering-exchange sequences with extended cones. 


Corollary 5.4 (Validation of periodic orbits with slow shadowing and m-cones). Consider (1.1). 


Assume that there exist p € N, cq > 0 such that (1.1) admits the following e (s [0, eo])-para77ieter 
family of sets in (see also Fig. [^; 

{iVh*}l“g’ ■: a sequence of fast-saddle type blocks which admits a sequence of faxt-exit faces 

|^yexitp=0. -.p jVi.exit ^ 

The collection {F)~^, ... C“ «, forms a covering-exchange sequence 

1 1 J j j 

with extended cones . Associating sets and F) are defined below. 
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Figure 8: Illustration o/ 


5.3 p = 2. 


in Theorem 

Rectangular parallelepipeds (colored red and yellow) are fast-saddle-type blocks with fast-exit faces 
(colored by green) satisfying the covering-exchange property and slow shadowing. Vertical black 
arrows show the slow vector field. Bold black curves outside blocks describe heteroclinic orbits 
{F-^l^^o curves or surfaces inside blocks correspond to limit critical manifolds. Their union 
generates a singular limit orbit Fq. Theorem \5.1\ claims that there exists a family of periodic orbits 
{rjee(o,eo] nearTo. 




= the unstable m^-cones of Identify ^ withCf^u. 


Cm” 0 = Or • • jP)= the stable m^-cone of Nf°. 

F{ {j = 0, • • • ,p): a fast-exit face ofClLu. Identify withFf. 

ILj / ILj 


Then all statements in Theorem 


5.1 


hold with S) = U Ufco U 


5.3 




.exit 


F) and respectively. 


Proof. Replace F) and in the proof of Corollary 
Here —)■ dC!)^^ is the Poincare map and C UC(^s) is the h-set corresponding 

j ,_L ' 

to Ml in Proposition 


4.30 


□ 


Remark 5.5. The key essence of our validation near slow manifolds consists of the following three 
pieces: 


• Construction of fast-saddle-type blocks; 


• Stable and unstable cone conditions; 

• Slow shadowing condition (if necessary). 

As long as validations of the above procedures pass, we can extend slow manifolds in an arbitrary 
range keeping the existence of points near slow manifolds which exits their neighborhoods after 
time T = 0(e), as stated in the Exchange Lemma. 


5.2 Existence of connecting orbits 

Similar settings and arguments to the previous subsection yield the existence of heteroclinic orbits 
near singular orbits. 


Theorem 5.6 (Existence of heteroclinic orbits). Consider {Op. Assume that there exist p G N 
and eo > 0 such that (1.1) admits the following e (g [Q^cq])- parameter family of sets in (see 

also Fig. 
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S^: j = 0; a fast-saddle-type block satisfying stable and unstable cone conditions. Moreover, it 
contains the nonempty maximal invariant sets Se,u- 

j = p: a fast-saddle-type block satisfying the covering-exchange property with respect to 
(defined below) except condition 3. Moreover, it contains the nonempty maximal invariant 
set S^^g. The set S^^g is contained in an attracting isolating block on the slow manifold in Sf. 

j = 1, • • • , p—1: a fast-saddle-type block satisfying the covering-exchange property with respect 
to defined below and for > 0. 

T): j = 0, ■■■, p — 1: a fast-exit face of S). 

If j = 0, the invariant set admits an isolating block B{Se^u) on the slow manifold in 
such that all assumptions in Lemma \4-3S\ holds. 

Then the following statements hold. 

1. When e G (0, Cq], 0 admits a heteroclinic orbit {re(t) | t G K} with a sequence of positive 
numbers 

0 < < • • • < t(,~^ < t^j~^ < oo 

such that 


dist(rE(t), —>■ 0 as t —> —00, 

r,((-oo, 0]) c c 5°, r,(o) e iF°, 

r,([tf \ti]) c Si, TM) e (j = 1 ,• • •,p- 1 ), 

oo)) C VF"(S',,,) C Sf and 
dist(re(t), S'e.s) —>■ 0 as t —>■-too. 


2. When e = 0, (1.1) admits the collection of heteroclinic orbits {F-^ = (t)}tGR}^=o 


lim Xy. (t) = pj G S( with f{pj,yj,0) = 0, 

t — y —C50 

= % G with f{qj,yj,0) = 0, 

,p-l) 

bridging 1-dimensional normally hyperbolic invariant manifolds C 5g and C 
/or j = 0, •• • ,p-l. 

Proof. Consider the case e G (0, eg]. We deal with the case that B{Sg^u) is attracting. By Propo¬ 
sitions |4.5[ |4.38| and our assumption we get the sequence of covering relations 




TiOb^ 




S 


1 


{Tp I ,■) 




for all sufficiently small t > 0 and arbitrary k G N. Here r> is the deformation retract given by 
(4.24) or (4.25), 6^,, is the horizontal disk at z G B(Sg^u) whose graph is the unstable manifold 
W’^(z) and P) : S) -G dS) is the Poincare map of S). The stable manifold W^{Se^g) is given by 
a vertical disk 6® in Sf with (6, y)-coordinates as the s(5f )-dimensional direction of the h-set Sf. 
Our statement is just the consequence of Propositions |4.38| and |2.14| See Dehnition |2.13| about 
horizontal and vertical disks. All arguments with respect to horizontal and vertical disks are valid 
from Corollaries 3.7 and 4.36 In the case that B{Sg^u) is repelling, the same arguments valid 

replacing the covering relation |H(S'g_„) B{Se^u)'^ by ^B{Se^u) H(S'e^„)| . 




The case e = 0 is the same as that of Theorem |5.1[ Remark that we do not need the sequence 
of covering relations 

to prove our statements in this case. □ 
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Theorem |5.6| can be generalized as stated below, replacing covering-exchange pairs by covering- 
exchange sequences. 


Corollary 5.7 (Validation of heteroclinic orbits with slow shadowing). Consider (1.1). Assume 
that there exist p € N and eg > 0 such that (1.1) admits the following e {g [0, Cq])- parameter family 
of sets in ; 


{Nf’‘}i^Q'(( with mg = 0. j = 1, • • • ,p — 1 .• a sequence of h-sets which forms a covering- 
exchange sequence with defined below. 

j = p : a sequence of h-sets which forms a covering-exchange sequence with defined 

below except the last assumption in Definition \4-.16 

Blocks iVP’° and contain nonempty maximal invariant sets and respectively. 

The invariant set is contained in an attracting isolating block on the slow manifold in 

j^p.mp 


J,^3 


j = 0, - ■ ■ , p — 1: a fast-exit face of 

V j = 0) ^^6 invariant set S^.u admits an isolating block B(Se,u) on the slow manifold in 
such that all assumptions in Lemma 4-39 holds. 


Then all statements in Theorem 


5.6 hold with S( = Ui=c 


Proof. We only consider the case e G (0,eg] and B{Se^u) is attracting. By Propositions 

,p 


4.14 


4.15 


and our assumptions, there is a sequence of h-sets 


wTuch admits a sequence of covering relations 


1 = 0 .- 
0 ,- 


4.38 


4.12 


with C 




riob* 








pl,0 


M. 


1,1 



vATi,-) 





^l.mi 





for arbitrary A: G N and sufficiently small f > 0. Here r\ is the deformation retract given by 
(4.24) or (|4.25), ^ is the horizontal disk at z G B{Se,u) whose graph is the unstable manifold 
^z). Also, Pr ■ is the Poincare m ap in with 

appropriate choices of y = ijj^i and h = hj^i following Propositions 


4.12 


4.14 


and 


4.15 


The stable 


manifold lV®(S'e^s) is given by a vertical disk in iVf’”*'’ with (6, j/)-coordinates as the s{Nf’™''’)- 
dimensional direction of the h-set . Our statement is just the consequence of Propositions 


4.38|and|2.141 All arguments with respect to horizontal and vertical disks are valid from Corollaries 
TTland OBl □ 

Corollary |5.7| can be further generalized as stated below, at least in the case u = s = 1, replacing 
covering-exchange sequences by covering-exchange sequences with extended cones. 


Corollary 5.8 (Validation of heteroclinic orbits with slow shadowing and m-cones). Consider 
(1.1). Assume that there exist p G N and eg > 0 such that (1.1) admits the following e (g [0,eg])- 
parameter family of sets in .• 


with mg =0. j = 1, • • • ,p — 1 .' a sequence of fast-saddle type blocks which admits 

a sequence of faxt-exit faces ,p j\ri,exit Moreover, the collection 

, {Nfp}i^Q ... m , A"/’®’'**, 0“ u, forms a covering-exchange sequence with extended 
cones . Associating sets Cf^-u., and Tl are defined below. 

j j j j 

j = p : a sequence of h-sets which forms a covering-exchange sequence with cones with Tp~^ 
and Ci,s defined below, except the last assumption in Definition 

P 

Blocks and contain nonempty maximal invariant sets and respectively. 

The invariant set is contained in an attracting isolating block on the slow manifold in 

j^p.irip 


4.16 
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C“ u (j = 1, • • • j P ~ 1)= the unstable m^-cones of . 

‘'j 

GLs (j = 1 , • • • ,p): the stable m^,-cone of 
H O' = 0 > • • • , P - 1 ): a fast-exit face (C^u of C" „. 

3 3 

V j = Oj the invariant set Se,u admits a n iso lating block B(Se,u) on the slow manifold in 
such that all assumptions in Lemma f.SQ holds with the fast-exit face of N^’^. 


Then all statements in Theorem 


hold with 


Si = {Cf^sUijTJoNf^ 

n. 


tfj = 0 , 

tfj = p, 

otherwise. 


Proof. Replace in the proof of Theorem |5.7| by and Mf^, respectively. 

Here —>■ 9(7“,, is the Poincare map and Mfo {y UC^s) is the h-set corresponding 

'‘'j .7 '‘'j 

to Ml in Proposition 4.30 □ 


Remark 5.9. Changing the choice of covering-exchange sequences and fast-saddle-type blocks 
containing nontrivial invariant sets, we can obtain various type of singularly perturbed global 
orbits near singular orbits. For example, in Theorem 5.6 further assuming 5° = Sf and = 
Se,s = {Pe}i an equilibrium, then there exists an e-family of homo clinic orbits {i7e}eg(o,eo] ofpe. 


We can replace Condition 4 in the covering-exchange property, 
of covering relations 




51+^, by a sequence 






M( 




Mi 


i vAn,-) 


Ml Sl+^ 


5.1 


for h-sets and positive numbers to prove the same statements as Theorem ^ 

and |5.6[ These are applications to Proposition |2.11| or |2.14| and such extensions are useful for 
validating trajectories with complex behavior. 


6 Sample validation results for FitzHugh-Nagumo system 

In this section we provide several examples of singularly perturbed orbits with computer assistance. 
Our sample system is the FitzHugh-Nagumo system given by 


u' = V 

v'= 5~^{cv — f {u)-\-w) (6.1) 

w' = ec~^(u — 7 U>), 


where a S ( 0 , 1 / 2 ), c, 7 and <5 are positive parameters and f{u) 
known as the system of traveling wave solutions ([/, W) = {ifu{x 
PDF: 

\Ut = C/.. + /([/) - W 
\Wt = e{U - 7 W) 


= u{u — a)(l 
-ct),ipw{x- 


— u). ( 6.11 is well- 
ct)) of the following 


( 6 . 2 ) 


6.1 Strategy and parameters 

By following arguments in previous sections, we validate global orbits. The following implemen¬ 
tation is basically common in our computations. In particular, we concentrate on construction of 
covering-exchange sequences associated with slow shadowing sequences. 
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Step 1. Fix eo > 0 and several parameters. 


Step 2. For constructing the j-th slow shadowing sequence (j = 0, • • • , p > 0) with 

TTy{N^’'^) = we set identical positive numbers ao,bo and h in (SS4) in advance. 


Before constructing each we compute approximate equilibria 0, for (6.1 )o 

(i.e., e = 0) along a (normally hyperbolic) branch of the nullcline. For simplicity, compute 
them with — Wj^i\ = h for all i so that is monotonously increasing (resp. 

decreasing) in the case g = -1-1 (resp. q = —1). 

4 -J.o 

Jj,eKit\ _ 


Also, let |?/+ - 


= H/2 for some H > 0 for further simplicity. Set 


7 ry(Afl’™‘) = + h,y+^ - d] in the case of g = -hi Similarly, set 7ry(A^®’^**) = -h 

(5, — h] in the case of g = — 1. Here d > 0 denotes an arbitrarily small number. For vali¬ 


dation of (SS5), we apply Proposition 4.32 Under these settings, verify the slow shadowing 
condition (4.5). 

In the case of validations of heteroclinic orbits, we additionally need to validate IF“(/e) for 
an isolating block A on slow manifolds. To this end, verify all assumptions in Lemma |4. 39 
We revisit this verification later in Lemma l6.11 


Step 3. If necessary, verify the unstable m-cone condition for an appropriate m > 1 in a block of 
the form ( 4.1l|) (Section [Fi] ). Also, construct the fast-exit face of unstable m-cones 

following (4.15[). Similarly, verify the stable m-cone condition of at the fast-entrance of blocks 


for an appropriate m > 1 if necessary. When we apply the predictor-corrector method for 


constructing fast-saddle-type block (Section 2.3.3), vertices of m-cones are slid in general. In 
this case, we cut the cone so that the fast-exit face is parallel to the y-axis. As a consequence, 
the length £u or £s of extended cones shorten at most lfx(x,y)~^fy(x,y)Hl, where H is the 
height of the block in the y-direction. See Fig. 

Step 4. Solve initial value problems of PDFs for setting a fast-exit-face of each block as an 

7-1,exit 


initial data and verify (CE4) in Definition 
operation consists of direct applications of 


4.2 


or Drop 




This 


interval arithmetics and ODE solver libraries such 


as CAPD [3]. Although our computations here are operated in full systems, slow dynamics 
can be regarded as the small error since our interest in this step is mainly the behavior of 
fast dynamics. 


Validation of assumptions in Lemma 4.39 in Step 2 can be done as follows. 


Lemma 6.1 (Validation of assumptions in Lemma 4.39). Let N be a fast-saddle-type bloek with 
'Ky{N) = [y)^,y)i)], and I be an isolating block I containing an equilibrium in the validated slow 
manifold with 0 < dist(S'e, V-^’”) < diam(7ra(V)) — where ra > 0, and dist(5'e, V-f’^) > ri, > 
0. Assume that N satisfies the unstable mu-cone condition with 


diam(7ro(V)) - 


ruu 


< n- 


Also assume that I can be chosen in 




mu 


(6.3) 


(6.4) 


and that any subset I G containing I is also an isolating block with the same isolating infor¬ 
mation as I. Then we can choose an isolating block C and the fast-exit face TV®’''* with 
7 ry(V®’''*) = [y",y+] and 




(6.5) 


so that assumptions in Lemma\4..3£\ holds. 
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(c) 


Figure 9: Cutting m-cones 


(a) : A fast-saddle-type block N with the basic form (Section 2.3.2} and associated unstable m- 
cone. The union of yellow and white regions represents N U C^. The yellow region represents 
{N U C'ff) D {y G Try{N)}. The x-components of block N are identical for y. 

(b) : A fast-saddle-type block with the predictor-correctorc form (Section 2.3 .k^ and associated 
unstable m-cone. The x-components of the block N as well as the cone C^ are slid, since the 
x-component of the center point also depends on y; (x,y) = (i -|- {dx/dy){y) ■ (y — y),y). The 
fast-exit face thus has non-trivial angle to the y-axis. 

(c) : Cut the edge of unstable m-cones. Resulting fast-exit face is parallel to the y-axis. The length 
iu of the extended cone shortens at most \fx{x,y)~^fy(x,y)H\, where H is the height of the block 
in the y-direction. 
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Proof. For eachp € I, any point q S W^{p) satisfies |7rh((7)—7rt,(p)| < (diam(7ra(-/V))—which 
is the consequence of properties of unstable cones. The same property holds for y-components. 
Therefore assumptions concerning with ( |6.3) an d (6.4) i mplie s (4.27). Similarly, it immediately 
holds that we can choose le and with (6.5) so that (4.28) is satisfied. □ 


In practical computations, the most difhcult part is Step 4. In general, the larger both the 
fast-exit face and the target h-set in Drop, i.e., _s, are, the easier validations of 

covering relations are. However, if we validate covering relations described in Step 4, the resulting 
h-sets is very close to the slow exit, i.e., + (resp. - ) in case that q = +1 

yj,i yj,i 

(resp. q = —1). In particular, the next fast-exit face becomes too thin to validate the next 

Drop. Procedures in Step 2 as well as Step 4 thus looks incompatible each other. Nevertheless, 
an appropriate choice of the slow shadowing ratio x avoids this inconsistency. 

Lemma 6.2 (Validity of Step 2). For each j = 1, • • • , p, consider the j-th slow shadowing sequence 
Assume that a covering relation holds. 

— ^3,0 

Let Xj be the slow shadowing ratio satisfying 


^ H / \wj^rnj - u!j,o\ - H\ H/mjh-H\ 

h ) ~ ~J h ) 


( 6 . 6 ) 


with mjh > H. 

Finally assume that is the slow shadowing sequence with the identical ratio Xj- Then 

we can choose a fast-exit face so that = [yj^. ~ (I'esp. 7ry(N™‘) = 

[yjmj + ^^y^rrij ~ ^]) casc of q = -\-l (rcsp. q = —1), where 6 > 0 is an arbitrarily small 

number. As a consequence, Step 2 and Step 4 are valid simultaneously. 


Proof. We only pr ove th e case g = +1. The case g = — 1 is similar. 
By Proposition 4.14 


pj.i 


with Xj, we can construct a covering relation where 

C and c 

The last equality follows from the choice of y^j^ and wj^i in Step 2. Repeating this argument, 

p3>'i--^ 

two h-sets describing the i-th covering relation Mp is located on (iV)'’*-i)p_^ and 

respectively, where 

h = y^j - {(* + 1 ) - iXj}h. 


We thus obtain = Vjm ~ {(™j + 1) ~ ^ the ratio Xj can be chosen satisfying (6.6), 

we obtain 


{{rrij + 1) - mjXjjh < < (m^ -h 1) - 1 - 


H - Wj.ol - H 


= 0 + 


rrijH (rridh — H 


-1' 


h. 


Consequently, 


ytm, - {{^1 + 1 ) - mjXj}^ < ylm, -h- 


m,i7 / m,7i — FI 


h 


-1 


< 2 /- 


-h- 


mjhH 
mjh — F[ 


< ytm. -h-H = y.^^ - h. 


Obviously, the slow shadowing pair with the ratio Xj satisfies the slow shadowing condition with 
the ratio x' for all x' € Therefore, arranging several Xj,i's, in the slow shadowing pair 
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we can take the h-set on . Statements of the lemma follow from 

the same arguments as the proof of Proposition |4.15[ Schematic pictures of the proof are shown 
in Fig. [T^ □ 


y 





(a) 




in 


h% 




(b) 


Figure 10: Schematic pictures of Lemma 6.2 (q = -1-1^ 

Dotted lines denote sections {Nf'^)y for each i. 

(a) : If the sequence of fast-saddle-type blocks is a slow shadowing sequence with the 

. . PJ’* . . 

ratio Xjj the section where the covering relation C holds gets 

lower than the usual version (xj = !)■ 

(b) : Repeating the procedure in (a) sufficiently many times, we can take the section {Ni’^^)y 


before Jump at the bottom of Nf ^; namely, {Ni 


As a consequence, we can take the 


fast-exit face large keeping the target h-set {Nf^)<:y in Drop large. 

Thanks to Lemma [6.2[ we replace Step 2 by the following, which enables us to verify assumptions 
of results in Section with large fast-exit faces and large target h-sets in Drop: 

Step 2’ Replace the statement “verify the slow shadowing condition” in Step 2 by “verify the 
slow shadowing condition with the slow shadowing ratio Xj” i where Xj < 1 is a given number 
satisfying (6.6). 


If Steps 1, 2’, 3 and 4 pass, then all assumptions of either Theorem |5.1[ Corollaries |5.3| or 
5.4 are satisfied in the case of periodic orbits for all e G (0, eg]. Similarly, all assumptions of either 
Theorem 15.61 Corollaries |5.7| or |5.8| are satisfied in the case of heteroclinic or homoclinic orbits. 


In Step 2, the predictor-corrector approach discussed in Section [2.3.3| is used for choosing local 
coordinates around (normally hyperbolic) invariant manifolds. A concrete form for validations is 
shown in Appendix [A} 

Note that all our examples below are cases u = 1. In such cases, o ne ca n directly verify covering 

relation ==A- by checking all assumptions in Proposition ! 

proposition to validating covering relations. 


2.9 


In Step 4, we apply this 


We gather parameters we deal with computations in our settings except ones which arise in 


(6.1), before moving to practical computation examples. These parameters are set for each branch 
of slow manifolds which we try to find. Let j G {0,1, • • • , p} be the number of branches. 


• h : The height in the w-direction for slow shadowing condition (4.5). A sequence of equilibria 
{{ui,Vi,Wi)} for (1.3) are set so that Iwi+i — Wi\ = h for all i. 


• H : The height of fast-saddle-type blocks in w-direction. Each w-interval of the length H 
corresponds to the set K in Section [2. 3. 2[ 


62 











































Kaname Matsue 


Rigorous Numerics for Fast-Slow Systems 


• wq : A given number such that jico — ii'ol < h, where wq is the w-component of an equilibrium 
{uo,vo,wo). The equilibrium (uo,vo,wo) is computed numerically and becomes the center of 
the fast-saddle-type block the target h-set of Drop. 

• Wruj '■ A given number such that Iwmj — Wm^ \ < h, where Wm, is the rc-component of 
an equilibrium (rtm,, Vmj, Wm ^) ■ The equilibrium {umj, f, Wm ^) is computed numerically 

jim - 

and becomes the center of the fast-saddle-type block ^; the last h-set of the j-th slow 
shadowing sequence containing a fast-exit face. 

• X ■ The slow shadowing ratio given by 


H f\wmj -wo\ - H - 2h\ 

X V h ) 


One can easily check that y satisfies ( 6 . 6 ). 


ra,ri, : The length of spaces in fast-saddle-type blocks introduced in (4.3). For simplicity, 
these numbers are identical for all blocks. Moreover, they are assumed to be identical each 
other. 


• da,dh : Positive numbers less than 1 introduced in (SS5). 


• m„,TOs : Positive numbers determining the sharpness of unstable and stable m-cones, re¬ 
spectively. 



All computations are done by MacBook Air 2011 model (1.6 GHz, Intel Core i5 Processor, 4GB 
Memory), GCC version 4.2.1 with -02 option and CAPD library [3] version 3.0. 


6.2 Demonstration 1 : slow shadowing sequences with the ratio y 


Validations of not only slow manifolds near critical manifolds consisting of equilibria for (1.3) 


but also the existence of trajectories which shadow slow manifolds are our starting points of whole 
considerations. Slow manifolds for (6.1) is now expected to be near the nullcline {u = 0, f{u) = w}. 


The aim of this subsection is to test how large slow manifolds can be validated in terms of slow 
shadowing sequence. Following Step 2 at the beginning of this section, we validate slow shadowing 
sequences. Note that validations in this section also verify Jump in Proposition 4.15| 

As a demonstration, we fix a = 0.3, 7 = 10.0, 5 = 9.0 and c G [0.799, 0.801]. These parameters 
are also used in Sections 16.71 and 16.81 


Computer Assisted Result 6.3. Consider (6.1) with a = 0.3, 7 = 10.0 and S = 9.0. Then 
for all c G [0.799,0.801] and e G [0,5.0 x 10“®], the branch of slow manifolds near the nullcline 
{v = 0, f{u) = w} in {-1.765629966434 x 10“® <u< 2.017612584956 x 10“®, -6.0 x 10“^ <w< 
0.099} is validated. In particular, the slow shadowing condition with q = —1 between blocks around 
this slow manifold is validated with parameters listed “First branch” in T'a&Ze[^ 


63 



























Kaname Matsue 


Rigorous Numerics for Fast-Slow Systems 


Similarly, the branch of slow manifolds near the nullcline {u = 0, f{u) = w} in {0.8504842978868 < 
u < 1.021440903396,-1.58 x 10“^ < w < 0.07} is validated. In particular, the slow shadowing 
condition with q = +1 between blocks around this slow manifold is validated with parameters listed 
“Third branch” in Table\^ 


Parameters 

First branch 

Third branch 

X 

0.8807339449541285 

0.8786764705882354 

h 

0.003 

0.003 

H 

0.0065 

0.0066 


0.75 

0.75 

db 

0.7 

0.75 

ra 

0.008 

0.008 

Tb 

0.0085 

0.0078 

mu 

100 

100 

ms 

100 

100 

computation time 

0.566 sec. 

0.467 sec. 


Table 1: Validation parameters of slow shadowing in Computer Assisted Result [6731 


This validation result implies that we have already validated trajectories with appropriately 
chosen initial data, say h-sets, which shadow slow manifolds with an arbitrary length for all 
e e (0,5.0 X 10“^], as long as slow shadowing are validated. Moreover, Jump has been also 
validated for any fast-exit face with an appropriate height from the bottom (for q = -fl). The 
height H can be explicitly determined from ( 6.61 in our setting. Remark that the range of our 
validating slow manifolds is not the limit of our verifications. Validations of slow shadowing 
sequences are just iterations of Step 2 and can be validated very fast, if we have fast solver of 
linear algebra. Notice that slow shadowing sequences for e G (0,eo] validate trajectories which 
shadow slow manifolds without solving any differential equations for all e G ( 0 , eo]- 

On the other hand, there is a trade-off for validating slow shadowing sequences. For example, 
if we raise the value of e, say 6.0 x 10“®, the slow shadowing condition (4.5) violates, since the 
slow speed becomes faster than expansion and contraction of h-sets in hyperbolic directions around 
slow manifolds. Factors determining (4.5) are h, H, da, d&, ra, r^ as well as eigenvalues and size 
of fast-saddle-type blocks. One expects that, the larger parameters, say h, ra, ri, are, the easier 
the validation of (|4.5|) will be. However, in such a case, the covering relation in Proposition 4.32 


is often violated. In particular. Assumption (SS5) is violated. This is mainly because the distance 
between two centers {ui,Vi, wf) and (wi+i, Ui+i, Wi+i) become larger and hence the affine map ra;,i 2 
moves h-sets larger, if we increase h, ra, ri,. This is also the case if we increase parameters da and 
db- 


6.3 Demonstration 2 : m-cones 

Next, we show a demonstration of m-cones. When we want to construct a covering-exchange 

sequence, we need to verify covering relation Tf ^ Se, where T} is a fast-exit face of a fast- 

saddle type block and is an other block. In order to verify this covering relation, we solve 
differential equation with the initial data J-f. On the other hand, we can replace iFf by a fast-exit 
face J-f of an extended cone, thanks to discussions in Section 
initial data, Ff and Ff, to see the following two points: 

1 . accuracy of solution orbits, and 

2 . verification of covering relations. 

As an example, we set a = 0.01, 7 = 0.0, <5 = 5.0, c G [0.495,0.505] and e G [0.0,5.0 x 10“®]. These 
are parameter values used in Section[6.4| In demonstrating computations, we used the ODE solver 


4.4 


Here we solve (6.1) with two 
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(c) (d) (e) 

Figure 11: Comparison of solution enclosures with and without m-cones. 

Horizontal axis: u. Vertical axis: v. Each figure represents the projection of trajectories and cones 
on {u, v)-plane. 

(a) Validation of solution orbits (green) with initial data Ff. Pink regions are the union of fast- 
saddle-type blocks and extended cones. Computation step of ODEs is 590. 


(b) The same computation result as (a). In this case the stable m-cone is not validated. The tiny 
pink region around {u,v) = (0.956721,0) is the validated fast-saddle-type block. Readers see that 
enclosures of trajectories are much bigger than the fast-saddle-type block. 

(c) Validation of solution orbits (red) with initial data F) (i.e. without unstable m-cones) and the 
same time steps as (a). Validated trajectories do not arrive at target region yet. 

(d) Validated trajectories with additional time step computations to (c). More precisely, computa¬ 
tion step of ODEs is 620. Enclosures become bigger and bigger. 


(e) Validated trajectories with additional time step computations to (d). More precisely, compu¬ 
tation step of ODEs is 650. Enclosures are already bigger than stable cones, which implies that 
validation of covering relations can be never done. 
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in CAPD library based on Lohner’s method discussed in [35]. The order of Taylor expansion is set 
p = 6 and time step size is set At = 0.0001. Computation result is shown in Fig. [m 


First we compare the case we use stable m-cones with the case we do not use stable m-cones. 


Fig. 11 (a) and (b) show the same computational results. The only difference is whether or not 
a stable m-cone is validated. In this example, we validate the stable 29-cone with the length 
£s = 0.108232 of a fast-saddle-type block around (0.956575,0,0.0392) G (see Fig. 11 (a)). We 
solved ODE with initial data 


n {w G [0.039249948844,0.03926994850296]} 

after dividing it into uniform 30 small pieces. Here denotes the fast-exit face of unstable 21- 
cone with the length £„ = 0.0247787 of a fast-saddle type block around the origin in K^. In this 
example, we solved ODE in 590 steps. 

In general, validated fast-saddle type blocks corresponding to are very small, as shown in 
Fig. [^(b) (colored by pink). On the other hand, validated trajectories are quite bigger than 
blocks. In our example, validated trajectories are already bigger than the block, which implies 

that we can never validate ') setting. A direct settlement of this problem would 

be a refinement of initial data, which leads to huge computational costs in many cases and is not 
realistic. Instead, we consider the problem with the help of stable m-cones, which is shown in Fig. 
11 (a). In this case, the target block corresponding to becomes big enough to validate covering 
relations. Thanks to Section |4.4[ we can discuss validation of trajectories with extended cones, 
which is much easier than verifications without cones. 


Next, we compare the case which we use unstable m-cones with the case which we do not use 
unstable m-cones. Computation result with unstable m-cone is Fig. 11 (a). Figs. 11 (c), (d) and 
(e) show enclosure of trajectories with initial data namely, a fast-exit face of a small fast- 
saddle-type block. If we do not use unstable m-cones, we need to solve ODEs for longer time steps 
than the case we use unstable m-cones (Fig. 11 (c)). Such extra computations cause additional 
computational errors and there is little hope to validate covering relations, as indicated in Figs. 


11 (d), (e). 


6.4 Periodic orbits 


We go to validations of global orbits for (6.1 1 . Our first example is validation of periodic orbits. 
As a demonstration we set a = 0.01, 7 = 0.0 and (5 = 5.0 and c « 0.5. All validation of covering- 
exchange sequences with extended cones yield the following computer assisted result. 


Computer Assisted Result 6.4. Consider (6.1) with a = 0.01, 7 = 0.0 and S = 5.0. Then for 
all cG [0.495,0.505] the following trajectories are validated. 


1. At e = 0, there is a singular heteroclinic chain Hq consisting of 


• heteroclinic orbits from pq to qq, and from pi to qi, 

• branches M°, of nullcline {u = 0 ,/(m) = wj. contains po and go. Similarly, 

contains pi and qi. 
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Equilibria po, go, pi and qi are validated by 

Uu,v{Po) - (-0.177098234,2.18166218 x lO"®)] < 2.49103628 x lO'^, 

\-Ky{po) - 0.03951 < 3.25 x 10-^ 

| 7 r„.„(gi) - (0.956125336,-5.98704406 x 10"®)| < 2.10571434 x 10■^ 

\TTy{qi) - 0.039321 < 3.3 x 10-^ 

Uu,v{Pi) - (0.850811351,-9.47162333 x lO"®)] < 2.44899980 x 10-^ 

\-Ky{pi) - 0.106021 < 3.3 X 10-^ 

k„,^(<?o) - (-0.282970990,1.85607735 x lO"®)! < 2.12552591 x 10 -^ 

\TTy{qo) - 0.106751 < 3.25 x 10"^. 

2. For all e G (0,5.0 x 10“®], there exists a periodic orbit near Hq. 

Parameters for validations are listed in Table 

We omit computation times for slow shadowing and Jump, since they take only a few seconds 
as stated in Section [R2] 


Parameters 

On {q = -1) 

On ]VE (g = +1) 

X 

0.8922056384742952 

0.8895212587880817 

h 

0.0025 

0.0023 

H 

0.0065 

0.0066 


0.8 

0.8 

db 

0.7 

0.8 

Ta 

0.008 

0.008 

n 

0.0085 

0.008 

rUu 

21 (around po) 

55 (around pi) 

uis 

28 (around go) 

29 (around gi) 

4 

0.0247787 (around po) 

0.0186724 (around Pi) 

4 

0.113375 (around go) 

0.108232 (around gi) 


jy0,exit y’e ^1,0 

_^l,exit ^0,0 

T, 

0.059 (with At = 1.0 x 10“^) 

0.0266 (with At = 2.0 x 10“®) 

computation time 

70 min. 6 sec. 

78 min. 53 sec. 


Table 2: Validation parameters of slow shadowing in Computer Assisted Result 



6.5 Heteroclinic cycles 

The second example is a family of heteroclinic cycles. Consider the cubic curve w = f(u). One 
can see that it is symmetric with respect to the inflection point 


(^inf : ^inf ) 


/1 -I- a (1-1- a)(l — 2a)(2 — a) 
V 3 ’ 27 


and we can choose (ua, ws) from the curve w = f{u) to be the point that is symmetric to the origin 
{u,w) = (0,0) with respect to the inflection point. Deng [5] shows that 70 := 9/(2 — a)(l — 2a) 
and Co := (1 — 2 a )/-\/2 with <5 = 1.0 admit a heteroclinic loop of ( 0 , 0 ) and ( 1 x 3 , 103 ) for sufficiently 
small e. Symmetry of the cubic curve w = f{u) with respect to (uinf, rcinf) implies that the vector 


field ( 6 . 1 ) with ( 7 , 0 ) = ( 70 , 09 ) is symmetric under the following transformation: 


2(1 -|- a) 

u- u, 


~ - 2 (l + a)(l- 2 a)( 2 -a) 
w - - - - -w. 
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Figure 12: Validation of covering-exchange sequences in Computer Assisted Result \K^ 
Horizontal axis: u. Vertical axis: v. Each figure represents the projection of trajectories and cones 
on (u, v)-plane. 

(a-1) : Validation o/(C“ J®”*) n {w S [0.0365499948149,0.0422098969893467]}. 

(a-2) : Validation o/(/5,(ro, (C" n {w e [0.0365499948149,0.0365699944692267]}. 

(a-3) : Validation o/(C“ n {w G [0.04218989733502,0.0422098969893467]}. 

(b-1) : Validation o/(p,(Ti, (C" J™*) n {w G [0.1037942789937,0.10913580929979]}. 

(b-2) : Validation o/(p,(Ti, (C” J™*) n {w G [0.1037942789937,0.103814744627057]}. 

(b-3) : Validation o/(p,(Ti, (C" J™*) n {w G [0.109115343666433,0.10913580929979]}. 
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It is thus sufficient to validate a heteroclinic orbit from the origin to ( 143 , 0 ,^ 3 ) if we want to 
validate heteroclinic cycles. Here we validate heteroclinic cycles for concrete e and specific a and 
6. All validations of covering-exchange sequence with extended cones yield the following computer 
assisted result. 


Computer Assisted Result 6.5. Consider (6.1 ) with a = 0.3, 7 = 70 , c = cq. Then, ford = 1.0, 
the following trajectories are validated. 


1. At € = 0, there is a singular heteroclinic chain Hq consisting of 

• heteroclinic orbits from Pq to qi, and from pi to qo, 

• branches M°, of nullcline {v = 0,f{u) = w}. contains po and qo. Similarly, 

contains pi and qi. 

Equilibria po and qi are validated by 

Uu,v{Po) - (3.34681997 x 10-'‘,-1.97213258 x lO"®)] < 1.35013586 x 10”®, 
k!/(Po)l < 1-5 X 10"^, 

Uu,v{<li) - (0.999441149,-1.32078036 x lO"^)] < 1.89700763 x 10"®, 

\TTy{qi) - 2.0 X 10-'‘| < 2.5 x 10"^ 

The equilibrium po admits an attracting isolating block on contained in {(u,v,w) G 
{u < 1.80738697 x lO"®} | w e [-1.200 x 10-^ 6.010 x IQ-^]}. 

Note that Pi is symmetric to po with respect to the inflection point (rtinf, 0, Winf). Similarly, 
qo is symmetric to qi with respect to the inflection point (rtinf, 0 , Winf). 

2. For all e G (0, 5.0 x 10“®], there exists a heteroclinic cycle near Ho. The cycle consists 
of heteroclinic orbits from po to qi and from pi to qo. 

Parameters for validations are listed in Table^ Finally, assumptions of Lemma \6.1\ are validated 
with the attracting isolating block {ru G [—1.208036756757 x 10“®, 1.208036756757 x 10“®]} on the 
slow manifold containing po and the unstable 37-cone condition. 


Parameters 

On {q = -1) 

On {q = +1) 

X 

0.9358974358974359 

0.9168053244592347 

h 

0.0001 

0.0001 

H 

0.003 

0.005 

da 

0.8 

0.8 

db 

0.7 

0.8 

ra 

0.003 

0.006 

n 

0.0045 

0.008 

mu 

14 (around po) 

- 

TOs 

- 

2.2 (around qi) 

2 (diam( 7 ra(A)) - ra)lmu 

4.271697297297 x 10“^ 

- 

4 

0.0151974 (around po) 

- 

4 

- 

0.0949713 (around qi) 


jyO,exit ^ 1,0 

jyl,ex.it jyO.O 

T. 

0.045 (with At = 1.0 X 10“^) 

- 

computation time 

43 min. 51 sec. 

- 


Table 3: Validation parameters of slow shadowing and isolating blocks in Computer Assisted Result 
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(a-1) 


Figure 13: Validation of covering-exchange sequences in Computer Assisted Result \C^ 
Horizontal axis: u. Vertical axis: v. Each figure represents the projection of trajectories and cones 
on (u, v)-plane. 

(a-1) : Validation ofipfiTo, (C” S [-0.00103331528021667,9.35989145343001 x 10-^]}. 

(a-2) : Validation o/(C““ n {w G [-0.00103331528021667,-0.001023982083887]}. 
(a-3) : Validation ofipfiTo, (C" G [9.26655949013333x 10-^ 9.35989145343001x 10-^]}. 


6.6 Heteroclinic orbits 


The third example is a family of heteroclinic orbits which are different from heteroclinic cycles. 

In this example we apply the method in Section [4.5| to isolating blocks on slow manifolds, which 
validates equilibria on slow manifolds. 

Let {u,0,w) zfi, (0,0,0) be a point in such that fifa) = w and that u = jw. Assume that 
fu(u) 7 ^ 0. Then there is a function u = h{w) which is unique in a small neighborhood of (u, 0, w) 
such that u = h{w), f(h{w)) = w and f{h{w)) = w hold in such a neighborhood. This is due to 
the Implicit Function Theorem. Assuming that fu(u) 0 holds for all m in a given neighborhood 
of u, the above implicit function representation holds in the given neighborhood. By using implicit 
function differential, the vector field g{u,v,w) = c~^{u — jw) is rewritten by 


g{u,v,w) 


- 

c \fu{u) 




u) 


w — w\ 

fu{u) ) 


near {h{w),0,'w) G Remark that u = jw. This expression leads to an effective estimate of 
vector fields on slow manifolds, as stated in Section [TSl 

All validation of slow shadowing sequences with extended cones and isolating blocks on slow 
manifolds yield the following computer assisted result. 


Computer Assisted Result 6.6. 


Consider (6.1) with a 


0.2, 7= 15.0, 5 = 5.0. 


1. For each c G [0.947,0.948], the following two kinds of trajectories are validated: 


(a) At e = 0, there is a singular heteroclinic chain consisting of a heteroclinic orbit from 
Po = (0,0,0) to qi and a branch of nullcline {w = 0,/(u) = uij connecting qi and 
Pi = Equilibria have the following estimates: 


huAPo) - (8.00733430 x 10-^,2.18164458 x lO"^)] < 2.77714227 x 10-^ 
liTyipo) - 3.3 X lO-'^l < 2.75 X 10"^, 

k„,„(gi) - (1.00015243,-1.21520324 x lO"®)! < 2.11148288 x 10-^ 
iTTyiqi) + 4.6 X 10-^1 < 3.0 x Wk 

K.viPi) - (0.882676299,6.72075379 x lO-^j < 2.64847328 x Wk 
l7rj,(pi) - 0.0596] < 3.0 x 10"^ 


(b) For all e G (0,1.0 X 10 ^], there exists a heteroclinic orbit H) from pq to pi near Hq. 
The equilibrium Pi admits an attracting isolating block on contained in {(u,v,w) G 
S', n {u > 0.848998902} ] w G [5.9080 x Wk 5.9988 x IQ-^]}. 
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Parameters for validations are listed in Table Finally, assumptions of Lemma \6.1\ are 
validated with the attracting isolating block {w € [—2.11442375 x 10“^, 2.77442375 x 10“^]} 
on the slow manifold containing po and the unstable 80-cone condition. 

2. For each c G [0.1995,0.2005], there exist the following two kinds of trajectories: 

(a) At e = 0, there is a singular heteroclinic chain 77q consisting of a heteroclinic orbit 
from Pi to qq and a branch of nullcline {u = 0,f{u) = ru} connecting qo and pq. 
Equilibria have the following estimates: 

KApi) - (0.876694738,2.76020695 x lO”®)! < 2.69669311 x 10-^ 
iTTyipi) - 0.06181 < 4.0 X 10-^ 

k„,„(go) - (-0.128335342,1.50674261 x 10"®)| < 1.52471761 x 10■^ 

|7rj^(qo) - 0.06221 < 4.0 x 10-^ 

k„,„(po) - (-2.67285871 x 10-^ 1.20922437 x lO"®)! < 2.32761788 x 10"^ 
iTTyipo) + 1.2 X 10-®| < 4.0 X 10"®. 


(b) For all e G (0,1.0 x 10 ®], there exists a heteroclinic orbit Hf from pi to po near TJq . 
The equilibrium Po admits an attracting isolating block on contained in {(u,v,w) G 
S', n {u < 1.06035231 x IQ-^} | w G [-1.040 x IQ-^, 1.264 x lO"®]}. 


Parameters for validations are listed in Table ^ Finally, assumptions of Lemma 6.1 are 
validated with the attracting isolating block {w G [5.800831058824 x 10“^,6.51916894117 x 
10“^]} on the slow manifold containing pi and the unstable 85-cone condition. 


In this case we have to care about the existence of an equilibrium on the slow manifold in Si. 
Easy calculations yield that Si possesses at most one equilibrium. 


Parameters 

On {q = —1) 

On Afi {q = -hi) 

X 

0.8949178448605274 

0.8877665544332212 

h 

0.00023 

0.00022 

H 

0.0055 

0.006 

da 

0.8 

0.75 

db 

0.8 

0.8 

Ta 

0.008 

0.008 

n 

0.0085 

0.008 

m„ 

50 (around po) 

- 

rus 

- 

25.5 (around gi) 

2(diam(7ro(lV)) - ra)/m„ 

4.169025 X 10-4 

- 

iu 

0.0136582 (around po) 

- 

4 

- 

0.117888 (around qi) 


ex.it jyljO 

jYl,exit ^0,0 

T, 

0.064 (with At = 1.0 X 10-4) 

- 

computation time 

38 min. 21 sec. 

- 


Table 4: Validation parameters of slow shadowing and isolating blocks in Computer Assisted Result 


6.6 with cG [0.947,0.948]. 
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Figure 14'- Validation of covering-exchange sequences in Computer Assisted Result \67E[ 
Horizontal axis: u. Vertical axis: v. Each figure represents the projection of trajectories and cones 
on (u, v)-plane. 

(a-1) : Validation o/(C“ J®”*) n {w S [-0.00142875124866067,0.00177996171276667]}. 
(a-2) : Validation o/(/5,(ro, (C" n {w e [-0.00142875124866067,-0.00141195170436]}. 

(a-3) : Validation o/(C“ fl {w G [0.001763162168466,0.00177996171276667]}. 

(b-1) : Validation o/(p,(Ti, (C" J™*) n {w G [0.058532027216316,0.06509312870908]}. 

(b-2) : Validation o/(p,(Ti, (C” J™*) n {w G [0.058532027216316,0.058557359654512]}. 

(b-3) : Validation o/(p,(Ti, (C" J™*) n {w G [0.065067796270884,0.06509312870908]}. 
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Parameters 

On {q = -1) 

On {q = -71) 

X 

0.8484848484848486 

0.8540145985401461 

h 

0.0001 

0.0002 

H 

0.008 

0.008 

da 

0.8 

0.75 

db 

0.8 

0.8 

Ta 

0.0042 

0.006 

n 

0.0045 

0.0068 

m„ 

- 

20 (around p{) 

rus 

26.5 (around qo) 

- 

2 (diam(7ro(A)) - ra)/m„ 

- 

4.083105882353 x 10-^ 

In 

- 

0.0539878 (around pi) 

Is 

0.113106 (around qo) 

- 


jY0,exit jyl,0 

jyl,exit r) ^0,0 

T, 

- 

0.055 (with At = 1.0 X 10-4) 

computation time 

- 

39 min. 41 sec. 


Table 5: Validation parameters of slow shadowing and isolating blocks in Computer Assisted Result 


6.6 with ce [0.1995,0.2005], 


6.7 Homoclinic orbits 


The final example is a family of homoclinic orbits. Set a = 0.3, 7 = 10.0, 5 = 9.0 and c « 0.8. 
All validation of covering-exchange sequences with extended cones yield the following computer 
assisted result. 


Computer Assisted Result 6.7. Consider (6.1) with a = 0.3, 7 = 10.0 and S = 9.0. Then for 
all cS [0.799,0.801] the following trajectories are validated. 


1. At € = 0, there is a singular heteroclinic chain Hq consisting of 


• heteroclinic orbits from Pq to qi, and from pi to Qq, 

• branches M°, of nullcline {v = 0,/(m) = w}. contains po and go. Similarly, 

contains qi and pi. 


Equilibria Po, qi, pi and qo are validated by 

|7r„,^(po) - (-4.22605959 x 10-^4.67157371 x lO"'^)! < 2.90338686 x 10-^ 
liTyipo) - 1.94 X 10-3| < 3.25 X 10-^ 

| 7 r„.„(( 7 i) - (0.996532931,-4.66925590 x lO"®)! < 2.13488628 x lO'^, 

\TTy{qi) - 2.02 X 10-^1 < 3.3 x 10-^ 

| 7 r„,^(pi) - (0.870020061,-2.91627199 x 10"®)| < 2.85183938 x 10■^ 
\-Ky{pi) - 0.063621 < 3.3 x 10-^ 

k„.„(<?o) - (-0.129665170,3.70642116 x lO"®)! < 2.20129903 x 10-^ 

|7ry(go) - 0.063351 < 3.25 x 10"^ 


2. For all e G (0,5.0 x 10 ®], there exists a homoclinic orbit 77^ ofpo near Ho. The equilibrium 
Po admits an attracting isolating block on contained in {{u,v,w) G S'en{M < 2.10010357 x 
10-2} I w G [-3.1550 X 10-^ 2.1545 x IQ-^]}. 

Parameters for validations are listed in Table^ Finally, assumptions of Lemma \6.1\ are validated 
with the attracting isolating block {w G [—1.09722 x 10-^,4.97722 x IQ-^]} on the slow manifold 
and the unstable 120-cone condition. 
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Parameters 

On {q = -1) 

On {q = +1) 

X 

0.8812568505663135 

0.8792535675082328 

h 

0.00023 

0.00022 

H 

0.0065 

0.0066 

da 

0.8 

0.8 

db 

0.75 

0.8 

ra 

0.008 

0.008 

n 

0.0085 

0.0078 

Ulu 

50 (around po) 

52 (around pi) 

iris 

45 (around qo) 

55 (around gi) 

2 (diam( 7 ra(V)) - ra)lmu 

2.9606 x 10“^ 

- 

4 

0.0441793 (around po) 

0.0370798 (around pi) 

Is 

0.12217 (around go) 

0.0984568 (around gi) 


jyO,exit jyl,0 

^l,exit ^ 0,0 

D 

0.067 (with At = 1.0 X 10“4) 

0.0288 (with At = 2.0 x 10“®) 

computation time 

61 min. 59 sec. 

73 min. 57 sec. 


Table 6 : Validation parameters of slow shadowing and isolating blocks in Computer Assisted Result 

o 


Remark 6 . 8 . In this case we easily know that the only equilibrium of (6.1 1 for e > 0 is (u, v, w) = 
( 0 , 0 , 0 ) and the dynamics around ( 0 , 0 , 0 ) can be easily understood to show that ( 0 , 0 , 0 ) is attract¬ 
ing on the slow manifold in 5°. 

Note that appropriate arrangements of computation schemes can improve the accuracy of val¬ 
idated trajectories and the validation range of parameters including e. 

In many articles, the uniqueness of global orbits is also discussed. From the geometrical view¬ 
point, typical arguments for the uniqueness require the transversality of locally invariant manifolds. 
In our cases, however, transversality in the sense of differential manifolds is not mentioned. In other 
words, we cannot prove the uniqueness of validated orbits in the current setting. 


6.8 Continuation of homoclinic orbits 


Finally, we verify the e-continuation of homoclinic orbits validated in the previous subsection. It is 
very hard to validate trajectories of (6.1 )e for all e S (0, eo] with large eo at a time. In most cases. 


covering relations in the fast dynamics fail. To overcome this difficulty, we divide a closed interval 
[ 0 , eo] into sub-intervals 


[0,eo]U[eJ,eg]U---U 


mi m _ 

1 ^0 J1 ^0 — £o 


and validate trajectories within each subinterval. 

In our demonstrating example, we divide the interval [0, 5.0 x 10“®] into 

[0,5.0 X 10“®] U [5.0 X 10“®, 1.2 x 10“®] U [1.2 x 10“®, 2.2 x 10“®] 

U [2.2 X 10“®,3.2 X 10“®] U [3.2 x 10“®,4.2 x 10“®] U [4.2 x 10“®,5.0 x 10“®]. 


Validating homoclinic orbits for each sub-interval and summarizing them, we obtain the fol¬ 
lowing. Let po, go,Pi and qi be points corresponding to those stated in Computer Assisted Result 

[RTl 


Computer Assisted Result 6.9. Consider (6.1) with a = 0.3, 7 = 10.0 and S = 9.0. Then for 
all c S [0.799,0.801] and for all e S (0,5.0 x 10“®], there exists a homoclinic orbit of po near 
Hq, where Hq is the singular homoclinic orbit for (6. Do obtained in Computer Assisted Result 6.7 
Several numerical data for validations of iLg are listed in Table 
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Figure 15: Validation of covering-exchange sequences in Computer Assisted Result^d^ 
Horizontal axis: u. Vertical axis: v. Each figure represents the projection of trajectories and cones 
on (u, v)-plane. 

(a-1) : Validation o/(C“ J®”*) n {w S [-0.001128757897204,0.00412596711158]}. 

(a-2) : Validation o/(/5,(ro, (C" n {w e [-0.001128757897204,-0.00110862485119333]}. 

(a-3) : Validation o/ (C“ fl {w G [0.00410583406556933,0.00412596711158]}. 

(b-1) : Validation o/(p,(Ti, (C" J™*) n {w G [0.06094493968776,0.0668173357334667]}. 

(b-2) : Validation o/(p,(Ti, (C” J™*) n {w G [0.06094493968776,0.0609654725410667]}. 

(b-3) : Validation o/(p,(Ti, (C" J™*) n {w G [0.06679680288016,0.0668173357334667]}. 
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e 

^mn 

h 

H 

X 

f 

[0.0,5.0 X 10-“] 

0.00175 

0.00023 

0.0065 

0.8812568505663135 

0.04417933 

[5.0 X 10"®, 1.2 X 10-5] 

0.00175 

0.00023 

0.0065 

0.8812568505663135 

0.0441793 

[1.2 X 10-5,2.2 X 10-5] 

0.00175 

0.00023 

0.0065 

0.8813652126300421 

0.0441725 

[2.2 X 10-5,3.2 X 10-5] 

0.00175 

0.00023 

0.0065 

0.8813652126300421 

0.0441528 

[3.2 X 10-5,4.2 X 10-5] 

0.00175 

0.00023 

0.0065 

0.8813652126300421 

0.0441429 

[4.2 X 10-5,5.0 X 10-5] 

0.0019 

0.00025 

0.0065 

0.8809523809523810 

0.0441474 


< h. 


Table 7: Main information of blocks around po = with 

The sharpness of unstable cones is an identical value : = 50. The ratio da,db in (4.5) are 

set as identical values : da = 0.8, db = 0.75. The space length parameters Va, Vb of fast-saddle-type 
blocks are set as identical values : ra = 0.008, r;, = 0.0085. 


e 


h 

H 

X 

4 

[0.0,5.0 X 10-«] 

0.0633 

0.00023 

0.0065 

0.8812568505663135 

0.12217 

[5.0 X 10-®,1.2 X 10-^] 

0.0633 

0.00023 

0.0065 

0.8812568505663135 

0.12217 

[1.2 X 10-^2.2 X 10-^] 

0.0635 

0.00023 

0.0065 

0.8813652126300421 

0.122164 

[2.2 X 10-^3.2 X 10-^] 

0.0635 

0.00023 

0.0065 

0.8813652126300421 

0.122145 

[3.2 X 10-^4.2 X 10-^] 

0.0635 

0.00023 

0.0065 

0.8813652126300421 

0.122141 

[4.2 X 10-5,5.0 X 10-5] 

0.0635 

0.00025 

0.0065 

0.8809523809523810 

0.122144 


Table 8: Main information of blocks around qo = {uq,Vq,Wq) with \w}j — ihol < h. 

The sharpness mg of stable cones is an identical value : mg = 45. The ratio da, db in (4.5) are set 
as identical values : da = 0.8, db = 0.75. The space length parameters ra,rb of fast-saddle-type 
blocks are set as identical values : = 0.008, r;, = 0.0085. 


Conclusion 

We have proposed a topological methodology for validating singularly perturbed periodic, homo¬ 
clinic and heteroclinic orbits for fast-slow systems with its applicability to concrete systems with 
computer assistance. The main features of our proposing methodology are the following. 

• Our central strategy consists of applications of well-known topological tools, covering rela¬ 
tions, cones and isolating blocks, with taking the singular perturbation structure of normally 
hyperbolic invariant manifolds into account. This is one of different points from preceding 
works such as [1]. 

• All our procedures can be validated via rigorous numerics (e.g. interval arithmetics). As a 
consequence, one can validate the continuation of singular limit orbits for all e € (0,eo] with 
a given Eq > 0 with computer assistance. 

A main concept for realizing the above points simultaneously is the covering-exchange; the 
singular perturbation’s version of covering relations. We have proposed not only its primitive 
form but also its generalization; a collection of local behavior named “slow shadowing”, “drop” 
and “jump”. The notion of covering-exchange describes the behavior of trajectories in the full 
system which shadow normally hyperbolic invariant manifolds as well as their stable and unstable 
manifolds, even for sufficiently small e. In particular, this notion enables us to validate such 
trajectories without solving any differential equations in practical computations. Moreover, as 
mentioned above, it exceeds the limit of the multiple timescale parameter range: “sufficiently 
small e”, to e with an explicit range in practical applications with computer assistance. 

Of course, further applications of topological tools such as covering relations and the Conley 
index enable us to prove the existence of trajectories with complicated behavior. We believe that 
all these ideas will overcome difficulties of a broad class of singular perturbation problems. 

We end this paper proposing further directions of our arguments. 
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e 

wL, 

h 

H 

X 

P 

[0.0,5.0 X 10“®] 

0.0636 

0.00022 

0.0066 

0.8792535675082328 

0.0370798 

[5.0 X 10“®,1.2 X 10“®] 

0.0636 

0.00022 

0.0066 

0.8792535675082328 

0.0370798 

[1.2 X 10“®,2.2 X 10“®] 

0.0636 

0.00022 

0.0066 

0.8792535675082328 

0.0370719 

[2.2 X 10“®,3.2 X 10“®] 

0.0636 

0.00022 

0.0066 

0.8792535675082328 

0.0370534 

[3.2 X 10“®,4.2 X 10“®] 

0.0636 

0.00022 

0.0066 

0.8792535675082328 

0.0370464 

[4.2 X 10“®,5.0 X 10“®] 

0.0634 

0.0003 

0.0065 

0.8786764705882354 

0.0370888 


Table 9: Main information of blocks around pi = (u^ , ) with 


mi ’ 

The sharpness mu of unstable cones is an identical value : m„ = 52. The ratio da,db in (4.5) 
are set as identical values : da = 0.8, db = 0.8 for e S (0.0,4.2 x 10“^], da = 0.75, db = 0.75 for 
e G [4.2 X 10“^, 5.0 x 10“®]. The space length parameters ra,rb of fast-saddle-type blocks are set 
as identical values : ra = 0.008, r;, = 0.0078. 


I 

I ^mi ^mi I 


< h. 


e 

Wq 

h 

H 

X 

Is 

[0.0,5.0 X 10“®] 

0.0019 

0.00022 

0.0066 

0.8792535675082328 

0.0984568 

[5.0 X 10“®,1.2 X 10“®] 

0.0019 

0.00022 

0.0066 

0.8792535675082328 

0.0984568 

[1.2 X 10“®,2.2 X 10“®] 

0.0019 

0.00022 

0.0066 

0.8792535675082328 

0.0984452 

[2.2 X 10“®,3.2 X 10“®] 

0.0019 

0.00022 

0.0066 

0.8792535675082328 

0.098424 

[3.2 X 10“®,4.2 X 10“®] 

0.0019 

0.00022 

0.0066 

0.8792535675082328 

0.0984201 

[4.2 X 10“®,5.0 X 10“®] 

0.0018 

0.0003 

0.0065 

0.8786764705882354 

0.0984281 


Table 10: Main information of blocks around qi = (ug,UQ, Wq) with 






< h. 


The sharpness of stable cones is an identical value : = 55. The ratio da,db in (4.5) 

set as identical values : da = 0.8, db = 0.8 for e G (0.0,4.2 x 10“"^], da = 0.75, df, = 0.75 
e G [4.2 X 10““^, 5.0 x 10“®]. The space length parameters ra,rb of fast-saddle-type blocks are 
as identical values : = 0.008, r;, = 0.0078. 


are 

for 

set 


Bridging validation results to “direct” approach 


Our methodology concentrates on validations of trajectories for e G (0, eo] on the basis of geometric 
singular perturbation theory. One of challenges concerning with continuation of trajectories is the 
integration of singularly perturbed trajectories with those by a direct approach; namely, validation 
of trajectories without taking account of singular perturbation structure of systems for positive e 
bounded away from 0 (e.g. mi In our validation examples, under a specific choice of parameter 
values, solution trajectories has been validated for e G (0, 5.0 x 10“®] in the case of homoclinic 


orbits for the FitzHugh-Nagumo system (6.1). In the case of (6.1), validations of each branch of 


slow manifolds are completed in less than one second. 

This parameter range is still far from validation ranges with direct approaches; e = 0.01 in [T], 
for example. Validations of slow shadowing for further range of e are not easy, which is mainly 


because there is a trade-off between the slow shadowing condition (4.5) and the covering relation 


in Assumption (SS5), as discussed in Section 6.2 Both conditions are essential to describing slow 
shadowing, and hence our task for overcoming this situation will be developments of an improved 
topological and numerical method, like multi-step methods in numerical initial value problems, 
keeping the essence of slow shadowing phenomena. Obviously, if one can validate branches of 
normally hyperbolic slow manifolds by one blocks, there is no problem in this direction, in which 
case we can apply the primitive form of covering-exchange. After improvements of validations 
of slow shadowing or blocks containing slow manifolds with computer assistance, there will be a 
possibility that validated trajectories can be further continued in e and, we hope, be connected 
to validated ones via direct approach under appropriate choice of various parameters or solvers. 
Further works in this direction from the viewpoint of topological and numerical analysis are ones 
of our future studies. 


Uniqueness and stability of validated orbits 
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e 

^0,exit ^1,0 

^l,exit at 0,0 

iVg ^ 

[0.0,5.0 X 10“®] 

61 min. 59 sec. 

73 min. 57 sec. 

[5.0 X 10“®,1.2 X 10“®] 

61 min. 27 sec. 

75 min. 57 sec. 

[1.2 X 10“®,2.2 X 10“®] 

61 min. 38 sec. 

74 min. 55 sec. 

[2.2 X 10“®,3.2 X 10“®] 

59 min. 48 sec. 

71 min. 45 sec. 

[3.2 X 10“®,4,2 X 10“®] 

58 min. 57 sec. 

71 min. 18 sec. 

[4.2 X 10“®,5.0 X 10“®] 

59 min. 31 sec. 

73 min. 53 sec. 

[0.0,5.0 X 10“®] 

(Total) 363 min. 20 sec. 

(Total) 441 min. 45 sec. 


Table 11: Computation times of Drop in Computer Assisted Result [67^ 

It has taken about 13.4 hours for validating Drop in {e G [0,5.0 x 10“®]} in our computation 
environments. As for the slow shadowing and Jump, it totally takes less than 30 seconds to 
validate. 


In our topological theorems (Theorems |5.1[[A^ Corollaries 5.3 5.4 5.7 and 5.8), only the existence 
of desiring orbits is stated. As for the uniqueness of such orbits, we have no answers yet. More 
precisely, these theorems only require topological transversality in terms of covering relations, 
which is weaker than the transversality in terms of differential manifolds. If we can even verify 
transversality in the sense of differential manifolds, then Exchange Lemma gives us the local 
uniqueness of singularly perturbed orbits. Transversality is also essential to discuss stability of 


connecting orbits. Jones m discusses stability of homoclinic orbits for ( |6.1[ ) as the traveling wave 
solutions of FitzHugh-Nagumo PDE (6.2) from geometric viewpoints. For this stability arguments, 


some information about the nature of transversality between two manifolds play a central role. 

One of standard approaches to transversality of manifolds is to use Melnikov integrals via 
solving variational equations. If we take rigorous numerics into account, the C^-Lohner method 
|38| provides the variational information of trajectories as well as solutions themselves and will 
give an answer to solve transversality problems. 

Covering-Exchange for fast-slow systems with multi-dimensional slow variables 

Our validation arguments can be applied only to fast-slow systems with one-dimensional slow 
variable at present. It is natural to question whether our concept is applicable to fast-slow systems 
with multi-dimensional slow variables. The key problem concerning this extension is how we should 
track true trajectories near slow manifolds in terms of covering-exchange properties. Systems with 
multi-dimensional slow variables have rich behavior even on slow manifolds. The validation of slow 
manifolds by one fast-saddle-type block as well as concrete vector fields is thus more difficult than 
systems with one-dimensional slow variables because of nonlinearity of manifolds. It is not thus 


realistic to extend Theorem 5.1 and 5.6 to systems with multi-dimensional slow variable directly 
from the viewpoint of practical computations. 

By the way, existence theorems of trajectories with covering relations imply that true trajec¬ 
tories shadow reference orbits. In the case of singular perturbation problems, reference orbits 
correspond to singular orbits consisting of heteroclinic orbits and trajectories on slow manifolds. 


Slow shadowing discussed in Section [4~^ reflects the aspect of shadowing trajectories in fast-slow 
systems with one-dimensional slow variable. This expectation will be valid even for systems with 
multi-dimensional slow variables. Such an extended sequence will yield the existence of true tra¬ 
jectories which shadow those on slow manifolds. Higher dimensional extension of slow shadowing 
condition or its analogue will give an explicit criterion for validations with rigorous numerics in 
reasonable processes. 


Extension of procedures for slow manifolds with non-hyperbolic points 


Our validation arguments are based on the assumption that limiting critical manifolds are normally 
hyperbolic everywhere in consideration. Our procedure in this paper can be never applied if there 
is a point which is non-hyperbolic, like fold points, inside limiting critical manifolds. Nevertheless, 
there are a lot of cases that the interesting phenomena in singular perturbation problems are 
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involved by such non-hyperbolic points. One of famous examples is canard solutions^ which was 
discovered and first analyzed by Benoit, Callot, Diener and Diener [2]. A canard solution is a 
solution of a singular perturbed system which is contained in the intersection of an attracting slow 
manifold and a repelling one. Canard solutions provide a rich phenomenon such as canard explosion, 
namely, a transition from a small limit cycle to a relaxation oscillation through a sequence of canard 
cycles [5]. Since canard cycles involve fold points as jump points, our current implementations 
can never validate canard cycles. Krupa and Szmolyan |21j discuss the extension of geometric 
singular perturbation theory for slow manifolds including fold and canard points via the technique 
of blow-up. Simultaneously, Liu |22] discusses the Exchange Lemma in case that the critical 
manifold contains loss-of-stability turning points, corresponding the phenomenon called the delay 
of stability loss. In |22j , the Fenichel-type coordinate which overcomes non-normal hyperbolicity is 
constructed by the other type of invariant manifold theorem proved by Chow-Liu-Yi . We also 
mention sequential works by Schecter and Szmolyan [23 [23 [23 EO] , which discuss the singularly 
perturbed Riemann-Dafermos solutions as the small perturbation of composite waves consisting 
of constant waves, shock waves and rarefaction waves in systems of conservation laws. There, 
non-hyperbolic points called gain-of-stability turning points naturally arise in invariant manifolds. 
They derive the General Exchange Lemma for describing behavior of tracking manifolds near such 
non-hyperbolic manifolds with the help of blow-ups. 

We have to take such mathematical techniques into account, if we deal with invariant manifolds 
which are not normally hyperbolic with computer assistance. 
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A Concrete terms for (6.1) according to Section 2.3.3 


Here we show concrete forms of error terms and Jacobian matrices calculated from the vector 


field (6.1) for validating fast-saddle-type blocks, and cone conditions with the predictor-corrector 
approach (Section 2.3.3). 


A.l Fast-saddle-type blocks 

We list various terms in Section [2.3.3| in the case of FitzHugh-Nagumo system for readers’ accessi¬ 
bility. Let (u, 0, w) be a numerical zero of —f(u) + w, where f{u) = u(u — o)(l — u), a G (0,1/2). 
Use the coordinate {u, v, w) given by the following: 

u = u{ufu{u)~^w), V = V, w = w-\-'w, (A.l) 
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where (u,0,w) is a (numerical) root of {v = 0,cv — f(u) -I- ic = 0} and /«(«) is assumed to be 
invertible. The transformed vector field as the extended system is then 

u' = u' - fu{u)~^w' 

= v- fu{u)~^ec~^{u - 'yw) 

= v - ec~^fuiu)~^{u -I- ifu{u)~^ - 'y)w + u- ytc), 
v' = 6~^{cv — f{u) + w) 

= {cv — f {u + U + fu{u)~^w) + W + tc} , 
w' = €C~^{u + {fu{u)~^ — 7)w -I- It — 7w), 
r,' = 0, 


where we introduced an auxiliary variable 77 to be e = £ 977 , rj € [ 0 , 1 ], with constant eq. 
The fast component of the vector field is 



fv-tc V«(m) +(/«(«) ^ -'y)w + u -'yw)\ 

\ 6~'^ {cv - f {u + u +fuiuy^w)+w+ w} , j 




f-ec ^fuiu) ^{u+{fuiu) ^-'y)w + u- 
V {-fiu) + f{u) + fu{u)u + w} 



Denote (6.1) by 


X = F{X,Y,e), 


where X = {u,v)'^ and Y = w; namely, 


Y = eG{X,Y,e), 


F{X,Y,e) 



V \ 

- f{u)+'w)J ’ 


The Jacobian matrices at (u,0,w) are 


G(X,Y,e)=c ^(u — jw). 


Fx |(s.o.») > {Fx 

Fy l(B.o,®)= (/i) , {Fx^ o Fy) 1(^.0.®)= . 


f 6-^c 
\S-^fu{u) 



Around the numerical equilibrium {X,Y) = {{u,0),w) = {u,0,w), we introduce the affine trans¬ 
formation T : {Z, W) I—>■ (A, Y) of the predictor-corrector form by 

(A, Y) = T{Z, W) := {PZ + X- {F^^ o Fy) W,W + Y) . 

where P is the nonsingular matrix such that P~^Fx |(fl,o,is) P = A = diag(Ai, A 2 ). Then, denoting 

Z = ( 2 : 1 , 22 )^, 



Ai 2 i\ 

A 2 ^ 2 / 

Ai2iA 

A22:2 / 

Ai 2 i\ 

MZ 2 ) 




0 


0 

-5-^W 


0 


-(/(w) - fu{u){u -u)) + W 


-G{X,Y,e)/Uu) 

0 


-eGiX,Y,e)/fuiu) 
-d-^{f{u) - fuiu){u - u)) 


Note that the linear order term of IF-variable in the error term is eliminated. 

Let A C be an h-set (in the (it, v, it;)-coordinate) containing {u, 0, w). Then the Mean Value 
Theorem implies that the higher order term f{u) — /„('u)(it — u) is included in the enclosure 


nu) 


(u — uY I It, u such that (u, u, ic), (it, v,w) G N 
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Thanks to the term (u — rt)^, this set is in general very small if N is small. In our example, 
f''{u) = —6u-|-2(a-|-1). Obviously, the first component —€G{X,Y,e)/fu{u) is very small if e > 0 is 
sufficiently small. Note that the denominator fu{u) is bounded away from 0 since we are focusing 
on validation of normally hyperbolic invariant manifolds. Finally, for given compact set N and 
eo > 0, we obtain enclosures of the error term 

p-i f fu{u) \ iu,v,w) G N, e e [0,eo]] \ ^ 

y [(5“^{3 m — (a -I- 1 )}(w — \ u,u such that {u, v, w), {u, v,w) G N^J ’ 

which can be computed by interval arithmetics and, in principle, returns very small enclosures for 
small N and cq. 


A.2 Jacobian matrices for cone conditions 

The Jacobian matrix in this coordinate at {u,v,w,r]) ignoring the differential of ry'-term is 
^-e(c/,i(u))“^ 1 e(c/„(u))“i(-/,,(w)“i + 7) -eo(c/„(?2))“i(M-f (/„(u)“i - 7)w-b u - yw) 

-S~^fuiu) <5“^c {Sfu{u))~^ {-fuiu) + fu{u)} 0 


€C 


-1 


0 


ec - 7 ) 


eoc '^(u+ifuiu) ^-j)w + u-jw) 

(A.2) 


where the variable u in the second row corresponds to (A.l). We additionally transform {u^v) 
linearly to (zi, Z 2 ) by (m, v) = P(zi, Z 2 ), where P is the nonsingular matrix such that 


p-i 


-S ^fu(,u) 




The Jacobian matrix (A.2) in the new coordinate ( 21 ,^ 21 ^^) is then written as 


/Ai 0 0 0\ 

0 A 2 0 0 

\ 0 0 0 0 / 

/ Fii Fi 2 e{cfuiu))~'^{-fuiu)~'^+-f) -eo{cfuiu))~^iu+{fu{u)~'^-j)w + u-jw)\ 

+ A21 F22 {Sfu{u))~^ {-fu{u) + fuiu)} 0 

\ec“^ 0 - 7) eoc~'^{u +{fu{u)~^ - j)w+ u-jw) ) 

(A.3) 

(Fii F^ 2 \ _ p_i ( -e(c/„(h))-i 0\ 

'v ^’21 F22) \S-\fu{u) - fu{u)) O)^- 


In order to verify cone conditions in a given h-set N containing (u, 0, w), we compute corresponding 
maximal singular values in N. All such calculations are done by interval arithmetics. As for 
calculations of the enclosure of fu{u) — fu{u) in N, it is useful to use the mean value form 


fu{u) - fu{u) G {-fuu{u){u - u) \u,u such that {u,v,w), {u,v,w) G N} 

= {(6u^ — 2(a + l))(u — u) \ u,u such that {u, v, w), {u, v, w) G N}. 

We apply the matrix ( |A.3 ) as the Frechet differential of the vector field in ( |3.1[ ) to validating 
cone conditions discussed in Section[3 The m-cone conditions can be treated in the similar manner. 
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